Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,415 T→C 5.6% intergenic (‑145/+93) 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,4150TC 5.6% 425.5 / 10.8 143intergenic (‑145/+93)rpsJ/gspB30S ribosomal subunit protein S10/part of gsp divergon involved in type II protein secretion
Reads supporting (aligned to +/- strand):  ref base T (63/72);  new base C (2/6);  total (65/78)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 148 total aligned reads displayed.

TGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTT  >  NC_000913/3453266‑3453561
                                                                                                                                                     |                                                                                                                                                  
 gCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTACCTGAGGTTCAGAtt                                                                                                                                                   >  2:1707114/1‑150 (MQ=255)
 gCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGAtt                                                                                                                                                   <  2:521833/150‑1 (MQ=255)
 gCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGAtt                                                                                                                                                   >  1:1469856/1‑150 (MQ=255)
  cATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTg                                                                                                                                                  <  1:635874/150‑1 (MQ=255)
   aTGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGa                                                                                                                                                 <  2:40219/150‑1 (MQ=255)
    tGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTACCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGaa                                                                                                                                                >  2:1056866/1‑150 (MQ=255)
     gagaCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGGAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAAc                                                                                                                                               <  1:332720/150‑1 (MQ=255)
      agaCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAAcc                                                                                                                                              >  2:397223/1‑150 (MQ=255)
       gaCCAGAGCTCCAATTATTTTATAAACGATAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACct                                                                                                                                             >  2:1505800/1‑150 (MQ=255)
        aCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTTGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctc                                                                                                                                            <  2:200277/150‑2 (MQ=255)
         ccAGAGCTCCAATTATTTTATAAACGAAAATGATACTCCTCAGACACTTTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgct                                                                                                                                          <  2:862872/150‑1 (MQ=255)
          cAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctcag                                                                                                                                          <  1:1150969/150‑4 (MQ=255)
           agagCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgctg                                                                                                                                         <  2:428949/150‑1 (MQ=255)
           agagCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctcagt                                                                                                                                         <  2:306723/150‑5 (MQ=255)
            gtgCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTCGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAAATCGGCTAACTGAGGTTCAGATTGAACCTGCTGt                                                                                                                                        <  1:139590/148‑1 (MQ=255)
             agCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTc                                                                                                                                       <  2:367078/150‑1 (MQ=255)
                 ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                   >  2:1222202/1‑150 (MQ=255)
                 ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                   >  1:994516/1‑150 (MQ=255)
                 ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                   >  2:335980/1‑150 (MQ=255)
                  cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTa                                                                                                                                  >  1:1514563/1‑150 (MQ=255)
                  cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTa                                                                                                                                  >  2:25011/1‑150 (MQ=255)
                  cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTACTGAGGTTCATATTGAACCTGCAGTCAACTaa                                                                                                                                 >  1:525411/1‑149 (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)
                          ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGTCAAATTCGGCTAACTGAGGTTCAATCTGAACctcag                                                                                                                                          <  2:1531299/134‑4 (MQ=255)
                             taAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCcg                                                                                                                       <  1:1656829/150‑1 (MQ=255)
                                 cGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCa                                                                                                                   <  2:850450/150‑1 (MQ=255)
                                  gAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCAt                                                                                                                  <  1:397223/150‑1 (MQ=255)
                                      aTGATTACTCCTCAGACCAATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATtata                                                                                                              <  2:1514563/150‑1 (MQ=255)
                                         aTTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTatacat                                                                                                           >  1:605665/1‑150 (MQ=255)
                                           tACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTatacatac                                                                                                         <  2:731854/150‑1 (MQ=255)
                                            aCTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACt                                                                                                        <  2:155444/150‑1 (MQ=255)
                                               cctcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGCCAAGCCAGCGCATTATACATACTTaa                                                                                                     >  2:148670/1‑150 (MQ=255)
                                               cctcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTaa                                                                                                     >  2:34215/1‑150 (MQ=255)
                                               cctcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTaa                                                                                                     >  1:161673/1‑150 (MQ=255)
                                                 tcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAacc                                                                                                   <  1:996331/150‑1 (MQ=255)
                                                 tcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAacc                                                                                                   <  1:943035/150‑1 (MQ=255)
                                                   aGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAaccac                                                                                                 <  1:394068/150‑1 (MQ=255)
                                                   aGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAaccac                                                                                                 <  2:1268489/150‑1 (MQ=255)
                                                      cccATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGa                                                                                              >  1:1334907/1‑150 (MQ=255)
                                                        cATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAAc                                                                                            >  1:817007/1‑150 (MQ=255)
                                                        cATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCAGCGCATTATACATACTTAACCACCGAAc                                                                                            >  2:1301603/1‑150 (MQ=255)
                                                          ttACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGaa                                                                                             <  1:354140/147‑1 (MQ=255)
                                                           tACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCa                                                                                         >  1:442433/1‑150 (MQ=255)
                                                            aCGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCaa                                                                                        >  2:1422587/1‑150 (MQ=255)
                                                            aCGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCaa                                                                                        >  1:368077/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:1263487/150‑1 (MQ=255)
                                                                      ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              <  1:335980/150‑1 (MQ=255)
                                                                      ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              <  2:662345/150‑1 (MQ=255)
                                                                      ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                              >  2:1629309/1‑150 (MQ=255)
                                                                        agagTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTag                                                                            <  1:916954/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)
                                                                                aCCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagtaaagtaa                                                                    <  2:994516/150‑1 (MQ=255)
                                                                                   gTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACa                                                                 >  2:1344001/1‑150 (MQ=255)
                                                                                     tCTTCCGTAGCACCCCGATTGGGAGCATTGTTAGGAAGCCAAATACGGCTAACTGAGGTGCAGATTCAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGc                                                               <  2:1577146/150‑1 (MQ=255)
                                                                                            tAGCTCCCCGATTTGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCACGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCggg                                                        <  1:244167/150‑1 (MQ=255)
                                                                                              gCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGc                                                               >  1:827827/1‑141 (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)
                                                                                                   ccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtctc                                                 >  1:671805/1‑150 (MQ=255)
                                                                                                    cGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtctct                                                <  1:1422587/150‑1 (MQ=255)
                                                                                                        tGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAg                                            <  2:817007/150‑1 (MQ=255)
                                                                                                            aGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcg                                        <  1:1505800/150‑1 (MQ=255)
                                                                                                             gCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgt                                       <  1:629565/150‑1 (MQ=255)
                                                                                                             gCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgt                                       <  1:491559/150‑1 (MQ=255)
                                                                                                              cATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgtt                                      <  1:1286740/150‑1 (MQ=255)
                                                                                                              cATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcgtt                                      >  1:1306830/1‑150 (MQ=255)
                                                                                                                ttGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCa                                    >  2:172937/1‑150 (MQ=255)
                                                                                                                 tGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAg                                   <  2:1079812/150‑1 (MQ=255)
                                                                                                                 cctcagttAGCCGAATTTGGCTAGCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGGATTGCTAGTAAAGTAAACAGCCTTCGGGGTATCTCTCAAGTTCGTTCAg                                   >  2:248394/8‑150 (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)
                                                                                                                            caatatCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTAAGACTAGCCCGCGCATTATACTTATTTAACCACCGAACGCTAGCGTTTTGCTAGTATAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTgt                        <  1:1077703/145‑2 (MQ=255)
                                                                                                                              aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                      <  1:1487866/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:1707114/150‑1 (MQ=255)
                                                                                                                                          tGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTaaa          >  2:706174/1‑150 (MQ=255)
                                                                                                                                          gcaGGTTCAATCTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTaaa          >  2:1296227/3‑150 (MQ=255)
                                                                                                                                           gAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTAAGTGCATTTTTGAAACGCTCTCTTAAAc         >  2:1420205/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                  
TGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTT  >  NC_000913/3453266‑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: