Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,610 C→G 7.8% 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 7.8% 341.2 / 17.0 100intergenic (+65/‑334)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (57/35);  new base G (3/5);  total (60/40)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 102 total aligned reads displayed.

GGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239461‑3239757
                                                                                                                                                     |                                                                                                                                                    
gggCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTc                                                                                                                                                      <  2:112893/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:1774402/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:449104/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  2:2085784/150‑1 (MQ=255)
   cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   >  1:1986708/1‑150 (MQ=255)
   cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   >  1:472358/1‑150 (MQ=255)
    gtggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                                                                                  >  2:517453/1‑150 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              >  1:2210047/1‑150 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              >  1:126855/1‑150 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCCTCCAGGgaga                                                                                                                                              >  1:1013094/1‑150 (MQ=255)
         gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:435909/150‑1 (MQ=255)
         gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGagag                                                                                                                                             <  1:1134258/150‑1 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTTTCGGCATGATAAGCAGGGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTTCCCTGTCATTAACATGTGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  1:4238/1‑150 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            <  1:1250994/150‑1 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            <  1:1311590/150‑1 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  1:1405577/1‑150 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            >  1:2202057/1‑150 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                                                                                                            <  1:1889694/150‑1 (MQ=255)
            tGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCAAAGCCGTGTACCCTGACAATAACATGAGCAACGATTTCTCCCTCACCCTTCCAGGGAGAGGGt                                                                                                                                          >  1:117312/1‑150 (MQ=255)
            tGGCAGTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAggga                                                                                                                                          <  2:2394697/150‑2 (MQ=255)
                 ttCTGGTGATGCCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgggagaaaa                                                                                                                                     <  2:1552944/150‑7 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGTGTCgggg                                                                                                                                     >  2:1687884/1‑150 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:1687899/1‑150 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  2:2107636/1‑150 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgggagaaaa                                                                                                                                     <  2:2204801/150‑7 (MQ=255)
                    tGGTGATGACGTTTATTATCAACGCCTGGGTTAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:536229/150‑1 (MQ=255)
                           gaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAt                                                                                                                          >  2:1358798/1‑150 (MQ=255)
                              gTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGAGTGAGGGTAAtttt                                                                                                                       >  2:138475/1‑150 (MQ=255)
                                     atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                >  2:322754/1‑150 (MQ=255)
                                     atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                                                                                                >  2:1182813/1‑150 (MQ=255)
                                              tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgac                                                                                                       >  2:2307328/1‑148 (MQ=255)
                                              tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                                                                                       >  2:912778/1‑150 (MQ=255)
                                              tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATATTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCCCCGATGCTGGc                                                                                                       >  2:912792/1‑150 (MQ=255)
                                                    aaTTATCGGCTTTATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTTTACCCTGCCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                 <  2:2359449/150‑1 (MQ=255)
                                                    aaTTAGCGGCTTGATACGCTGGGGGTTGTATATTTTTTAATCTGCCTACGCCGTGTACCCTTTATTTAACTTGAGCACCGTTTTCTCCCTCCCCCTTCCTGGGTGCGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                 <  2:189502/150‑1 (MQ=255)
                                                       tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                              >  1:1230680/1‑150 (MQ=255)
                                                       aaTCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTATCTCCCTCTCCCTACCAGCGAGAGGGTCGCGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccca                                                                                              >  1:743780/2‑149 (MQ=255)
                                                            gcctGATCAGCCGCGGGGTGGATCCTTTTTAAGCTGCCTACGCCGTGTCCCCTGTCCTTCACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTGTCGCACCGCTGCTGGCGTGGTCCCCGCAcc                                                                                         >  2:2159088/4‑150 (MQ=255)
                                                                   aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc                                                                                  >  2:2103953/1‑150 (MQ=255)
                                                                   aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc                                                                                  >  2:77554/1‑150 (MQ=255)
                                                                       cggggggTGGTTATTTTTTAATCTGCCTACGCGGGGTACCCTGTCTTTAACTTGAGCACGGTTTTCTCCCTCTCCCTGCCTGGGTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc                                                                              <  2:2016658/147‑1 (MQ=25)
                                                                          gggtgggATAATTTTTACTCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGGAGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcccc                                                                           <  2:1363946/145‑1 (MQ=37)
                                                                             ttttatatTTTTTAATGTGCCAACGCGTTTTACCCTTTATTTAACTTGAGCACGTTTTTCCCCCACACCCTACCTGGTTGTGGTTGGGGTTTGGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                        <  2:524381/146‑1 (MQ=12)
                                                                                 ttattttttAATCTGCCTATGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    <  2:1987159/146‑1 (MQ=37)
                                                                                 ttAATTTTTAATCTGCCTATGCCGTGTACCCTGTCATTAACATGAGCACGGTTTTCTCCCTCTCCGTTCCGGGGAGAGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    <  2:1987145/149‑1 (MQ=37)
                                                                                 aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTTAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    <  2:1222401/150‑1 (MQ=255)
                                                                                 aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCTCCGATGCTGGCCTGTTCCCCTCACCCTATCCCTCTCCCCAAACggg                                                                    >  1:382166/1‑150 (MQ=255)
                                                                                 aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                    >  1:269841/1‑150 (MQ=255)
                                                                                      ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                               <  2:818222/150‑1 (MQ=255)
                                                                                       tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAgg                                                              >  1:878358/1‑150 (MQ=255)
                                                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                        >  1:1994279/1‑150 (MQ=255)
                                                                                             cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                        >  1:1994249/1‑150 (MQ=255)
                                                                                              ggcatacgCCGTGTACCCGGTCATTAACATGAGCACGGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                       <  1:2331456/143‑1 (MQ=25)
                                                                                              cgCCTAAGCCTTGTACCCTGTCATTAACATGAGCACCGTTTTCCCCCACTCCCTTCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                       <  1:2331472/149‑1 (MQ=25)
                                                                                                ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                     >  2:1813340/1‑150 (MQ=255)
                                                                                                ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                     >  2:131102/1‑150 (MQ=255)
                                                                                                  tAAGCCGTGTCCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGTGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                   <  2:1396635/150‑1 (MQ=255)
                                                                                                  tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                   >  2:981560/1‑150 (MQ=255)
                                                                                                          gtACCCTGTATTTAACATGCGCACCTTTTTCTCCCGCGCCCTACCTGGGTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttt                                           <  2:2116783/150‑1 (MQ=21)
                                                                                                              cctttaaTTAACAGGGGCACGTTTTTCTCCCTCGCCCTCCCTGGTTTTGGGTGGGGTGGTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                       <  2:689702/144‑1 (MQ=18)
                                                                                                              ccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                       >  2:1294735/1‑150 (MQ=255)
                                                                                                                tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCGGGCCTGTTCCCCGCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTGTTTTGAtg                                     >  2:942237/1‑150 (MQ=255)
                                                                                                                        aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatc                             >  2:2133777/1‑150 (MQ=255)
                                                                                                                          aTGGGCACCGGTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGCGTCCTGGTTTGATGGTGTCCTCAg                           >  2:39723/1‑150 (MQ=255)
                                                                                                                          aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAGCGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGGTGTCATCAg                           >  2:2043774/1‑150 (MQ=255)
                                                                                                                          aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAg                           >  2:1389572/1‑150 (MQ=255)
                                                                                                                               cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTAATCAGTCTgg                      >  1:2317280/1‑150 (MQ=255)
                                                                                                                                  ggTTTCCTCCCTCGCCCTTCCGGGGTGCGGGGGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   <  2:1928637/149‑1 (MQ=25)
                                                                                                                                  cGTTTTCTCCCTCTCCCTTCCCGGGTGCGGGGGGGGGTGCGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   <  1:1270213/150‑1 (MQ=37)
                                                                                                                                  cGTTTTCTCCCTCTCCCTTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   <  1:962323/150‑1 (MQ=255)
                                                                                                                                  cGTTTTCGCCCTCCCCCTTCCGGGGTGGGGGGGGGGGTGTGGTCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAg                   <  1:1635415/150‑1 (MQ=25)
                                                                                                                                    ccccgaccctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAAACCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGcc                 >  2:628408/7‑150 (MQ=255)
                                                                                                                                       cgaccctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCa              >  2:1519272/4‑150 (MQ=255)
                                                                                                                                       cgaccctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCa              >  2:1519277/4‑150 (MQ=255)
                                                                                                                                         tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTGTTTGATGTTGTCATCAGTCTGGAAGCCGCACg            >  1:1109713/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCCGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACCGGGGCGAGGGTACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  2:1425360/150‑1 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGCTGCTGGCCTGGTCCCCTCACCCAAACCCGCTCCCCCAACGGGGCGAGGGGCCGGCCCGCGTCCTTTTTTGCTGGGGTCATCAGGCTGGAAGCCGCACGGTgg        >  1:230477/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  2:1442869/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  1:347030/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  1:347031/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCCAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  2:1537139/150‑1 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        >  1:541517/1‑150 (MQ=255)
                                                                                                                                             tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCATTTTTTATTTTGTCATCAGTATTGAAGCCGCACGTTGGc       >  1:1938778/1‑150 (MQ=25)
                                                                                                                                             tcgccctccctGGGTGTGGGGGGGGGTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  1:1248755/147‑1 (MQ=37)
                                                                                                                                             cctccctGCCCGGGAGAGGGTCGGGGTGAGGTTAATTTTTCGCACCGATGCTGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  1:1149690/149‑1 (MQ=255)
                                                                                                                                             ccgccctCCCGGGGTGTGGGTGGGGGTGCGGTTCATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  1:1357330/147‑1 (MQ=37)
                                                                                                                                             cccccctGCCCGGGTGTGGGTGGGGGGGTGGGTCATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg        <  1:1357198/147‑1 (MQ=37)
                                                                                                                                              cgccctCCCGGGGTGTGGGTGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGc       <  1:73754/148‑1 (MQ=37)
                                                                                                                                               tccctTCCAGGGAGAGGGTCGGGGTGATGGTAATATTAAGCACCGATTCTGGCCTGTTACCCTCACCCTAAACCTCACCCCAAACGGGGCGAGGGGACTGACCGAGTACTTTTTTTATGTTGTCAAAAGTCTGGAAACCGCACGTAGGCt      >  2:25558/1‑150 (MQ=255)
                                                                                                                                               tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGGCCGAGTCCATTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt      >  1:1136961/1‑150 (MQ=255)
                                                                                                                                               tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCTCTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTGTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGGTGGCt      >  2:2036879/1‑150 (MQ=255)
                                                                                                                                               gccctTCCGGGGTGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACCGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt      <  2:2155586/149‑1 (MQ=255)
                                                                                                                                                ccctTCCTGGGTGTGGGTCGGGTTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCtt     <  1:1536650/150‑1 (MQ=38)
                                                                                                                                                ccctTCCCGGGTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCtt     <  2:1986708/150‑1 (MQ=255)
                                                                                                                                                ccctTCCCGGGAGTGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCtt     <  1:1003685/150‑1 (MQ=255)
                                                                                                                                                   ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGTTGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt  >  1:517060/1‑150 (MQ=255)
                                                                                                                                                   ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt  >  1:186973/1‑150 (MQ=255)
                                                                                                                                                   ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt  >  1:226961/1‑150 (MQ=255)
                                                                                                                                                   ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt  >  1:975857/1‑150 (MQ=255)
                                                                                                                                                   ttCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAt  >  1:186982/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAT  >  NC_000913/3239461‑3239757

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: