Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,996,377 T→A 7.6% L152F (TTA→TTT sdiA ← quorum‑sensing transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,996,3770TA 7.6% 257.0 / 8.7 91L152F (TTA→TTTsdiAquorum‑sensing transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base A (2/5);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

GCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGC  >  NC_000913/1996231‑1996521
                                                                                                                                                  |                                                                                                                                                
gCTATCTCTGCTGATGTTTTCCCTTCCGCCGGCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTaa                                                                                                                                                <  2:1518927/149‑1 (MQ=255)
   aTCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTg                                                                                                                                             <  1:1021474/149‑1 (MQ=255)
   aTCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTg                                                                                                                                             <  2:643686/149‑1 (MQ=255)
     ctctGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCa                                                                                                                                           <  1:1043625/149‑1 (MQ=255)
     ctctGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCa                                                                                                                                           <  2:996670/149‑1 (MQ=255)
      tctGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCGTGCTGAAATTCATCTCTGGCGTCATCACTAGTTCATCATTAAACCGCATCAGAGCCATCAGACTTTCGCGCACCAGTATCTCCATTTTTAATTGCAg                                                                                                                                          >  2:195838/1‑149 (MQ=255)
       ctgctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGt                                                                                                                                         <  1:123490/149‑1 (MQ=255)
       ctgctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGt                                                                                                                                         <  2:340907/149‑1 (MQ=255)
         gctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGtta                                                                                                                                       <  1:1493150/149‑2 (MQ=255)
          ctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTtca                                                                                                                                      >  1:1240979/1‑149 (MQ=255)
          ctgATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGttac                                                                                                                                      <  2:227223/149‑3 (MQ=255)
          ctgATGTTTTCCCTTACGCCTTACACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGTCTTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGttac                                                                                                                                      <  2:667421/149‑3 (MQ=255)
             aTGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTtcatca                                                                                                                                   >  2:373476/1‑149 (MQ=255)
             aTGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCAATTAAAAATGCAGttactgg                                                                                                                                   <  1:396484/149‑6 (MQ=255)
                ttttCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTa                                                                                                                                >  1:576034/1‑149 (MQ=255)
                    cccTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGaa                                                                                                                            <  1:545946/149‑1 (MQ=255)
                      cTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATg                                                                                                                          >  1:1362105/1‑149 (MQ=255)
                      cTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATg                                                                                                                          <  2:621340/149‑1 (MQ=255)
                       tttctgcgTGCAACTCAGAATTTCTTTTTCGCTCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCATAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTTCAGTTCATCACTAAGAATgg                                                                                                                         <  2:351746/143‑1 (MQ=255)
                       ttCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATgg                                                                                                                         >  1:1044015/1‑149 (MQ=255)
                       ttCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATgg                                                                                                                         <  2:249379/149‑1 (MQ=255)
                         ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCCTCAGACTTTCGCGAACCAGTAACTGCATTTATAATTGCAGTTCATCACTAAGAATGGGt                                                                                                                       >  1:1082549/1‑149 (MQ=255)
                         ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGt                                                                                                                       >  1:1416421/1‑149 (MQ=255)
                         ccgccgTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGt                                                                                                                       >  2:528722/1‑149 (MQ=255)
                              gTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTc                                                                                                                  <  2:1386835/149‑1 (MQ=255)
                                     tCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTcgcgcgcg                                                                                                           >  2:1584483/1‑149 (MQ=255)
                                     tCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTcgcgcgcg                                                                                                           >  2:869674/1‑149 (MQ=255)
                                       aGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCt                                                                                                         <  2:364197/149‑1 (MQ=255)
                                        gAATTTCTTTTTCGCGCTGGCTGAAATTCATCTCTGGCGGCATCACTATATCATCATTTAAACGCATCAGCGCCATAAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTg                                                                                                        <  2:320325/149‑1 (MQ=255)
                                         aaTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGc                                                                                                       >  1:1197014/1‑149 (MQ=255)
                                         aaTTTCTTTTTCGCGCTTGCTGAAAATCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCCTCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGc                                                                                                       <  1:119606/149‑1 (MQ=255)
                                          aTTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCa                                                                                                      >  1:1312839/1‑149 (MQ=255)
                                             tCTTTTTCGAGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCGGCAACg                                                                                                   <  1:818594/149‑1 (MQ=255)
                                              cTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACgg                                                                                                  <  1:1055301/149‑1 (MQ=255)
                                                       gcTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGAc                                                                                         >  1:418236/1‑149 (MQ=255)
                                                            cTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACaaaaa                                                                                    <  1:1551039/149‑1 (MQ=255)
                                                             tGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAg                                                                                   >  2:719421/1‑149 (MQ=255)
                                                               aaaTTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCACCGGGAAAAGGACAAAACGcc                                                                                 >  2:991846/1‑149 (MQ=255)
                                                                 aTTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGTGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCa                                                                               <  2:901383/149‑1 (MQ=255)
                                                                       ctctGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGccc                                                                         >  1:7722/1‑149 (MQ=255)
                                                                         ctGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCgg                                                                       <  1:551982/149‑1 (MQ=255)
                                                                          tGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGCCAAAAAGCCCAGCGCCCGGt                                                                      >  2:6436/1‑149 (MQ=255)
                                                                          tGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGt                                                                      >  2:540392/1‑149 (MQ=255)
                                                                          tGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTCCATTTTTAATTGCAGTTCATCACAAAGAATAGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCgga                                                                      >  2:1083857/1‑148 (MQ=255)
                                                                           ggCGTCATCACTATTTCATCATTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGttt                                                                    <  1:490940/149‑1 (MQ=255)
                                                                           ggCGTCATAACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTAATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGtt                                                                     <  2:1001169/149‑1 (MQ=255)
                                                                            gCGTCATCAATATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTAAATTGCAGATCATCACTAAGACTGCGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGttt                                                                    <  1:259201/149‑1 (MQ=255)
                                                                             cGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTg                                                                   <  2:1416421/149‑1 (MQ=255)
                                                                              gTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTgg                                                                  <  2:750180/149‑1 (MQ=255)
                                                                              gTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGATCATCACTAAGAATGCGGATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCgggtggg                                                                  <  2:1095465/149‑6 (MQ=255)
                                                                                catcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgca                                                                <  2:1312839/149‑1 (MQ=255)
                                                                                 atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgcag                                                               >  1:304987/1‑149 (MQ=255)
                                                                                 atcaCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGgcag                                                               <  2:891116/149‑1 (MQ=255)
                                                                                  tcaCTATTTCATAATTTAAACTCCTCAGCGCCATCAGACATTAGCGCGCCACTAACTGCATTTTTAATTGCAGTTCAACACTAAGAATGGGTATTTCTCCCTCGCACCAACGGGAAAAGCCCAAAAAGCCCAGCGACCGGATTCgcagc                                                              >  2:1289672/1‑149 (MQ=255)
                                                                                    aCTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGTAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCAt                                                            >  2:733630/1‑149 (MQ=255)
                                                                                     cTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCAtt                                                           >  2:1269489/1‑149 (MQ=255)
                                                                                        tttCATCATTTCTACGCATCAGAGCCATCAGACTTTCGCGAACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCGGCAACGGGAAAAGGACAAAACGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                        <  1:927448/149‑1 (MQ=255)
                                                                                        tttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                        <  1:686130/149‑1 (MQ=255)
                                                                                        tttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                        <  1:373476/149‑1 (MQ=255)
                                                                                        tttCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTaaa                                                        <  1:1058236/149‑1 (MQ=255)
                                                                                            atcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACt                                                    >  1:1434976/1‑149 (MQ=255)
                                                                                             tcatTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTg                                                   >  1:1439658/1‑149 (MQ=255)
                                                                                              catTTAAACGCATCAGAGCCATCAGACTTTCGCGCAGCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGa                                                  <  2:233470/149‑1 (MQ=255)
                                                                                               atTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAATGCCACGCGCCCGGTTTGGCAGCATTAAATACTGAg                                                 >  2:143013/1‑149 (MQ=255)
                                                                                               atTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACAGCAATTTTAATTGCAGTTCAACACTAAGAATGGGAATTTCGCGCGCGCTGCAACGGGAAAACCACAAAAAGCCAAGCGCCCGCTTTGGCAGCATTAAAAACTGAg                                                 >  1:486846/1‑149 (MQ=255)
                                                                                                tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                >  2:1035170/1‑149 (MQ=255)
                                                                                                tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                <  1:829425/149‑1 (MQ=255)
                                                                                                tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                <  1:698772/149‑1 (MQ=255)
                                                                                                tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                <  2:1362105/149‑1 (MQ=255)
                                                                                                tttAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCAATAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGt                                                >  2:305881/1‑149 (MQ=255)
                                                                                                   aaaaGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGAAATTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGac                                             <  1:1462651/149‑1 (MQ=255)
                                                                                                     aCGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGacac                                           >  2:826559/1‑149 (MQ=255)
                                                                                                         aTCAGAGCCATCAGACTTTCGAGCACCAGTAACTGCATTTTTAATTGCAGTTAATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACcgcg                                       >  2:1518722/1‑149 (MQ=255)
                                                                                                                ccATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTaaa                                >  1:1094040/1‑149 (MQ=255)
                                                                                                                 cctCAGACTTTCGCGCACCAGTAACTGCATTTATAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAAc                               <  1:775512/147‑1 (MQ=255)
                                                                                                                 cATGAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAAc                               <  1:419064/149‑1 (MQ=255)
                                                                                                                  aTCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAAcc                              >  2:664564/1‑149 (MQ=255)
                                                                                                                   tCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCa                             >  2:589382/1‑149 (MQ=255)
                                                                                                                      gACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgt                          <  2:1240979/149‑1 (MQ=255)
                                                                                                                       aCTTTCGCTCACCAGTAACTGGATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGGGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCAtgtg                         <  2:1082549/149‑1 (MQ=255)
                                                                                                                         tttCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTgcg                       >  1:1420913/1‑149 (MQ=255)
                                                                                                                                 aCCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGc               >  1:1440773/1‑149 (MQ=255)
                                                                                                                                    atgaaCTGCAATTAAAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAACGCCCAGCGCCCGGTTTGGCAGCAGTAAATACTGAGTGACACCGCGGCGTCACCCCTGTGCGCGCGCGGCTTc            >  2:137839/4‑149 (MQ=255)
                                                                                                                                    atgaaCTGCAATTAAAAATGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTc            >  2:1353680/4‑149 (MQ=255)
                                                                                                                                    aGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTc            <  1:869674/149‑1 (MQ=255)
                                                                                                                                      tAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTccc          <  2:1044015/149‑1 (MQ=255)
                                                                                                                                        aCTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCAt        >  2:490521/1‑149 (MQ=255)
                                                                                                                                         cTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATa       >  2:1260829/1‑149 (MQ=255)
                                                                                                                                          tGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGACAAGCGCCAGGTTTGGCAGCATTAAATACTGAGTGACAACGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATaa      >  1:58431/1‑149 (MQ=255)
                                                                                                                                           gCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAAc     >  1:690770/1‑149 (MQ=255)
                                                                                                                                            cATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCACCGGGAAAAGGCCAAAAAGCCCAGCGCCCGGTTTGGCAGCAGTAAAGACTGCGTGACACCGCGGCGTAGACCATGTGCGCGCGCGGCCTccaaacgcg    >  2:223664/1‑142 (MQ=255)
                                                                                                                                              tttttAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACTAAAAGCCCAGCGCCGGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGGACACCAGGTGCGCGCGCGGCTTCCAATAACggg  >  2:1055357/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTATCTCTGCTGATGTTTTCCCTTCCGCCGTCCACCTCAGAATTTCTTTTTCGCGCTTGCTGAAATTCATCTCTGGCGTCATCACTATTTCATCATTTAAACGCATCAGAGCCATCAGACTTTCGCGCACCAGTAACTGCATTTTTAATTGCAGTTCATCACTAAGAATGGGTATTTCGCGCGCGCTGCAACGGGAAAAGGACAAAAAGCCCAGCGCCCGGTTTGGCAGCATTAAATACTGAGTGACACCGCGGCGTAAACCATGTGCGCGCGCGGCTTCCCATAACGGC  >  NC_000913/1996231‑1996521

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