Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,610 C→G 5.1% intergenic (+65/‑334) 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,6100CG 5.1% 528.3 / 15.2 153intergenic (+65/‑334)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (86/59);  new base G (2/6);  total (88/65)
Fisher's exact test for biased strand distribution p-value = 7.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 155 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:3773513/1‑150 (MQ=255)
         tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:2176401/1‑150 (MQ=255)
         tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  2:3104831/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)
               aTTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                                                                                                                      >  1:757017/1‑150 (MQ=255)
                ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:2977159/1‑150 (MQ=255)
                ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:1143240/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                                                                                                                            >  1:2991373/1‑150 (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)
                                              gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTGGGGGTGAGGGTAATTTTTCGCTCCGATGCTGGcg                                                                                                       >  2:1176071/1‑149 (MQ=255)
                                                 tGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                                                                                    >  2:756128/1‑150 (MQ=255)
                                                          ggCATTATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTATGCCGTGTACCCTGTCATTAACATGAGCAGCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           <  2:2645139/150‑1 (MQ=255)
                                                           gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:1148521/1‑150 (MQ=255)
                                                           gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  1:1920543/1‑150 (MQ=255)
                                                                    gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCGAAcccac                                                                                 >  1:2334125/1‑148 (MQ=255)
                                                                      tggggggTGTATATTTTTTAATCTGCCTACGCCGTGTACCCTGTCATTAACATGTGAACCGTTTTCCCCCTCTCCCTTCCCGGTAGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2315174/147‑1 (MQ=25)
                                                                      cggggggTGTATAATTTTTAATCTGCCTAAGCGGTTTACCCTTTATTTAACTTGAGCACGGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTGGGTGTGAGGGTAAATTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2434783/147‑1 (MQ=21)
                                                                      cggggggTGGATAATTTTTAATCTGCCTATGCGGTGTACCGTGTCTTTAACTTGTGCACGGTTTTCTCCCTCTCCCTTCCCGGGTGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:3841792/147‑1 (MQ=25)
                                                                      cgcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCGCTCCCTCCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  2:2434830/149‑1 (MQ=255)
                                                                      agcGGGGTGTATAATTTTTAATCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCTTTTTCTCCCTCTCCCTCCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                               <  1:2789954/150‑1 (MQ=38)
                                                                       ggggggggTATAATTTTTAATGTGCCTACGCGGTGTACCCTTTATTTAAAATGAGCACCGTTTTCACCCCCTCCCTGCCTGGGTGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                              <  1:2827780/148‑1 (MQ=25)
                                                                             tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        >  1:2859743/1‑150 (MQ=255)
                                                                             tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc                                                                        >  1:3280626/1‑150 (MQ=255)
                                                                               gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  2:3889112/1‑150 (MQ=255)
                                                                                aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                     >  1:3704475/1‑150 (MQ=255)
                                                                                  tttttttATTCTGCCTACGCGTTGTACCGTGAAATTAACATGCGCACGGTTTTCTCCCTCACCCTTCCCGGTAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   <  2:3147424/148‑1 (MQ=21)
                                                                                  aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   >  1:2625126/1‑150 (MQ=255)
                                                                                      tttAATGTGCCTACGCCGTGTCCCCGGTCTTTAACAGGCGCCCGGTTTTCGCCCTCCCCTTCCCGGGTTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  2:1274261/150‑1 (MQ=21)
                                                                                      tttAATCGGCCTACGCCGTGTCCCCGGTCATTAACAGGCGCACCGTTTTCTCCCTCGCCCTTCCTGGGTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                               <  1:2426596/150‑1 (MQ=25)
                                                                                       ttATTCTGCCTATGCGTTTTACCCTTTCTTTAACTTGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGTGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAggg                                                              <  2:1647342/150‑1 (MQ=21)
                                                                                       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)
                                                                                                       gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtg                                              >  2:1835484/1‑149 (MQ=255)
                                                                                                        ggtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttt                                             <  2:631979/149‑1 (MQ=255)
                                                                                                           aCCCTTTCATTAACTTGAGCACCTTTTCCTCCCTCGCCCTTCCCGGGTGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:2004407/150‑1 (MQ=25)
                                                                                                           aCCCTGTCATAAACAGCCCCCCCGTTTTCACCCACGCCCTCCCCGGGAGTGGGGCGGGGTTTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                          <  1:2381880/150‑1 (MQ=18)
                                                                                                            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)
                                                                                                              ctgaaaTTAACATCCGCCCGGTTTTCCCCCCCCCCCTCCCTGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  1:3426296/145‑1 (MQ=21)
                                                                                                                  aaTTAACATGAGCACCTTTTACCCCCCCCCCCCGCCCGGGTGTGGGGGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtgttg                                   <  1:3182307/149‑1 (MQ=21)
                                                                                                                       aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatc                              >  2:1988272/1‑150 (MQ=255)
                                                                                                                       aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatc                              >  1:3157876/1‑150 (MQ=255)
                                                                                                                         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)
                                                                                                                              cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGGGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATTTTGTCATCAGTCTgg                       >  1:1858887/1‑150 (MQ=255)
                                                                                                                               aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATAAGTCTGGa                      >  2:942355/1‑150 (MQ=255)
                                                                                                                                cggTTTTCTCCCCCTCCCTTCCCGGGTGTGGGGGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     <  1:440248/148‑1 (MQ=37)
                                                                                                                                    tttCTCCCTCTCCGTGCCCGGTTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCg                 <  2:2737757/150‑1 (MQ=37)
                                                                                                                                    tttCTCCCTCGCCCTTCCCGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCg                 <  2:2737764/150‑1 (MQ=37)
                                                                                                                                       ctccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCAc              >  2:2901824/1‑150 (MQ=255)
                                                                                                                                        tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTGGTTGTCATCAGTCTGGAAGCCGCACg             >  1:1239594/1‑150 (MQ=255)
                                                                                                                                          cctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:3290697/1‑150 (MQ=255)
                                                                                                                                          cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  2:1577903/1‑150 (MQ=255)
                                                                                                                                          cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:2541137/1‑150 (MQ=255)
                                                                                                                                          cctctcccgACCGGGGTGTGGGGCGGGGTGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACCGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           <  1:539645/150‑1 (MQ=37)
                                                                                                                                           ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  2:2238483/1‑150 (MQ=255)
                                                                                                                                           ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg          >  1:2611834/1‑150 (MQ=255)
                                                                                                                                           ctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTTTTGTTTTCAACATTCAGGAAGCCGCACGTTg          >  1:2333125/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)
                                                                                                                                             ctccctTCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc        <  1:2930297/150‑1 (MQ=255)
                                                                                                                                             ccccctCCCGGGGCGGGGGGGGGGTTGAGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc        <  2:3704475/148‑1 (MQ=37)
                                                                                                                                              tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:3004824/1‑150 (MQ=255)
                                                                                                                                              tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:2543457/1‑150 (MQ=255)
                                                                                                                                              tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  1:2551379/1‑150 (MQ=255)
                                                                                                                                                  ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt   >  2:3797078/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: