Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,757 T→C 15.9% intergenic (‑140/+24) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7570TC15.9% 351.5 / 56.0 153intergenic (‑140/+24)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (66/62);  new base C (15/9);  total (82/71)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 158 total aligned reads displayed.

GATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478658‑478856
                                                                                                   |                                                                                                   
gATTGTGGAAAAAGGCTGGCTGGAGATAAGGGGGTAGGGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGGAGAGCtt                                                                                                     >  1:3615330/1‑100 (MQ=255)
    gtgGAAACAGCCTGGATGGAAATAAGTAGGTAGATCTGGTATTTCATGTAGTGATGGGTGCTAACGATCATTGTCATTTTACTAGTCGCAGAGCTTTACg                                                                                                 <  1:1315206/100‑1 (MQ=255)
    gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                 >  2:2090600/1‑100 (MQ=255)
    gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                 >  1:3577715/1‑100 (MQ=255)
    gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                 >  2:1556207/1‑100 (MQ=255)
    gtgGAAAAAGGCTGGATGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                 <  2:3759282/100‑1 (MQ=255)
      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               >  1:1765879/1‑100 (MQ=255)
      ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGATGGTATTTCATGAAGTGATGGGTGGAAACGATCATTCTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                               >  1:2299506/1‑100 (MQ=255)
       gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTGTTTCATGAAGTGATGGGTGGAAACGATGATTGTGATTTGACAAGGCGCCGAGATTTACGCCt                                                                                              >  1:3672410/1‑100 (MQ=255)
       gAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCt                                                                                              >  2:1354132/1‑100 (MQ=255)
       gAAAAAGACTGGCTGGAGATAAGGAGGTCGAGCTGGTATTTAATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                 >  1:1309714/1‑97 (MQ=255)
        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             >  2:498011/1‑100 (MQ=255)
        aaaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTg                                                                                             >  2:1666150/1‑100 (MQ=255)
         aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                            <  2:410118/100‑1 (MQ=255)
          aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                           <  2:1361773/100‑1 (MQ=38)
           aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                          <  2:618487/100‑2 (MQ=38)
           aaGGCTGGCTAGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGAGCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                          <  1:2531970/100‑2 (MQ=37)
            aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                         <  2:1633150/100‑3 (MQ=38)
            aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                         <  2:1980857/100‑3 (MQ=38)
            aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                         <  2:648653/100‑1 (MQ=255)
            aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                         <  2:918174/100‑1 (MQ=255)
                 ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                    <  2:2970398/100‑1 (MQ=255)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                   >  2:2257961/1‑100 (MQ=255)
                  gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                    <  2:3064257/99‑1 (MQ=255)
                   cttccgatctggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                  <  1:1189073/90‑1 (MQ=255)
                   ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                  >  1:176773/1‑100 (MQ=255)
                   ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                  <  2:2047040/100‑1 (MQ=255)
                   ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                  >  2:2458541/1‑100 (MQ=255)
                      gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                               >  1:1001193/1‑100 (MQ=255)
                      gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTaa                                                                               >  2:2479682/1‑100 (MQ=255)
                           aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                          <  2:3352817/100‑1 (MQ=255)
                             ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcccagatcggaag                                                                        >  2:1189080/1‑90 (MQ=255)
                                 gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                    <  2:3082192/100‑1 (MQ=255)
                                 gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                    <  2:1236163/100‑1 (MQ=255)
                                 gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAAAGATCATTGTCATTTGAGCAGGCGCAGAGGTTTACGCCTGGTATGTAACCCTAACTCAAGATGt                                                                     >  1:2406054/1‑99 (MQ=255)
                                  tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                   <  2:3615352/100‑1 (MQ=255)
                                   agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                                  <  1:2179294/100‑1 (MQ=255)
                                       cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGTGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAt                                                                 <  1:980315/97‑1 (MQ=255)
                                         ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                            >  1:3192929/1‑100 (MQ=255)
                                         ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                            >  1:359890/1‑100 (MQ=255)
                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         >  1:666263/1‑100 (MQ=255)
                                            aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGcgcg                                                         <  1:784490/100‑1 (MQ=255)
                                              ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                       >  1:604468/1‑100 (MQ=255)
                                                   gAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                  >  1:2943510/1‑100 (MQ=255)
                                                    aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                 <  1:1502932/100‑1 (MQ=255)
                                                     aGTGATGGGTGGAAACTATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                <  2:315268/100‑1 (MQ=255)
                                                      gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                               >  1:927151/1‑100 (MQ=255)
                                                       tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaa                                              >  2:1988870/1‑100 (MQ=255)
                                                          tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                           <  1:1695461/100‑1 (MQ=255)
                                                          tGGGTGGAAACGATCATTGTCATTTGACCAGGCGAAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                           <  1:3787215/100‑1 (MQ=255)
                                                             gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  2:2888230/1‑100 (MQ=255)
                                                             gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                        >  1:1916377/1‑100 (MQ=255)
                                                              tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                       >  1:1131666/1‑100 (MQ=255)
                                                               ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                      <  2:1097257/100‑1 (MQ=255)
                                                               ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                      <  2:3616916/100‑1 (MQ=255)
                                                               ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                      <  2:2982289/100‑1 (MQ=255)
                                                                    cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                 >  2:43564/1‑100 (MQ=255)
                                                                     gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                 >  1:2810666/1‑99 (MQ=255)
                                                                      aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                               <  1:944710/100‑1 (MQ=255)
                                                                           tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                          >  2:2954120/1‑100 (MQ=255)
                                                                           gggttacataCCAGGCGTAAAGCTCTGCGCCTGGTAAGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                          >  2:2689138/10‑100 (MQ=25)
                                                                            ggttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGAGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                         >  1:2197932/9‑100 (MQ=25)
                                                                            ggttacataCCAGGCGTAAAGCTCTGCGCCTGGGATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGGGAGGGTAAAAATATCGTAAAAAGACAGtt                         >  1:2292040/9‑100 (MQ=25)
                                                                            gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                         >  1:1273398/1‑100 (MQ=255)
                                                                            gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                         >  1:2093034/1‑100 (MQ=255)
                                                                                  tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                   >  1:1839074/1‑100 (MQ=255)
                                                                                  tGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                   <  1:1744055/100‑1 (MQ=255)
                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  2:1734169/2‑100 (MQ=38)
                                                                                   taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  1:2844801/2‑100 (MQ=38)
                                                                                   taCCAGGCGTAAAGCTATGCGCTTGGAATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAGAAGACAGtcaatcttt                  >  1:3700881/2‑92 (MQ=255)
                                                                                   gACCAGGAGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGGGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  1:2005808/1‑100 (MQ=255)
                                                                                    aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                  >  1:2840669/1‑99 (MQ=38)
                                                                                    aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                 >  2:3648659/1‑100 (MQ=255)
                                                                                     ccAGGCGTAAAGCTCTGTGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGATCACCTTTAc                >  1:1939165/1‑100 (MQ=255)
                                                                                     ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  1:1609765/1‑100 (MQ=255)
                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  2:3355594/1‑100 (MQ=255)
                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                >  2:208016/1‑100 (MQ=255)
                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                <  2:1928356/100‑1 (MQ=255)
                                                                                     ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                <  1:3798760/100‑1 (MQ=255)
                                                                                     ccAGACGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAGTATCGTAAAACGACAGTTCACCttt                  >  1:2880789/1‑98 (MQ=255)
                                                                                      cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATTTTGATTCGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGACAGTTCACCTTTACg               >  2:1254592/1‑100 (MQ=255)
                                                                                      cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg               <  1:2791363/100‑1 (MQ=255)
                                                                                      cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg               >  2:3205055/1‑100 (MQ=255)
                                                                                        ggCGTAGAGCTTTACGCCTGGTATGGAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg               <  1:1667080/98‑1 (MQ=255)
                                                                                        ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa             >  2:1405977/1‑100 (MQ=255)
                                                                                        ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa             >  2:832959/1‑100 (MQ=255)
                                                                                          cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           >  2:1909528/1‑100 (MQ=255)
                                                                                           gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa          >  2:3126269/1‑100 (MQ=255)
                                                                                           gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa          <  2:3466975/100‑1 (MQ=255)
                                                                                           gcAGAGCTTTACCCCTGGTATGTAACCCTAACGTAAGATGTTCATTCGGGTCGAAAAAAAAAAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg           <  1:2737813/99‑1 (MQ=255)
                                                                                            cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         <  1:2985429/100‑1 (MQ=255)
                                                                                            cAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         <  2:2299517/100‑1 (MQ=255)
                                                                                             agagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         >  2:928811/1‑99 (MQ=255)
                                                                                              gtgCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGATCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt         <  1:2109052/96‑1 (MQ=255)
                                                                                              gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa       >  2:1972977/1‑100 (MQ=255)
                                                                                               agCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAt      <  2:1905603/100‑1 (MQ=255)
                                                                                                gCTTTACGCCTGGTATGTAAACCTAACTTCAGATGTTGATTAGGCGCGAGCGAGGGTAAAACTATCGTAAAAAGAAAGTTCACCTTTACggcag           >  1:3056082/1‑91 (MQ=255)
                                                                                                  tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt   <  2:1803763/100‑1 (MQ=255)
                                                                                                  tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt   >  2:1704294/1‑100 (MQ=255)
                                                                                                   ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt  <  2:1068855/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTT  >  NC_000913/478658‑478856

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: