Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,996,376 T→A 5.4% K153* (AAA→TAA)  sdiA ← quorum‑sensing transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,996,3760TA 5.4% 338.5 / 6.5 111K153* (AAA→TAA) sdiAquorum‑sensing transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (49/56);  new base A (4/2);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 116 total aligned reads displayed.

CATCGCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGCTGT  >  NC_000913/1996227‑1996524
                                                                                                                                                     |                                                                                                                                                    
catcGCTATCTCTGCTGATGTTTTCCCTTCAGCCGTCCACCTCAGAATTTCTTTTTCGCGATTGATGAAATTCATCTCTGGCGTCATCACGATTTCATCATATAAACGCATCAGAGACAGCAGACTTTCGCGCCCCAGTAACTGCAgttt                                                                                                                                                      >  2:1450027/1‑150 (MQ=255)
aatcGCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAtttt                                                                                                                                                      <  2:2002041/149‑1 (MQ=255)
 atcGCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCACTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAttttt                                                                                                                                                     <  1:1149433/150‑1 (MQ=255)
    gCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCAGATCAGAATTTCTTTTTCGAGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAAt                                                                                                                                                  >  1:311027/1‑150 (MQ=255)
      tATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTg                                                                                                                                                <  2:593850/150‑1 (MQ=255)
      tATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAAATTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGATCCATAAGACTTTAGCTCTCAATTAACTGCAATTTTAATTg                                                                                                                                                >  1:379934/1‑150 (MQ=255)
       aTCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGc                                                                                                                                               <  1:902413/150‑1 (MQ=255)
           ctgctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGtt                                                                                                                                           <  2:906203/150‑1 (MQ=255)
           ctgctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGtt                                                                                                                                           <  2:742163/150‑1 (MQ=255)
             gctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAACCGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTtca                                                                                                                                         >  1:553884/1‑150 (MQ=255)
              ctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTtcat                                                                                                                                        <  2:2125296/150‑1 (MQ=255)
              ctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGttact                                                                                                                                        <  1:1675833/150‑4 (MQ=255)
                gATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCTTGTTTAATTGCAGTtcatca                                                                                                                                      >  1:1069015/1‑150 (MQ=255)
                      ttCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGa                                                                                                                                >  2:1422606/1‑150 (MQ=255)
                         ccTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATg                                                                                                                             >  1:2102838/1‑150 (MQ=255)
                             ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTa                                                                                                                         <  1:1543257/150‑1 (MQ=255)
                             ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTa                                                                                                                         <  1:1262319/150‑1 (MQ=255)
                               gccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTAtt                                                                                                                       <  2:1074294/150‑1 (MQ=255)
                                 cgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTc                                                                                                                     >  2:184308/1‑150 (MQ=255)
                                  gTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTcg                                                                                                                    <  2:1040247/150‑1 (MQ=255)
                                        cTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTcgcgcgcg                                                                                                              >  1:1837066/1‑150 (MQ=255)
                                          cAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCATAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTAATCACTAAGAATGTGTATTTAGCGCGCGCt                                                                                                            >  1:788364/1‑150 (MQ=255)
                                           aGAATTTCTTTTTCGCGCTTGCTGAAATTAATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTg                                                                                                           <  2:784615/150‑1 (MQ=255)
                                            gAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGc                                                                                                          <  1:371988/150‑1 (MQ=255)
                                             aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCa                                                                                                         >  2:836609/1‑150 (MQ=255)
                                              aTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCaa                                                                                                        >  2:275054/1‑150 (MQ=255)
                                                 tCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACgg                                                                                                     >  1:2179677/1‑150 (MQ=255)
                                                      ttCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGaaaa                                                                                                <  2:2196520/150‑1 (MQ=255)
                                                       tCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAg                                                                                               <  1:828214/150‑1 (MQ=255)
                                                       tCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAg                                                                                               <  1:76347/150‑1 (MQ=255)
                                                          cgcTTGCGGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGAc                                                                                            <  1:1997217/150‑1 (MQ=255)
                                                              tGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACaaaa                                                                                        >  1:1544672/1‑150 (MQ=255)
                                                              tGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACaaaa                                                                                        >  1:2137753/1‑150 (MQ=255)
                                                                cTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACggg                                                                                                    <  2:149517/136‑1 (MQ=255)
                                                                cTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACggg                                                                                                    >  1:149517/1‑136 (MQ=255)
                                                                cTGAAATTCATCTCTGGCGTCAGCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACATTCCCGCACCAGTAACTGCATTTATAATTGCAGTTCATCACTAAGAATGGGTAATACGCGCGCGATGCAACGGGAAAAGGACAAAAAg                                                                                      >  2:51833/1‑150 (MQ=255)
                                                                 tGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGc                                                                                     <  1:1152883/150‑1 (MQ=255)
                                                                 tGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGc                                                                                     <  1:1160550/150‑1 (MQ=255)
                                                                   aaaaTCATCTCAGGCGTCATCACTATTTCAACATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATAGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCGGCAACGGGAAAAGGACAAACAGccc                                                                                   <  2:385740/150‑1 (MQ=255)
                                                                    aaTTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCa                                                                                  >  1:894779/1‑150 (MQ=255)
                                                                     aTTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAg                                                                                 <  1:1024013/150‑1 (MQ=255)
                                                                        cATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAgcgc                                                                              >  1:16937/1‑150 (MQ=255)
                                                                           ctctGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCg                                                                           >  2:2033528/1‑150 (MQ=255)
                                                                             ctGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGt                                                                         >  2:1318987/1‑150 (MQ=255)
                                                                                gCGTCATCACTATTTCATCATTTAAACGCATCAGAGCGATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAAGGGTAAAAGGACAAAAAGCCCAGCGCCCGGTTTg                                                                      <  1:2010240/150‑1 (MQ=255)
                                                                                  gTCATCACTTTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgc                                                                    <  1:778571/150‑1 (MQ=255)
                                                                                  gTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgc                                                                    <  1:206778/150‑1 (MQ=255)
                                                                                  gTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgc                                                                    <  1:952848/150‑1 (MQ=255)
                                                                                    catcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgcag                                                                  <  2:1030752/150‑1 (MQ=255)
                                                                                     atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgcagc                                                                 <  2:1214297/150‑1 (MQ=255)
                                                                                      tcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgcagca                                                                <  1:2033528/150‑1 (MQ=255)
                                                                                        aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCAtt                                                              >  1:1764495/1‑150 (MQ=255)
                                                                                         cTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTa                                                             <  2:2201931/150‑1 (MQ=255)
                                                                                           aTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                           >  1:663620/1‑150 (MQ=255)
                                                                                           aTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                           >  1:124372/1‑150 (MQ=255)
                                                                                           aTTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                           >  2:1708857/1‑150 (MQ=255)
                                                                                             ttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATa                                                         >  1:2179844/1‑150 (MQ=255)
                                                                                               catcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACt                                                       <  1:1285188/150‑1 (MQ=255)
                                                                                                atcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAATTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTg                                                      >  2:818690/1‑150 (MQ=255)
                                                                                                atcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGTGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTg                                                      <  1:564004/150‑1 (MQ=255)
                                                                                                atcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTg                                                      >  1:1884188/1‑150 (MQ=255)
                                                                                                   atTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                   <  2:605179/150‑1 (MQ=255)
                                                                                                   atTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTATTAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGAGCCCGATTTGTCAGCATTAAAAACTGAga                                                   >  2:1335066/1‑149 (MQ=255)
                                                                                                   atTTAAACGCATCAGAGCCATCAGACTATCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGAATTTCGCGCGCGCTGCATCGGGAAAAGGACAAAAAGCCCAGCGCCCGGGTTGGCAGCTTTAAATACTGAGt                                                   <  2:56109/150‑1 (MQ=255)
                                                                                                    tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTg                                                  <  2:898562/150‑1 (MQ=255)
                                                                                                    tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGAAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTg                                                  <  2:1971502/150‑1 (MQ=255)
                                                                                                      tAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGac                                                >  1:696413/1‑150 (MQ=255)
                                                                                                      tAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGac                                                >  1:466976/1‑150 (MQ=255)
                                                                                                             aTCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCgcgg                                         <  2:1684486/150‑1 (MQ=255)
                                                                                                              tCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCgcggc                                        >  2:1134483/1‑150 (MQ=255)
                                                                                                              tCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCgcggc                                        <  1:2023554/150‑1 (MQ=255)
                                                                                                                agagCCATCAGACTTTCGCGCACCAGTAACTGCATGTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGCCACAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCTGCGt                                      >  2:1302159/1‑150 (MQ=255)
                                                                                                                 gagCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTa                                     >  2:1642801/1‑150 (MQ=255)
                                                                                                                 gagCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTa                                     <  1:184308/150‑1 (MQ=255)
                                                                                                                  agCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTaa                                    >  1:1277043/1‑150 (MQ=255)
                                                                                                                   gCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAACGCCCAGCGCCCGGTTTGGCAGCATTATATACTGAGTGACACAGCGGCGTaaa                                   >  1:1220906/1‑150 (MQ=255)
                                                                                                                      aTCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCa                                <  2:1837066/150‑1 (MQ=255)
                                                                                                                      aTCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCa                                <  1:1630471/150‑1 (MQ=255)
                                                                                                                        cAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtg                              >  2:1868279/1‑150 (MQ=255)
                                                                                                                         aGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgt                             >  2:1159925/1‑150 (MQ=255)
                                                                                                                          gACTTTCGCGGACCAGTAACTGCATTTTTAATAGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCACCGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgtg                            <  2:1478939/150‑1 (MQ=255)
                                                                                                                          gACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgtg                            <  2:2011105/150‑1 (MQ=255)
                                                                                                                          gACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgtg                            <  1:275054/150‑1 (MQ=255)
                                                                                                                          gACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTGCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgtg                            <  1:109683/150‑1 (MQ=255)
                                                                                                                                cgcgCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTgcgcgcg                      >  1:494451/1‑150 (MQ=255)
                                                                                                                                cgcgCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATAATGAGTGACACCGCGGCGTAAACCATGTgcgcgcg                      >  2:558528/1‑150 (MQ=255)
                                                                                                                                    cACCAGTAACTGCATTTTTAATTGCAGTTCAACACTAAGAATGGGTATTTCGCGCGCGCGGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTTCGCGCGCgga                  >  1:1353509/1‑149 (MQ=255)
                                                                                                                                     aCCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCt                 <  2:1646332/150‑1 (MQ=255)
                                                                                                                                     aCCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCt                 >  2:1197155/1‑150 (MQ=255)
                                                                                                                                      tgatgaaCTGCAATTAAAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCtt                >  1:422386/6‑150 (MQ=255)
                                                                                                                                      ccAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCtt                >  1:1122838/1‑150 (MQ=255)
                                                                                                                                        atgaaCTGCAATTAAAAATGCAGTTAATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTTACACCGCGGCGTAAACCATGTGCGCTCGCGGCTTcc              >  1:2134443/4‑150 (MQ=255)
                                                                                                                                        aGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTcc              >  2:414200/1‑150 (MQ=255)
                                                                                                                                        aGTAACTGAATTTTTAATTGCAGTTAATCACTATCAATGGGTATTTCGCGAGCGCTGCAAGTGGAAAAGGACAAAAAGCCAAGAGCCCGGTATGGAAGCATTAAATACTGAGTGAAACCGAGCCGTATACAATGTGCGCGCGCGTCGTcc              >  2:463050/1‑150 (MQ=255)
                                                                                                                                          gaaCTGCAATTAAAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCa            >  1:15035/2‑150 (MQ=255)
                                                                                                                                                cATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACgg      >  1:1492267/1‑150 (MQ=255)
                                                                                                                                                cATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACgg      <  2:1752918/150‑1 (MQ=255)
                                                                                                                                                 aTTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGGAACGGGAAACGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGGAAACGAAGTGCGCGCGCGGCTTCCCATAACGGc     >  2:1767961/1‑150 (MQ=255)
                                                                                                                                                    tttAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGCtgt  >  1:149823/1‑150 (MQ=255)
                                                                                                                                                    tttAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGCtgt  <  1:2042529/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATCGCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGCTGT  >  NC_000913/1996227‑1996524

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: