Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,453,413 G→A 5.4% intergenic (‑143/+95) 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,4130GA 5.4% 481.7 / 8.5 148intergenic (‑143/+95)rpsJ/gspB30S ribosomal subunit protein S10/part of gsp divergon involved in type II protein secretion
Reads supporting (aligned to +/- strand):  ref base G (63/76);  new base A (2/6);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 150 total aligned reads displayed.

TCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTTC  >  NC_000913/3453264‑3453562
                                                                                                                                                     |                                                                                                                                                     
 cTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGa                                                                                                                                                      <  1:961056/150‑1 (MQ=255)
   tcATGAGACCAGAGGTCCATTAATTTTATAAACGAATATGATTACTCCTCAGACCCATTACGATTGTTGGGAGAGTGTAACCGTGCTTACGTAGGTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGAtt                                                                                                                                                    <  1:145687/149‑1 (MQ=255)
   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)
        agaCCAGAGCTCCAATTATTTTATAAAAGAAAATGATACTCATCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTACCAGAGTGGCAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACct                                                                                                                                              >  1:862872/1‑150 (MQ=255)
         gaCCAGAGCTCCAATTATTTTATAAACGATAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACct                                                                                                                                              >  2:1505800/1‑150 (MQ=255)
          aCCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTTGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctc                                                                                                                                             <  2:200277/150‑2 (MQ=255)
           ccAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgc                                                                                                                                            <  2:1702526/150‑1 (MQ=255)
           ccAGAGCTCCAATTATTTTATAAACGAAAATGATACTCCTCAGACACTTTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACctgct                                                                                                                                           <  2:862872/150‑1 (MQ=255)
             agagCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAATCTGAACctcagt                                                                                                                                          <  2:306723/150‑5 (MQ=255)
              gtgCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTCGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAAATCGGCTAACTGAGGTTCAGATTGAACCTGCTGt                                                                                                                                         <  1:139590/148‑1 (MQ=255)
              cagCTCCAATTATTTTAAAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTTGGAGCATTGTTAGGTAGCCAAATTCGCCTAACTGAGGTTCAATCTGAACctcagtt                                                                                                                                         <  1:1685271/149‑6 (MQ=255)
               agCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTc                                                                                                                                        <  2:367078/150‑1 (MQ=255)
                 cTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCaa                                                                                                                                      >  1:1337286/1‑150 (MQ=255)
                   ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                    >  1:994516/1‑150 (MQ=255)
                   ccAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACt                                                                                                                                    >  2:1222202/1‑150 (MQ=255)
                    cAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTACTGAGGTTCATATTGAACCTGCAGTCAACTaa                                                                                                                                  >  1:525411/1‑149 (MQ=255)
                      attattTTATAAACGAAAATGATTACTCCTCAGACCCTTTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACg                                                                                                                                 <  1:643250/150‑1 (MQ=255)
                         attTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACa                                                                                                                              <  2:668123/150‑1 (MQ=255)
                            ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGTCAAATTCGGCTAACTGAGGTTCAATCTGAACctcag                                                                                                                                           <  2:1531299/134‑4 (MQ=255)
                            ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGTCAAATTCGGCTAACTGAGGTTCAATCTGAACctcag                                                                                                                                           >  1:1531299/1‑131 (MQ=255)
                            ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGc                                                                                                                           >  2:1592398/1‑150 (MQ=255)
                            ttATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGc                                                                                                                           >  2:1544402/1‑150 (MQ=255)
                             tataAACGAAACTGATTCCTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGcc                                                                                                                          <  1:320415/150‑1 (MQ=255)
                                aaaCGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCcgc                                                                                                                       >  1:918841/1‑150 (MQ=255)
                                 aaCGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCcgcg                                                                                                                      >  1:1498365/1‑150 (MQ=255)
                                  aCGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAGCTGAGGTTCAGATTGAACCTGCTGTCAACTACGCCAAGCCCgcgc                                                                                                                     >  2:1376526/1‑150 (MQ=255)
                                   cGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCa                                                                                                                    <  2:850450/150‑1 (MQ=255)
                                    gAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCAt                                                                                                                   <  1:397223/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:943035/150‑1 (MQ=255)
                                                   tcAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAacc                                                                                                    <  1:996331/150‑1 (MQ=255)
                                                     aGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAaccac                                                                                                  <  2:1268489/150‑1 (MQ=255)
                                                         ccATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCGAACTGAGGTTCAGATTGAACCTGCTGTCATCTACGACAAGCCGGCGCAGTAGACAGACTTAGGCACCGaa                                                                                              >  2:423857/1‑150 (MQ=255)
                                                          cATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCAGCGCATTATACATACTTAACCACCGAAc                                                                                             >  2:1301603/1‑150 (MQ=255)
                                                            ttACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGaa                                                                                              >  2:354140/1‑147 (MQ=255)
                                                            ttACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGaa                                                                                              <  1:354140/147‑1 (MQ=255)
                                                             tACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCa                                                                                          >  1:442433/1‑150 (MQ=255)
                                                              aCGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCaa                                                                                         >  1:368077/1‑150 (MQ=255)
                                                               cGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAg                                                                                        >  2:937780/1‑150 (MQ=255)
                                                                 aTTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCg                                                                                      <  1:25011/150‑1 (MQ=255)
                                                                        ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                               >  2:1629309/1‑150 (MQ=255)
                                                                        ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                               <  2:662345/150‑1 (MQ=255)
                                                                        ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                               <  1:1263487/150‑1 (MQ=255)
                                                                        ggAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCt                                                                               <  1:1222202/150‑1 (MQ=255)
                                                                          agagTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTag                                                                             <  1:916954/150‑1 (MQ=255)
                                                                           gagTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTagt                                                                            <  2:918841/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)
                                                                                     gTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACa                                                                  >  2:1344001/1‑150 (MQ=255)
                                                                                       tCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATAGATACTTACCCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGc                                                                >  1:1415087/1‑150 (MQ=255)
                                                                                            cGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTATTAAAGTAAAAAGCCTTCg                                                           >  2:229841/1‑150 (MQ=255)
                                                                                            cGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCg                                                           <  1:727532/150‑1 (MQ=255)
                                                                                                gCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGc                                                                >  1:827827/1‑141 (MQ=255)
                                                                                                  tCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAt                                                     >  1:1079812/1‑150 (MQ=255)
                                                                                                    cccGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTAtct                                                   <  2:766412/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)
                                                                                                        aTTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCa                                               <  2:1337286/150‑1 (MQ=255)
                                                                                                          tGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAg                                             <  2:817007/150‑1 (MQ=255)
                                                                                                              aGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAgttcg                                         <  1:1505800/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:1522756/1‑150 (MQ=255)
                                                                                                                   tGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAg                                    <  1:1301603/150‑1 (MQ=255)
                                                                                                                   tGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAg                                    <  2:1079812/150‑1 (MQ=255)
                                                                                                                   cctcagttAGCCGAATTTGGCTAGCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGGATTGCTAGTAAAGTAAACAGCCTTCGGGGTATCTCTCAAGTTCGTTCAg                                    >  2:248394/8‑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)
                                                                                                                            gCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCAttttt                           <  2:368077/150‑1 (MQ=255)
                                                                                                                             ccAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTg                          <  2:1498365/150‑1 (MQ=255)
                                                                                                                              caatatCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTAAGACTAGCCCGCGCATTATACTTATTTAACCACCGAACGCTAGCGTTTTGCTAGTATAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTgt                         <  1:1077703/145‑2 (MQ=255)
                                                                                                                                aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                       <  1:1592398/150‑1 (MQ=255)
                                                                                                                                aaTTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGaca                       <  1:414948/150‑1 (MQ=255)
                                                                                                                                          aCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTa             <  1:1707114/150‑1 (MQ=255)
                                                                                                                                          aCTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTa             >  1:381452/1‑150 (MQ=255)
                                                                                                                                            tGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTaaa           >  2:706174/1‑150 (MQ=255)
                                                                                                                                             gAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTAAGTGCATTTTTGAAACGCTCTCTTAAAc          >  2:1420205/1‑150 (MQ=255)
                                                                                                                                               ggatCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCg        <  1:349622/147‑1 (MQ=255)
                                                                                                                                               ggTTCAGATTGCACCTGCTGTCAAGTACGACACGCCAGCGCATTATAGATACTTAACCACCGAACGAAAGCGTAATGCTAGTTATGTAAACAGCCTTCAGGTTATCTCTCAAGTTCGTTCTGTGGCTTTTTGACACGCTGTCGTaaccgg        >  2:252613/1‑146 (MQ=255)
                                                                                                                                                gTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCgg       <  2:525411/150‑1 (MQ=255)
                                                                                                                                                 ttCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGc      >  1:1536838/1‑150 (MQ=255)
                                                                                                                                                     gATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATAATTAACCACCGAACGCAAGCTTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTATAAAGTTAGTACAGTGCATTTTTTACACGCTCTCTTAAACCTGCttaa  >  1:248000/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCTTACGTAGCTCCCCGATTGGGAGCATTGTTAGGTAGCCAAATTCGGCTAACTGAGGTTCAGATTGAACCTGCTGTCAACTACGACAAGCCCGCGCATTATACATACTTAACCACCGAACGCAAGCGTATTGCTAGTAAAGTAAACAGCCTTCGGGTTATCTCTCAAGTTCGTTCAGTGCATTTTTGACACGCTCTCTTAAACCGGCTTTC  >  NC_000913/3453264‑3453562

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