Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,610 C→G 5.2% 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.2% 1591.7 / 63.1 460intergenic (+65/‑334)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base C (267/168);  new base G (11/13);  total (279/181)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 463 total aligned reads displayed.

GGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239461‑3239758
                                                                                                                                                       |                                                                                                                                                    
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:7304201/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:1673269/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:4291189/150‑1 (MQ=255)
  gCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    >  1:3704965/1‑150 (MQ=255)
   cGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                                                                                   <  2:1797091/150‑1 (MQ=255)
     tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 <  2:6631948/150‑1 (MQ=255)
     tggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                 >  1:7673302/1‑150 (MQ=255)
      ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                <  2:9700476/150‑1 (MQ=255)
      ggtgTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                <  1:2279062/150‑1 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              >  1:124689/1‑150 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                              <  2:3314316/150‑1 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgaga                                                                                                                                              <  2:2366001/150‑1 (MQ=255)
        tgtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCGTCCAGGgaga                                                                                                                                              <  2:5745378/150‑1 (MQ=255)
          tttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCACCCTCTCCCTGGAAGGGAGAgg                                                                                                                                            <  1:9971012/150‑1 (MQ=255)
          tttGGCATGCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgg                                                                                                                                            <  2:7815489/150‑1 (MQ=255)
              gCATTCTGGTGTTGACGTTTATTATCAACGCCTGGGTGAATTATCGGAATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGGAGAgggaga                                                                                                                                        <  1:937350/150‑4 (MQ=255)
               cATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                       >  1:1403486/1‑150 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  1:3975666/1‑150 (MQ=255)
                 ttCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                                                                                                                     >  1:2719248/1‑150 (MQ=255)
                    tGGTGATGACGTTTATTATCAACGCGTGGGTGAATTATCGGCATGATAAGCAGGGGGGTGGATAATTTTTAATCTGCCTACGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  1:5877643/150‑1 (MQ=255)
                    tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:390525/150‑1 (MQ=255)
                    tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  <  2:73998/150‑1 (MQ=255)
                    tGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGa                                                                                                                                  >  1:6959842/1‑150 (MQ=255)
                       tgatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAggg                                                                                                                               >  1:964427/1‑150 (MQ=255)
                          tgaCGTTTATTATCACCGCTGGGTTTATTTTTGGGATTTATAAGCTGGGGGGTGGATAATTTTTAATCTGCCTACGCGGTGTACCCTTACTTTAACTTGAGCACCTTTTTCTCCCTCTCCCTTCCTGGGTGTGGGTCGGGGTGTGGGTaa                                                                                                                            <  1:8224673/150‑1 (MQ=255)
                          tgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGAGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTaa                                                                                                                            <  2:6508997/150‑1 (MQ=255)
                               tttATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                                                                                                                       >  1:6034834/1‑150 (MQ=255)
                                    tatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAc                                                                                                                  >  2:4494294/1‑150 (MQ=255)
                                            ccTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                                                                                                          >  2:2629770/1‑150 (MQ=255)
                                            ccTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                                                                                                          >  1:421081/1‑150 (MQ=255)
                                                  tGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                                                                                    >  1:8495670/1‑150 (MQ=255)
                                                    aaTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                  >  2:4343312/1‑150 (MQ=255)
                                                      ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                >  2:4750785/1‑150 (MQ=255)
                                                       tATCGGCATGATAAGCTGCGGGGTGTATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGTTGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                               <  2:1103662/150‑1 (MQ=255)
                                                           ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           >  1:4480136/1‑150 (MQ=255)
                                                           ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                                                                                           >  1:3441597/1‑150 (MQ=255)
                                                            gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGTTGAGGGTAATTATTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:8939212/1‑150 (MQ=255)
                                                            gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                          >  2:6756800/1‑150 (MQ=255)
                                                                tGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa                                                                                       >  1:6268369/1‑150 (MQ=255)
                                                                       ctggggggTGTATATTTTTTAATGTGCCTCCGCCGTGTACCCTGACATTAACATGCGCACCGTTTTCTCCCTCGCCCTTCCCGGGTGTGGGTGGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                                                <  2:8756544/146‑1 (MQ=25)
                                                                        tggggggTGTATAATTTTTATTCTCCCAATCCCGTGTACCGTGTATTTAACTTGAGCCCGGTTTTCTCCCTCCCCCTTCCCGGTTGTGGGTGGGGGTGCGGGTAAATTTTCGCACCGATGCTGGCCTGGTCCCCTCACCCTAACCCtctc                                                                               <  1:3165780/147‑1 (MQ=16)
                                                                         gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                              >  1:8785825/1‑150 (MQ=255)
                                                                              gtttaTAATTTTAAATCTGCCTACGCCGTGTCCCTTGTATTTAAATTGAGCACCTTTTTCCCCCTCCCCCTCCCCGGGTGTGGGTGGGGGTTTGGTAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                         <  2:1962880/146‑1 (MQ=16)
                                                                                 tatatTTTTTATTGTCCCTACGCCTTGTACCCTTTATTTAAATTGCGCCCGTTTTCCCCACCCTCCTTGCCGGGTTTTGGGTGGGGTTGGGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      <  2:3611384/149‑1 (MQ=12)
                                                                                 gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                    <  1:2363731/72‑1 (MQ=255)
                                                                                 gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg                                                                      >  2:8112326/1‑150 (MQ=255)
                                                                                  aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGCGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACggg                                                                     >  1:5689855/1‑150 (MQ=255)
                                                                                    attttttACTGGGCATAAGCGGTTTACCCTGACTTTAAATTGTGCACGGTTTGCTCCCTCGCCATACCTGGGTGTGGGTGGGGTTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   <  2:6068736/148‑1 (MQ=16)
                                                                                    attttttAATCTGCATAAGCCGTGTACCCTGTCTTTAACATGCGCACCTTTTTCTCCCTCGCCCTTCCTGGGTGTGGGTGGGGTTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGc                                                                   <  2:4878080/148‑1 (MQ=25)
                                                                                      tttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGa                                                                 >  1:944947/1‑150 (MQ=255)
                                                                                      tttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGa                                                                 <  2:8360942/150‑1 (MQ=255)
                                                                                       ttttAATCTCCCTACGCGGTTTACCCGGTAATTAACATGCGCACCGTTTACGCCCTCTCCCTTCCTGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAg                                                                <  1:7391153/150‑1 (MQ=25)
                                                                                           aaccTGCTTACGCGGTGTCCCGTTTAATTAACTTGCGCCCGTTTTTCCCCCCCCCCCGTCCGGGGTGTGGGGGGGGGTGCGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGa                                                            <  1:8879751/147‑1 (MQ=12)
                                                                                            aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGAc                                                           >  2:9845635/1‑150 (MQ=255)
                                                                                             tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACt                                                          >  2:6314785/1‑150 (MQ=255)
                                                                                               tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGa                                                        >  2:7283647/1‑150 (MQ=255)
                                                                                                 ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGAcc                                                      >  2:7097058/1‑150 (MQ=255)
                                                                                                 ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGGCGAGGGGCCTGAcc                                                      >  2:4268519/1‑150 (MQ=255)
                                                                                                 ccTAAGCCGTGTACCCTGAAATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCTGGGAGAGGGTCGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCACACGGGGCGAGGGGACTGAcc                                                      <  2:6571575/150‑1 (MQ=37)
                                                                                                  cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCg                                                     >  2:9507887/1‑150 (MQ=255)
                                                                                                   tAAGCGGTGTACCCTGTCATTAACTTGAGCACCTTTTTCTCCCTCTCCCTTCAGGGGAGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    <  2:1387611/150‑1 (MQ=25)
                                                                                                   tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                    >  2:5926908/1‑150 (MQ=255)
                                                                                                      gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTc                                                 <  1:8347915/150‑1 (MQ=255)
                                                                                                       ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTcc                                                >  2:2150780/1‑150 (MQ=255)
                                                                                                          tgaaCCCTGTCATTAACATGCCCCCGTTTTTCCCCCTCTCCCGGCCTGGTTGTGGGGGGGGTTGCGGTCAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCttt                                             <  1:4974761/147‑1 (MQ=21)
                                                                                                           ttACCCTGTAATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCCGGGTGTGGTTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttt                                            <  2:830538/149‑1 (MQ=37)
                                                                                                           gtccccTGTCATTAACATGTGCACCGTTGTCGCCCGCTCCCTGCCCGGGTGAGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttt                                            <  1:6514943/147‑1 (MQ=25)
                                                                                                              cccTGTCTTTAACCGGCGCCCCGTTTTCTCCCTCGCCCTGCGTGGTTGTGGGTGGGGGTGTGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTg                                         <  2:7921435/150‑1 (MQ=21)
                                                                                                               ccTGTCATAAACATGCGCCCGGTTTTCTCCCACTCCCTTCCGGGGTGTGGGTGGGGGTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                        <  1:5412779/150‑1 (MQ=25)
                                                                                                                cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  1:9855676/150‑1 (MQ=255)
                                                                                                                cTGTCATAAAAATGAGCCCCGTTTTCTCCCTCTCCCTACCCGTGAGCGGGTCGGGGTGCGTTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAt                                       <  1:4831800/150‑1 (MQ=25)
                                                                                                                 tGTCATTAACATGAGCACCGTTTTCTCACTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGCTGCTGGCCTGTTCCCCTCACCCTCACCCACTCCCCCAACGGGGCGAGGGGACTGACCGAGTCCTTGTTTTAtg                                      >  1:2751131/1‑150 (MQ=255)
                                                                                                                      ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATTTTGtc                                 >  2:7589832/1‑150 (MQ=255)
                                                                                                                        aaCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcat                               >  1:5921098/1‑150 (MQ=255)
                                                                                                                         aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcatc                              >  1:1026413/1‑150 (MQ=255)
                                                                                                                         aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCATTTTTGATGTTGtcatc                              >  1:8667641/1‑150 (MQ=255)
                                                                                                                             gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCGAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTc                          >  1:2540201/1‑150 (MQ=255)
                                                                                                                                  ccGTTTTCTCCCTATCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGaa                     >  2:5467408/1‑150 (MQ=255)
                                                                                                                                    gTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGATGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGc                   <  2:1655990/150‑1 (MQ=255)
                                                                                                                                    gTTTTCCCCCTCTCCCTTCCAGGTTGTGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGc                   <  2:5387069/150‑1 (MQ=25)
                                                                                                                                     ttttCTCACTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCACTACCCAAACGCGGCGAGGGGACTGGCAGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGcc                  >  2:190035/1‑150 (MQ=25)
                                                                                                                                         ctccctctccctTCCAGGGAGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGATGCCGCAc              <  2:1552341/150‑1 (MQ=255)
                                                                                                                                         ccccctccccgtTCCGGGGTGTGGGTGGGGGTGGGGTTAATTTTTCGCACCGATGCAGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCAc              <  2:1329128/148‑1 (MQ=25)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAATCGGGGCGAGGGGACTGACCGAGACCTTTTTTGAGGTGGTCATCCGGCTGGAAGCCGCACg             >  2:8541414/1‑150 (MQ=25)
                                                                                                                                          tccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAAGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTATGGTAGCCGCACg             >  1:7060979/1‑150 (MQ=255)
                                                                                                                                           ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGt            >  2:6265562/1‑150 (MQ=255)
                                                                                                                                            cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           >  1:7884012/1‑150 (MQ=255)
                                                                                                                                            cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCCTCAGTCTGGAAGCCGCACGtt           >  1:5133982/1‑150 (MQ=255)
                                                                                                                                            cctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTACCCTCACCCTAACCCTCTCCCCAAAAGGGGCGAGGGGACTGACCGAGTACTTTTTTTATGTTTTGATCAGTCTTGAAGCCGCACGtt           >  1:1982952/1‑150 (MQ=21)
                                                                                                                                            cctccccctTCATGGGATTGGTTCGGGGTTAGGCTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGtt           <  2:3200162/150‑1 (MQ=37)
                                                                                                                                              ctctaactTCCAGGGAGACGTTCGGGGTTAGGGTAATGCTTCGCACCGAAGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGCGGGGACTCACCCTGTCCTTGTTTGATGTTGTCCTCAGTCTTGCAGCCGCACGTTg          >  2:9281220/7‑150 (MQ=21)
                                                                                                                                               tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGCGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         >  1:3886262/1‑150 (MQ=255)
                                                                                                                                               tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAACGGGGCGCGGGGCCTGACCGAGTCCTTTTTTGATGGTGTCATCAGTCTGGAAGCCGCCCGTTgg         >  2:3679834/1‑150 (MQ=255)
                                                                                                                                               tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         >  2:3104853/1‑150 (MQ=255)
                                                                                                                                               cccccccTCCGGGGCGTGGGTGGGGGTGTGGGAAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         <  1:310866/147‑1 (MQ=37)
                                                                                                                                               actccctTCCGGGGTGTGGGTCGGGGTGTGGTTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         <  2:7056509/149‑1 (MQ=37)
                                                                                                                                               actccctTCCCGGGAGTGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTgg         <  2:9113575/149‑1 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  2:5279802/1‑150 (MQ=255)
                                                                                                                                                 tccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTGGTCATCAGTCTGGAAGCCGCACGTTGGCt       >  2:864475/1‑150 (MQ=255)
                                                                                                                                                      tCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAtt  >  1:4744116/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATT  >  NC_000913/3239461‑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: