Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,925 A→C 5.7% Y815D (TAT→GAT)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,9250AC 5.7% 274.9 / 7.9 88Y815D (TAT→GAT) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base C (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACC  >  NC_000913/1141776‑1142073
                                                                                                                                                     |                                                                                                                                                    
ttcttcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATa                                                                                                                                                      <  1:363783/150‑1 (MQ=255)
    tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcga                                                                                                                                                  >  1:704339/1‑150 (MQ=255)
    tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcga                                                                                                                                                  >  2:471534/1‑150 (MQ=255)
    tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcga                                                                                                                                                  <  1:1622160/150‑1 (MQ=255)
       tGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcgacga                                                                                                                                               >  1:1621945/1‑150 (MQ=255)
          tCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcgacgacga                                                                                                                                            >  1:971569/1‑150 (MQ=255)
           cgcgCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCgacgacgac                                                                                                                                           >  2:993245/1‑150 (MQ=255)
           cgcgCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCgacgacgac                                                                                                                                           <  2:1201867/150‑1 (MQ=255)
             cgcTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGGTAACGCTCGTCGCTATCGtgacgagcgtta                                                                                                                                         <  1:906170/150‑7 (MQ=255)
             cgcTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGc                                                                                                                                         <  1:1647301/150‑1 (MQ=255)
               ctgctCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTg                                                                                                                                       <  2:845131/150‑1 (MQ=255)
                     ttCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACt                                                                                                                                 <  1:993245/150‑1 (MQ=255)
                      tCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACtt                                                                                                                                <  2:971569/150‑1 (MQ=255)
                          ccTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTAcgc                                                                                                                            <  2:348523/150‑1 (MQ=255)
                           cTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTAcgcg                                                                                                                           <  1:1520240/150‑1 (MQ=255)
                           cTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTAcgcg                                                                                                                           >  2:789339/1‑150 (MQ=255)
                           cTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGAT‑ACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACgcgc                                                                                                                          <  2:672213/150‑1 (MQ=255)
                             gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                         <  2:515195/150‑1 (MQ=255)
                             gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                         <  2:46438/150‑1 (MQ=255)
                             gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                         >  1:1117302/1‑150 (MQ=255)
                             gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTGGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                         <  2:894157/150‑1 (MQ=255)
                                cATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGt                                                                                                                      >  1:1657284/1‑150 (MQ=255)
                                 aTCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTg                                                                                                                     <  2:1621945/150‑1 (MQ=255)
                                     tGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGa                                                                                                                 >  2:544720/1‑150 (MQ=255)
                                       cGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAgg                                                                                                               >  2:382987/1‑150 (MQ=255)
                                        ggACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGc                                                                                                              <  1:523309/150‑1 (MQ=255)
                                            gTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgc                                                                                                          >  2:771147/1‑150 (MQ=255)
                                             tACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcg                                                                                                         >  1:404725/1‑150 (MQ=255)
                                               cAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGTATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcggc                                                                                                       >  1:802571/1‑150 (MQ=255)
                                               cAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcggc                                                                                                       <  1:542865/150‑1 (MQ=255)
                                               cAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcggc                                                                                                       >  1:1084496/1‑150 (MQ=255)
                                               cAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcggc                                                                                                       >  1:1312753/1‑150 (MQ=255)
                                                  ttGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCgag                                                                                                    <  2:630995/150‑1 (MQ=255)
                                                  ttGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCgag                                                                                                    <  2:1084496/150‑1 (MQ=255)
                                                    ggATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGaa                                                                                                  >  2:1327706/1‑150 (MQ=255)
                                                         gCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAacgacg                                                                                             <  1:576334/150‑1 (MQ=255)
                                                          cGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAac                                                                                                 >  2:416756/1‑145 (MQ=255)
                                                          cGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAac                                                                                                 <  1:416756/145‑1 (MQ=255)
                                                           ggATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGAcgcg                                                                                           <  2:1496911/150‑1 (MQ=255)
                                                              tCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCa                                                                                        >  1:956419/1‑150 (MQ=255)
                                                                cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                      >  1:688856/1‑150 (MQ=255)
                                                                cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                      >  2:892266/1‑150 (MQ=255)
                                                                cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                      <  1:771147/150‑1 (MQ=255)
                                                                 aGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGc                                                                                     >  1:920374/1‑150 (MQ=255)
                                                                    cTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCAc                                                                                  >  2:2875/1‑150 (MQ=255)
                                                                    cTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCAc                                                                                  >  2:261430/1‑150 (MQ=255)
                                                                       tGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGt                                                                               >  1:1255140/1‑150 (MQ=255)
                                                                       tGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGt                                                                               <  2:1403899/150‑1 (MQ=255)
                                                                         ccAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTg                                                                             >  1:1213447/1‑150 (MQ=255)
                                                                           agagGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTc                                                                           >  1:686582/1‑150 (MQ=255)
                                                                           agagGCCAGTTACGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTc                                                                           >  2:1227731/1‑150 (MQ=255)
                                                                            gagGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCa                                                                          >  2:328291/1‑150 (MQ=255)
                                                                                 cAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGtt                                                                     <  1:977876/150‑1 (MQ=255)
                                                                                  aGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTa                                                                    >  2:380090/1‑150 (MQ=255)
                                                                                       cGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGc                                                               >  1:935034/1‑150 (MQ=255)
                                                                                        ggAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCa                                                              >  2:1215292/1‑150 (MQ=255)
                                                                                        ggAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCa                                                              >  2:682507/1‑150 (MQ=255)
                                                                                         gagaCGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCAt                                                             >  2:1489946/1‑150 (MQ=255)
                                                                                           gaCGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTg                                                           <  2:956419/150‑1 (MQ=255)
                                                                                            aCGCTCACGCTACGGTCATCGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAAGGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGt                                                          <  2:806080/150‑1 (MQ=255)
                                                                                              gcgcACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTtc                                                        >  1:925762/1‑150 (MQ=255)
                                                                                                gcACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTtctc                                                      <  2:1180671/150‑1 (MQ=255)
                                                                                                    gCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGGGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCtcttct                                                  <  2:704339/150‑1 (MQ=255)
                                                                                                    gCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCtcttct                                                  <  2:716346/150‑1 (MQ=255)
                                                                                                          gTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGt                                            <  1:328291/150‑1 (MQ=255)
                                                                                                             aaCGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCt                                         <  2:1117302/150‑1 (MQ=255)
                                                                                                              aCGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTg                                        >  2:1069883/1‑150 (MQ=255)
                                                                                                               cGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTgg                                       >  1:1299443/1‑150 (MQ=255)
                                                                                                               cGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTgg                                       <  2:1657284/150‑1 (MQ=255)
                                                                                                                ggCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGt                                      <  1:502436/150‑1 (MQ=255)
                                                                                                                   aTTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCa                                   <  2:272162/150‑1 (MQ=255)
                                                                                                                         gACTGGGTTGGTTTATGCTCGTCACGATAGCGACGACGACGGTGACCAGTTACGCGCTGGTGACGAGTCTAGCGTCGAGAACGACGCTGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTgc                             <  1:709342/150‑1 (MQ=255)
                                                                                                                         gACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTgc                             <  1:161754/150‑1 (MQ=255)
                                                                                                                          aCTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTgcg                            <  1:892266/150‑1 (MQ=255)
                                                                                                                           cTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACTCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGTTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTgcgc                           <  2:404725/150‑1 (MQ=255)
                                                                                                                              ggTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTAc                        >  2:1505562/1‑150 (MQ=255)
                                                                                                                                    aTAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTgg                  <  1:382987/150‑1 (MQ=255)
                                                                                                                                       gtcgtcgtcgCTATCGTGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTcttctt                                                 >  1:694567/5‑116 (MQ=255)
                                                                                                                                       gtcgtcgtcgCTATCGTGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTcttctt                                                 <  2:694567/112‑1 (MQ=255)
                                                                                                                                        tcgtcgtcgCTATCGTGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGc              >  2:140033/4‑150 (MQ=255)
                                                                                                                                        cGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGc              <  1:2875/150‑1 (MQ=255)
                                                                                                                                          gtcgtcgCTATCGTGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCgg            >  1:594317/2‑150 (MQ=255)
                                                                                                                                            cGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGa          <  1:471534/150‑1 (MQ=255)
                                                                                                                                            cGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTATTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGa          >  2:1417721/1‑150 (MQ=255)
                                                                                                                                               cACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACtt       >  2:1281617/1‑150 (MQ=255)
                                                                                                                                                 cGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTc     <  2:1299443/150‑1 (MQ=255)
                                                                                                                                                 cGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTc     <  2:920374/150‑1 (MQ=255)
                                                                                                                                                  gATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGAGGGACTTTCa    >  2:507386/1‑150 (MQ=255)
                                                                                                                                                    tAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCAcc  <  1:789339/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCACC  >  NC_000913/1141776‑1142073

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: ‑