D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	61,065	G→A	G48G (GGC→GGT)	SAUSA300_RS00265 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	61,065	0	G	A	99.2%	367.4 / NA	127	G48G (GGC→GGT)	SAUSA300_RS00265	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (7/119); total (7/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

ATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATAA  >  NC_007793/61000‑61160

tgctcttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTCTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATcc                                                               <  1:1447598/90‑1 (MQ=255)

 tcttccgatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCa                                                              <  1:2455282/90‑1 (MQ=255)

 gctcatccgaTCTGCCTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCa                                                              <  1:220221/91‑1 (MQ=255)

   tcttccgaTCTGCCTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATa                                                            <  1:1142155/93‑1 (MQ=255)

    cttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATAc                                                           <  1:3360591/94‑1 (MQ=255)

     ttccggtCTGCCTATACAAAAACTAATCCATGTCAGAATATGATTTCCTAAAGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa                                                          <  1:4409056/94‑1 (MQ=255)

     ttccgaTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa                                                          <  1:5224148/95‑1 (MQ=255)

     ccgatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACa                                                          <  1:2334305/94‑1 (MQ=255)

       gatctCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACaaa                                                        <  1:2776804/96‑1 (MQ=255)

          aTCTGCCTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATaa                                                     <  1:3021366/100‑1 (MQ=255)

              gcctgcacaaaaACTAATCCATGTCATAATATGATTCCCTAATGCCTTTTTACCACGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 >  1:312784/9‑100 (MQ=255)

              gccTATACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 <  1:49580/99‑1 (MQ=255)

              gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 >  1:4033749/2‑100 (MQ=255)

              gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 >  1:5478884/2‑100 (MQ=255)

              gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 >  1:3628319/2‑100 (MQ=255)

              gccTACACAAAAACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGt                                                 <  1:4712518/99‑1 (MQ=255)

              gccTACACAAAAACTAATCCATGTCATAATATGAGTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGGTCATCAGCAAGTATCCATACAAATAAAGGt                                                 <  1:3935916/99‑1 (MQ=255)

               cctccacaaaaACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTa                                                <  1:5467683/93‑1 (MQ=255)

               cctccacaaaaACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTa                                                <  1:2432831/93‑1 (MQ=255)

                cggcacaaaaACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGCCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTaa                                               <  1:3248254/94‑1 (MQ=255)

                  ataCAAAAACTAATCCATGTCATAATATGATTCCCTAAGGCGTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTagca                                             <  1:2640796/100‑4 (MQ=255)

                      aaaaACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGAc                                         <  1:5555394/100‑1 (MQ=255)

                      aaaaACTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGAc                                         <  1:4624969/100‑1 (MQ=255)

                        aaaCTAATCCATGTCATAGTATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCTAGTATCCATACAAATAAAGGTAAAATGACAg                                       <  1:5038059/98‑1 (MQ=255)

                          aCTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGt                                     <  1:249309/100‑1 (MQ=255)

                           cTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGta                                    <  1:5611630/100‑2 (MQ=255)

                           cTAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTg                                    <  1:4085334/100‑1 (MQ=255)

                            tAATCCATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGc                                   <  1:377691/100‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATGCAAATAAAGGTAAAATGACAGGTGCg                                  <  1:2132753/97‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGATAGGTGCg                                  >  1:2692673/4‑100 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGcgcg                                  >  1:1708366/4‑100 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGct                                  <  1:420304/97‑2 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCg                                  <  1:5503437/97‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCg                                  <  1:4201379/97‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCg                                  <  1:4386860/97‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCg                                  <  1:4207953/97‑1 (MQ=255)

                             aatccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCg                                  <  1:3301640/97‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:4822449/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:1385751/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:5436286/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:5034526/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:2341241/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:2345561/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGa                                 <  1:4870723/98‑1 (MQ=255)

                              atccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCAGCAGCAAGTATCCATACAAATAAAGGTAAGATGACAGGTGCGa                                 <  1:2585414/98‑1 (MQ=255)

                               tccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaa                                <  1:916508/99‑1 (MQ=255)

                               tccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaa                                <  1:3704873/99‑1 (MQ=255)

                               tccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaa                                <  1:4538667/99‑1 (MQ=255)

                               tccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaa                                <  1:620046/99‑1 (MQ=255)

                               tccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaa                                <  1:2565557/99‑1 (MQ=255)

                                ccATGTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaa                               <  1:5650104/100‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:740667/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:3515593/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:1214596/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:3450811/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:1860485/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:3409549/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:1993825/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:4337639/97‑1 (MQ=255)

                                 catgTCATAATATGATTCCCTAATGCTTTTTTACCATGTCCTGATGTTGGCTCACCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGaaaa                              <  1:100971/97‑1 (MQ=255)

                                        taataTGATTCCCTAATGC‑TTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACtt                      <  1:15794/100‑1 (MQ=255)

                                                   cctaatgCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTACGAAAAAGATACTTAAATAACTCa            <  1:2923527/97‑1 (MQ=255)

                                                       atgCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGc        <  1:4804277/100‑1 (MQ=255)

                                                        tgCTTTTTTGCCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCt       <  1:3391975/100‑1 (MQ=255)

                                                        tgCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCt       <  1:2621623/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCta      <  1:2012041/100‑2 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:5403035/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:5089558/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:4960887/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:3639657/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:878671/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:2502986/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:4053516/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:3788418/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:479860/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:2900481/100‑1 (MQ=255)

                                                         gCTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTg      <  1:3923573/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTgt     <  1:5124025/100‑2 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTgg     <  1:2806003/100‑2 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:59921/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:1437549/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:517184/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:2852262/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:3777347/100‑1 (MQ=255)

                                                          cTTTTTTACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGa     <  1:4661375/100‑1 (MQ=255)

                                                          cTTTTTTACCAAGTCCTGATGTTGGCTCATCAGCAAGTACCCATACAAATAAAGGTAAAATGACAGGTACGAAAAAGATACTTAAATAACTCAATGCTgg     <  1:873352/100‑2 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGag    <  1:650317/100‑2 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGAt    <  1:1071903/100‑1 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGAt    <  1:1988241/100‑1 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGAt    <  1:5372861/100‑1 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGAt    <  1:4383470/100‑1 (MQ=255)

                                                           ttttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGAt    <  1:252052/100‑1 (MQ=255)

                                                           ttttttACCATGCCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTgcc    <  1:3981727/100‑3 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:3589640/100‑1 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:4007409/100‑1 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:179022/100‑1 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:1952097/100‑1 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:3040482/100‑1 (MQ=255)

                                                            tttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATa   <  1:2472069/100‑1 (MQ=255)

                                                             ttttACCATGTCCTGATGTTGGCTCATCAGCAAGTATCCGTACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATaa  <  1:4937263/100‑1 (MQ=255)

ATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAATGACAGGTGCGAAAAAGATACTTAAATAACTCAATGCTGATAA  >  NC_007793/61000‑61160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`