Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,414 A→T 5.4% intergenic (‑144/+94) rpsJ ← / ← gspB 30S ribosomal subunit protein S10/part of gsp divergon involved in type II protein secretion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,453,4140AT 5.4% 440.1 / 10.3 148intergenic (‑144/+94)rpsJ/gspB30S ribosomal subunit protein S10/part of gsp divergon involved in type II protein secretion
Reads supporting (aligned to +/- strand):  ref base A (65/75);  new base T (2/6);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 149 total aligned reads displayed.

CTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTT  >  NC_000913/3453265‑3453561
                                                                                                                                                     |                                                                                                                                                   
cTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGa                                                                                                                                                     <  1:961056/150‑1 (MQ=255)
  gCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGAtt                                                                                                                                                   <  2:521833/150‑1 (MQ=255)
  gCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGAtt                                                                                                                                                   >  1:1469856/1‑150 (MQ=255)
    aTGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGa                                                                                                                                                 <  2:40219/150‑1 (MQ=255)
      gagaCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGGAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAAc                                                                                                                                               <  1:332720/150‑1 (MQ=255)
       agaCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAAcc                                                                                                                                              >  2:397223/1‑150 (MQ=255)
       agaCCAGAGCTCCAATTATTTTATAAAAGAAAATGATACTCATCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTACCAGAGTGGCAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACct                                                                                                                                             >  1:862872/1‑150 (MQ=255)
          ccAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgc                                                                                                                                           <  2:1702526/150‑1 (MQ=255)
          ccAGAGCTCCAATTATTTTATAAACGAAAATGATACTCCTCAGACACTTTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgct                                                                                                                                          <  2:862872/150‑1 (MQ=255)
           cAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgct                                                                                                                                          >  1:991968/1‑150 (MQ=255)
           cAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctcag                                                                                                                                          <  1:1150969/150‑4 (MQ=255)
            agagCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgctg                                                                                                                                         <  2:428949/150‑1 (MQ=255)
             gtgCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTCGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAAATCGGCTAACTGAGGTTCAGATTGAACCTGCTGt                                                                                                                                        <  1:139590/148‑1 (MQ=255)
              agCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTc                                                                                                                                       <  2:367078/150‑1 (MQ=255)
                  ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                   >  2:335980/1‑150 (MQ=255)
                  ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                   >  2:1222202/1‑150 (MQ=255)
                   cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTa                                                                                                                                  >  2:25011/1‑150 (MQ=255)
                   cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTa                                                                                                                                  >  1:1514563/1‑150 (MQ=255)
                     attattTTATAAACGAAAATGATTACTCCTCAGACCCTTTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACg                                                                                                                                <  1:643250/150‑1 (MQ=255)
                       tattTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGAc                                                                                                                              <  2:13356/150‑1 (MQ=255)
                       tattTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATGGGTGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGAc                                                                                                                              >  1:573177/1‑150 (MQ=255)
                        attTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACa                                                                                                                             <  2:1532560/150‑1 (MQ=255)
                        attTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACa                                                                                                                             <  2:668123/150‑1 (MQ=255)
                           ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGc                                                                                                                          >  2:1544402/1‑150 (MQ=255)
                           ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGc                                                                                                                          >  2:1592398/1‑150 (MQ=255)
                              taAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCcg                                                                                                                       <  1:1656829/150‑1 (MQ=255)
                               aaaCGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCcgc                                                                                                                      >  1:918841/1‑150 (MQ=255)
                                 aCGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAGCTGAGGTTCAGATTGAACCTGCTGTCAACTACGCCAAGCCCgcgc                                                                                                                    >  2:1376526/1‑150 (MQ=255)
                                  cGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCa                                                                                                                   <  2:850450/150‑1 (MQ=255)
                                       aTGATTACTCCTCAGACCAATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATtata                                                                                                              <  2:1514563/150‑1 (MQ=255)
                                            tACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTatacatac                                                                                                         <  2:731854/150‑1 (MQ=255)
                                                cctcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTaa                                                                                                     >  2:34215/1‑150 (MQ=255)
                                                cctcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTaa                                                                                                     >  1:161673/1‑150 (MQ=255)
                                                  tcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAacc                                                                                                   <  1:996331/150‑1 (MQ=255)
                                                    aGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAaccac                                                                                                 <  1:394068/150‑1 (MQ=255)
                                                       cccATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGa                                                                                              >  1:1334907/1‑150 (MQ=255)
                                                         cATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAAc                                                                                            >  1:817007/1‑150 (MQ=255)
                                                           ttACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGaa                                                                                             <  1:354140/147‑1 (MQ=255)
                                                             aCGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCaa                                                                                        >  2:1422587/1‑150 (MQ=255)
                                                              cGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAg                                                                                       >  2:937780/1‑150 (MQ=255)
                                                               gATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGc                                                                                      >  2:1354248/1‑150 (MQ=255)
                                                               gATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGc                                                                                      <  2:354843/150‑1 (MQ=255)
                                                                aTTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCg                                                                                     <  1:25011/150‑1 (MQ=255)
                                                                       ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              <  1:335980/150‑1 (MQ=255)
                                                                       ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              <  2:662345/150‑1 (MQ=255)
                                                                       ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              <  1:1263487/150‑1 (MQ=255)
                                                                          gagTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagt                                                                           <  2:918841/150‑1 (MQ=255)
                                                                            gtgtAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaa                                                                         >  1:1405910/1‑150 (MQ=255)
                                                                                aaCCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaaagta                                                                     <  1:887992/150‑1 (MQ=255)
                                                                                aaCCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaaagta                                                                     >  2:184851/1‑150 (MQ=255)
                                                                                aaCCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaaagta                                                                     <  1:465900/150‑1 (MQ=255)
                                                                                 aCCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaaagtaa                                                                    <  2:994516/150‑1 (MQ=255)
                                                                                   cGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAAc                                                                  <  2:1469856/150‑1 (MQ=255)
                                                                                      tCTTCCGTAGCACCCCGATTGGGAGCATTGTTAGGAAGCCAAATACGGCTAACTGAGGTGCAGATTCAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGc                                                               <  2:1577146/150‑1 (MQ=255)
                                                                                           cGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTATTAAAGTAAAAAGCCTTCg                                                          >  2:229841/1‑150 (MQ=255)
                                                                                           cGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCg                                                          <  1:727532/150‑1 (MQ=255)
                                                                                             tAGCTCCCCGATTTGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCACGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCggg                                                        <  1:244167/150‑1 (MQ=255)
                                                                                               gCTCCCCGCTTTGGAGCATTGTTAGGTAGCTAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGGCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGtt                                                      <  2:1210038/150‑1 (MQ=255)
                                                                                                 tCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAt                                                    >  1:1079812/1‑150 (MQ=255)
                                                                                                  ccccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGtatct                                                   >  1:440557/1‑146 (MQ=255)
                                                                                                   cccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTACCCACCGAACGCAAGCGTATTGCTAGTAAAGTAACAGCGTTCGGGTTAtctc                                                 >  2:168502/1‑150 (MQ=255)
                                                                                                   cccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtct                                                  <  2:766412/150‑1 (MQ=255)
                                                                                                   cccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtct                                                  <  2:226273/150‑1 (MQ=255)
                                                                                                   cccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCAGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGTTTAtct                                                  >  2:613601/1‑150 (MQ=255)
                                                                                                     cGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtctct                                                <  1:1422587/150‑1 (MQ=255)
                                                                                                       aTTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCa                                              <  2:1337286/150‑1 (MQ=255)
                                                                                                         tGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAg                                            <  2:817007/150‑1 (MQ=255)
                                                                                                             aGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcg                                        <  1:1505800/150‑1 (MQ=255)
                                                                                                              gCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgt                                       <  1:491559/150‑1 (MQ=255)
                                                                                                              gCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgt                                       <  1:629565/150‑1 (MQ=255)
                                                                                                               cATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgtt                                      <  1:1286740/150‑1 (MQ=255)
                                                                                                               cATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgtt                                      >  1:1306830/1‑150 (MQ=255)
                                                                                                                aTTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgttc                                     <  2:1296272/150‑1 (MQ=255)
                                                                                                                aTTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgttc                                     <  1:559834/150‑1 (MQ=255)
                                                                                                                 ttGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCa                                    >  2:172937/1‑150 (MQ=255)
                                                                                                                 ttGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCa                                    >  2:1522756/1‑150 (MQ=255)
                                                                                                                  tGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAg                                   <  2:1079812/150‑1 (MQ=255)
                                                                                                                   gTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGt                                  >  2:1487866/1‑150 (MQ=255)
                                                                                                                    ttAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTg                                 <  2:655052/150‑1 (MQ=255)
                                                                                                                     tAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGc                                <  2:991968/150‑1 (MQ=255)
                                                                                                                        gTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTACAtt                             <  2:442433/150‑1 (MQ=255)
                                                                                                                         tAGCCAAATTAGGCTAACTGAGCTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCATCGTACGCAAGCGTATTGCTATTAATGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCAttt                            <  1:60718/150‑1 (MQ=255)
                                                                                                                            ccAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTg                         <  2:1498365/150‑1 (MQ=255)
                                                                                                                             caatatCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTAAGACTAGCCCGCGCATTATACTTATTTAACCACCGAACGCTAGCGTTTTGCTAGTATAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTgt                        <  1:1077703/145‑2 (MQ=255)
                                                                                                                               aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                      <  1:998100/150‑1 (MQ=255)
                                                                                                                               aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                      <  1:1487866/150‑1 (MQ=255)
                                                                                                                               aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                      <  1:414948/150‑1 (MQ=255)
                                                                                                                                      cTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGctctc               <  1:1544402/150‑1 (MQ=255)
                                                                                                                                       tAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGctctct              >  1:26275/1‑150 (MQ=255)
                                                                                                                                       tAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCACGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGctctct              <  1:650396/150‑1 (MQ=255)
                                                                                                                                         aCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTa            >  1:381452/1‑150 (MQ=255)
                                                                                                                                         aCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTa            <  1:1707114/150‑1 (MQ=255)
                                                                                                                                           tGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTaaa          >  2:706174/1‑150 (MQ=255)
                                                                                                                                           gcaGGTTCAATCTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTaaa          >  2:1296227/3‑150 (MQ=255)
                                                                                                                                              ggTTCAGATTGCACCTGCTGTCAAGTACGACACGCCAGCGCATTATAGATACTTAACCACCGAACGAAAGCGTAATGCTAGTTATGTAAACAGCCTTCAGGTTATCTCTCAAGTTCGTTCTGTGGCTTTTTGACACGCTGTCGTaaccgg       >  2:252613/1‑146 (MQ=255)
                                                                                                                                              ggTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCg       >  1:1447575/1‑150 (MQ=255)
                                                                                                                                               gTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCgg      <  2:525411/150‑1 (MQ=255)
                                                                                                                                                ttCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGc     >  1:1536838/1‑150 (MQ=255)
                                                                                                                                                ttCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGc     <  1:937780/150‑1 (MQ=255)
                                                                                                                                                   aGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCttt  >  1:356407/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTT  >  NC_000913/3453265‑3453561

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