Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 12.1% intergenic (+64/‑335) alx → / → sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6090TG12.1% 369.0 / 23.0 149intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (83/37);  new base G (2/14);  total (85/64)
Fisher's exact test for biased strand distribution p-value = 1.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239462‑3239758
                                                                                                                                                     |                                                                                                                                                     
 gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:266083/150‑1 (MQ=255)
 gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:153932/150‑1 (MQ=255)
 gCGTGGTGTTTGGCATTCTGGTGATGACGATTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:2626902/150‑1 (MQ=255)
    tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 >  1:3052053/1‑150 (MQ=255)
    tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 >  1:3773513/1‑150 (MQ=255)
    tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 >  1:811123/1‑150 (MQ=255)
        gtgtGGCATTCTGGTGATGACGTGTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCTCTCCCTGGAAGGGagag                                                                                                                                             <  2:2497902/147‑1 (MQ=255)
        gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                             <  2:173269/150‑1 (MQ=255)
        gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                             <  1:2363584/150‑1 (MQ=255)
        gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                             >  1:1774553/1‑150 (MQ=255)
        gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:51207/150‑1 (MQ=255)
         tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:3104831/1‑150 (MQ=255)
         tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:2766271/1‑150 (MQ=255)
          ttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGTGAggg                                                                                                                                           <  1:157686/150‑1 (MQ=255)
           tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAggga                                                                                                                                          <  2:268094/150‑2 (MQ=255)
             gCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgggaga                                                                                                                                        <  2:569737/150‑4 (MQ=255)
              cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                       <  2:1974312/150‑1 (MQ=255)
              cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                       <  1:2255001/150‑1 (MQ=255)
                ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:2977159/1‑150 (MQ=255)
                   tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  >  1:2965057/1‑150 (MQ=255)
                    ggTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                                                                                                                 >  2:498854/1‑150 (MQ=255)
                    ggTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAg                                                                                                                                 <  2:431249/150‑1 (MQ=255)
                       gttgatttttATTATCAACGCCTGGGTGATTTTTCGGCATGATATGCAGCGGGGTGGATATTTTTTAATCTGCCTACGCCGTGTACCCTGTCTTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTGGGGGTTAGGGt                                                                                                                              <  2:2871900/143‑1 (MQ=255)
                         tgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTaa                                                                                                                            >  2:3549064/1‑150 (MQ=255)
                            cGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                                                                                                         >  1:2003200/1‑150 (MQ=255)
                             gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  1:1210665/1‑150 (MQ=255)
                             gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtttt                                                                                                                        >  1:1999870/1‑150 (MQ=255)
                              tttATTATCATCGCCTGGGTGTATTTTGGGCATGTTAAGCAGCGGGGTGGATTATTTTTAATCTGCCTAAGCGTTGTACCCTTTCATTAACTTGAGCACCGTTTTCTCACTCTCCCTTCCAGGTAGTGGGTCGGGGTGAGGGTAAttttt                                                                                                                       <  2:1971953/150‑1 (MQ=255)
                                    atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                 >  2:3925339/1‑150 (MQ=255)
                                     tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                                >  1:2184858/1‑150 (MQ=255)
                                     tCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                                                                                                >  1:2184842/1‑150 (MQ=255)
                                             tGGGTTAATTTTGGGATTTATACGCGGGGGGGTGTATATTTTTTAATGTGCATACGCGGTGTACCCTGTAATTAACTTGAGCACGGTTTTCTCCCGCCCCCTTCCGGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGc                                                                                                        <  2:3354589/150‑1 (MQ=255)
                                                 ttatttATGGGAATGATAAGCTGGGGGGTTTAAATTTTTTAATCTGCATAAGCGTTGTACCCTTTAATTAACTTCACCACCTTTTCCCCCCTCACCCTCCCGGGGTGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                                                                                    <  1:2871576/146‑1 (MQ=255)
                                                 tGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                                                                                    >  2:756128/1‑150 (MQ=255)
                                                          ggCATTATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTATGCCGTGTACCCTGTCATTAACATGAGCAGCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           <  2:2645139/150‑1 (MQ=255)
                                                           gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  1:1920543/1‑150 (MQ=255)
                                                           gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:1148521/1‑150 (MQ=255)
                                                           gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCAcc                                                                                          >  2:1148451/1‑150 (MQ=255)
                                                                    gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctc                                                                                 >  1:2524731/1‑150 (MQ=255)
                                                                    gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCGAAcccac                                                                                 >  1:2334125/1‑148 (MQ=255)
                                                                      cggggggTGTATAATTTTTAATCTGCCTAAGCGGTTTACCCTTTATTTAACTTGAGCACGGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTGGGTGTGAGGGTAAATTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2434783/147‑1 (MQ=21)
                                                                      cggggggTGGATAATTTTTAATCTGCCTATGCGGTGTACCGTGTCTTTAACTTGTGCACGGTTTTCTCCCTCTCCCTTCCCGGGTGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:3841792/147‑1 (MQ=25)
                                                                      agcGGGGTGTATAATTTTTAATCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCTTTTTCTCCCTCTCCCTCCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  1:2789954/150‑1 (MQ=38)
                                                                       ggggggTGTATAATTTTTAATCTGCCTAAGCGGTGTACCCTTTAATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                              <  1:2827778/148‑1 (MQ=37)
                                                                       cgggggTGGTTAATTTTTTAATCTCCAAAGGCGGGTTACCTTTTATTTAACTTGAGCCCCTTTTACTCCAGCCCCCTGCCGGGTTGTGGGTGGGGGTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2354610/148‑1 (MQ=255)
                                                                             tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        >  1:3280626/1‑150 (MQ=255)
                                                                             tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        >  1:3280643/1‑150 (MQ=255)
                                                                             tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        >  1:2859743/1‑150 (MQ=255)
                                                                               gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  2:3889112/1‑150 (MQ=255)
                                                                                  aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   >  1:2625126/1‑150 (MQ=255)
                                                                                  aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGc                                                                   >  1:3383338/1‑150 (MQ=255)
                                                                                   tttttttATTCTGCCTACGCGTTGTACCGTGAAATTAACATGCGCACGGTTTTCTCCCTCACCCTTCCCGGTAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   <  2:3147424/148‑1 (MQ=21)
                                                                                       tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACTTGCGCACCGTTTTCGCCCTCTCCATTCCCGGGTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  2:2953048/150‑1 (MQ=37)
                                                                                       tttAATCGGCCTACGCCGTGTCCCCGGTCATTAACAGGCGCACCGTTTTCTCCCTCGCCCTTCCTGGGTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  1:2426596/150‑1 (MQ=25)
                                                                                        ttACTCGGCCGACGCGGTGTACCCTGCCTTTAACTGGGGCCCGGTTTTCTCCCTCGCCGTGCCGGGTTGCGGGGGGGGGGGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                              <  2:3510822/150‑1 (MQ=12)
                                                                                        ttAATCTGCATACGCCGTGTACCCTGTAATTACCAGTCGAACGTTTTTCCCCCTCTCCCGCCCTGGTTGTGGGTGGGGGTTTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                              <  2:1561297/150‑1 (MQ=16)
                                                                                            tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGatt                                                          >  2:774623/1‑148 (MQ=255)
                                                                                               gCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAc                                                       >  2:1272470/1‑150 (MQ=255)
                                                                                                ccaacgCCGTGTACCTTGAAATTAACTGGAGCACGTTTTTCACCCCCCCCCTTCCCGGGAGTGGTGGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      <  1:3494061/145‑1 (MQ=18)
                                                                                                ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      >  2:1759379/1‑150 (MQ=255)
                                                                                                ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTCCCCGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      <  1:1644147/150‑1 (MQ=37)
                                                                                                      ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTcc                                                >  1:126876/1‑150 (MQ=255)
                                                                                                       ccggggaCCCTGTCAGTACCATGCGCACCGTTTGCTCACTCCCCCGACCCGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCACCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTcc                                                <  2:1609020/144‑1 (MQ=21)
                                                                                                        gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtg                                              >  2:1835484/1‑149 (MQ=255)
                                                                                                         ggtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttt                                             <  2:631979/149‑1 (MQ=255)
                                                                                                             aCCCTTTCATTAACTTGAGCACCTTTTCCTCCCTCGCCCTTCCCGGGTGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:2004407/150‑1 (MQ=25)
                                                                                                             aCCCTGTCATTAACATGCGCACCGTTTTCTCCCTCTCCCTGCCAGGGAGAGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  2:3822967/150‑1 (MQ=255)
                                                                                                              cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTg                                         >  1:2006856/1‑150 (MQ=255)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                        <  1:3602416/150‑1 (MQ=255)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACATACCGAGTCCtgtttttt                                        >  1:1754352/1‑148 (MQ=255)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTACCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTgt                                        >  1:1019435/1‑149 (MQ=255)
                                                                                                               ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCACACCCAAACCCTCTCCCCAAACGGGGCGAGGGCACTGACCGAGTCCTTTTTTGa                                        >  1:2310624/1‑150 (MQ=255)
                                                                                                                    aaTTAACATGAGCACCTTTTACCCCCCCCCCCCGCCCGGGTGTGGGGGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttg                                   <  1:3182307/149‑1 (MQ=21)
                                                                                                                           aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                            >  2:1617895/1‑150 (MQ=255)
                                                                                                                             gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTc                          >  2:2512244/1‑150 (MQ=255)
                                                                                                                              cgCACGGTTTTCACCCCCTCCCTGCCCGGGTGTGGGTGGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                         <  1:100371/149‑1 (MQ=25)
                                                                                                                               gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTGATGTTGTCATCAGTCTg                        >  2:2095527/1‑150 (MQ=255)
                                                                                                                                cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTGATGTTGTCATCAGTCTgg                       >  1:2955607/1‑150 (MQ=255)
                                                                                                                                cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTgg                       >  2:3325512/1‑150 (MQ=255)
                                                                                                                                cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTgg                       >  2:3626937/1‑150 (MQ=255)
                                                                                                                                 aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATAAGTCTGGa                      >  2:942355/1‑150 (MQ=255)
                                                                                                                                      tttCTCCCTCTCCGTGCCCGGTTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCg                 <  2:2737757/150‑1 (MQ=37)
                                                                                                                                        tctccctcacccTGCCCGGGAGTGGGTGGGGGTGAGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGAACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCa               <  2:1248259/150‑1 (MQ=37)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTGGTTGTCATCAGTCTGGAAGCCGCACg             >  1:1239594/1‑150 (MQ=255)
                                                                                                                                            cctctcccgACCGGGGTGTGGGGCGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACCGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           <  1:539645/150‑1 (MQ=37)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:2238483/1‑150 (MQ=255)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:2238312/1‑150 (MQ=255)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  1:2611834/1‑150 (MQ=255)
                                                                                                                                             ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:2032791/1‑150 (MQ=255)
                                                                                                                                             ctccccctTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          <  1:3293555/150‑1 (MQ=255)
                                                                                                                                              tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGGTGTCATCAGTCTGGAAGCCGCACGTTgg         >  2:3407577/1‑150 (MQ=255)
                                                                                                                                              tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         >  1:143113/1‑150 (MQ=255)
                                                                                                                                               ccccctCCCGGGGCGGGGGGGGGGTTGAGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc        <  2:3704475/148‑1 (MQ=37)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGACGCCGCACGTTGGCt       >  1:3732789/1‑150 (MQ=255)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:3004824/1‑150 (MQ=255)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:2551379/1‑150 (MQ=255)
                                                                                                                                                tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:2543457/1‑150 (MQ=255)
                                                                                                                                                    ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:149998/1‑150 (MQ=255)
                                                                                                                                                    ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  1:2727302/1‑150 (MQ=255)
                                                                                                                                                     tCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAtt  <  1:1279652/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239462‑3239758

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