Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,145 G→C 16.2% intergenic (‑196/+33) ibsC ← / ← serA toxic membrane protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1450GC16.2% 645.2 / 109.2 236intergenic (‑196/+33)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (87/110);  new base C (16/22);  total (103/133)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 238 total aligned reads displayed.

CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCA  >  NC_000913/3057046‑3057244
                                                                                                   |                                                                                                   
  ttGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                                                                                                   <  2:2116484/100‑1 (MQ=255)
    gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                                 <  1:3213824/100‑1 (MQ=255)
      tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGActt                                                                                               >  2:922075/1‑98 (MQ=255)
       cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCATCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTg                                                                                              >  2:1788396/1‑100 (MQ=255)
       cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTg                                                                                              >  1:1919826/1‑100 (MQ=255)
         gACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                            <  2:3771676/100‑1 (MQ=255)
         gACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGCGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                            >  1:2223199/1‑100 (MQ=255)
          aCTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                                                                                           <  1:3843191/100‑1 (MQ=255)
          aCTTTGAGTTATTTTCTGGAGCAAGCTCGCAAAGTAGAATGCGCAACGCGGCAACGTTGTGGAGAAAGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                <  2:3362295/95‑1 (MQ=255)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                      <  2:2620218/100‑1 (MQ=255)
                atttatttTCTGGAGCAGAATCGCAAAGTAGAACGCTCAACGCGGAAACGATGTGGAGACGTGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                      <  2:3944820/99‑1 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                    <  1:343923/100‑5 (MQ=38)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                    <  2:3458708/100‑5 (MQ=38)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                    <  1:922075/100‑5 (MQ=38)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                    <  2:904040/100‑5 (MQ=38)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                    >  2:1581412/1‑100 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                    >  2:1373800/1‑100 (MQ=255)
                  ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                                   <  1:3385931/100‑1 (MQ=255)
                    atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttttat                                                                                 <  1:1620329/100‑8 (MQ=37)
                    atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                 >  1:991209/1‑100 (MQ=255)
                    atttTCTGGAGCAGACTCGCAAAGTAGAAAGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                 <  2:206171/100‑1 (MQ=255)
                     ttttCTGGAGCATACTCCCGAAGTAGGACGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAGACACGGGCAGGGCACTGACTTGCCCGtttttttatc                                                                                <  2:3620261/100‑9 (MQ=16)
                     ttttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttttatc                                                                                <  1:300254/100‑9 (MQ=37)
                     ttttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                <  2:2149298/100‑1 (MQ=255)
                      tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               <  2:1119425/100‑1 (MQ=255)
                      tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               <  1:849320/100‑1 (MQ=255)
                      tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               <  1:759625/100‑1 (MQ=255)
                      tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               >  1:2659250/1‑100 (MQ=255)
                        tCTGGAGCAGCCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGAGGAGAAGGGATAAAAAAACGGGTAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                             >  2:4023790/1‑100 (MQ=255)
                         cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                            >  1:181281/1‑100 (MQ=255)
                            gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                         <  1:195026/100‑1 (MQ=255)
                              gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                       <  1:3175469/100‑1 (MQ=255)
                              gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                       >  2:1655007/1‑100 (MQ=255)
                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                   <  1:967809/100‑1 (MQ=255)
                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                   <  2:646639/100‑1 (MQ=255)
                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                                 <  1:3128365/100‑1 (MQ=255)
                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                                 >  2:1203064/1‑100 (MQ=255)
                                     cGCAAAGTAGAATGCTCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                <  1:4123917/100‑1 (MQ=255)
                                       cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                               <  1:400877/99‑1 (MQ=255)
                                       cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                              <  1:2078579/100‑1 (MQ=255)
                                         aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                            <  1:125288/100‑1 (MQ=255)
                                         aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                            <  1:1604059/100‑1 (MQ=255)
                                           gTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGACTTAGTAcagca                                                          >  1:2922144/1‑100 (MQ=255)
                                           gTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagca                                                          >  1:1044076/1‑100 (MQ=255)
                                              gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagcaact                                                       >  1:3349941/1‑97 (MQ=255)
                                               aaTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACg                                                      <  1:3550606/100‑1 (MQ=255)
                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                  >  2:3853962/1‑100 (MQ=255)
                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                  >  1:666350/1‑100 (MQ=255)
                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                  >  1:3901111/1‑100 (MQ=255)
                                                          cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATgn                                           <  2:1974930/100‑2 (MQ=255)
                                                          cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATg                                            >  2:2350239/1‑99 (MQ=255)
                                                                 ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                                    >  2:3236873/1‑100 (MQ=255)
                                                                   tgtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGaa                                  >  1:3732024/1‑100 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                 <  1:1751802/100‑1 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                 <  1:2876246/100‑1 (MQ=255)
                                                                      ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                               >  1:3784791/1‑100 (MQ=255)
                                                                       gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                              >  2:749399/1‑100 (MQ=255)
                                                                       gagaAGGCATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                              <  1:644867/100‑1 (MQ=255)
                                                                        agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                             >  1:2707141/1‑100 (MQ=255)
                                                                        agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                             >  1:2136736/1‑100 (MQ=255)
                                                                         gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCtt                            <  1:3561485/100‑1 (MQ=255)
                                                                          aaGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCttt                           <  2:3085625/100‑1 (MQ=255)
                                                                           aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                          <  1:3202278/100‑1 (MQ=255)
                                                                           aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                          >  1:1422343/1‑100 (MQ=255)
                                                                           aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                          <  2:3200919/100‑1 (MQ=255)
                                                                           aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                          <  1:550970/100‑1 (MQ=255)
                                                                               aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTg                      <  2:18638/100‑1 (MQ=255)
                                                                                tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                     >  2:1578199/1‑100 (MQ=255)
                                                                                 aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                    <  1:408961/100‑1 (MQ=255)
                                                                                 aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                    >  1:51638/1‑100 (MQ=255)
                                                                                  aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                   >  1:3377696/1‑100 (MQ=255)
                                                                                    tcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                 >  1:2629730/3‑100 (MQ=38)
                                                                                    aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                 <  1:150549/100‑1 (MQ=255)
                                                                                    aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                 <  2:306425/100‑1 (MQ=255)
                                                                                    aaaaCGGGCAAGTCAGGTACCTGCCCGTTGATTTTCAGAGAAGTGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                 <  2:2192495/100‑1 (MQ=255)
                                                                                      aaCGGGCAGGGCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAGGGTACCCGGAATAGCCTTCATTGCCTGCAg               >  2:832309/1‑100 (MQ=255)
                                                                                      aaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg               <  1:2163310/100‑1 (MQ=255)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  1:4065118/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  1:114866/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:2867950/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  1:1000886/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCCTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:47980/1‑100 (MQ=255)
                                                                                       aCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              <  1:1314333/100‑1 (MQ=255)
                                                                                        cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:459307/1‑99 (MQ=39)
                                                                                         gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc            >  1:441893/1‑100 (MQ=255)
                                                                                         gggCAAGTCAGTGACCTGCCCGTTGAGTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc            <  2:3642211/100‑1 (MQ=255)
                                                                                          ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt           <  2:2951229/100‑1 (MQ=255)
                                                                                          ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt           >  2:3412722/1‑100 (MQ=255)
                                                                                          ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt           <  2:2707134/100‑1 (MQ=255)
                                                                                           gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt          >  1:3635142/1‑100 (MQ=255)
                                                                                           gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt          <  1:3811052/100‑1 (MQ=255)
                                                                                           gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt          <  2:1422340/100‑1 (MQ=255)
                                                                                              aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt       >  2:4083061/1‑100 (MQ=255)
                                                                                                taaGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg     <  1:2163787/98‑1 (MQ=255)
                                                                                                tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg     <  2:1069125/100‑1 (MQ=255)
                                                                                                tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg     <  2:1770692/100‑1 (MQ=255)
                                                                                                 cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg    <  1:1544527/100‑1 (MQ=255)
                                                                                                 cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg    <  2:3187167/100‑1 (MQ=255)
                                                                                                  aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGc   >  1:859881/1‑100 (MQ=255)
                                                                                                   gTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCa  >  2:3818345/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCA  >  NC_000913/3057046‑3057244

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: