Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,922 G→T 5.7% R816S (CGT→AGT)  rne ← ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,9220GT 5.7% 298.7 / 5.4 88R816S (CGT→AGT) rneribonuclease E
Reads supporting (aligned to +/- strand):  ref base G (39/44);  new base T (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCA  >  NC_000913/1141775‑1142071
                                                                                                                                                   |                                                                                                                                                     
cttcttCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGAt                                                                                                                                                     <  2:1514126/150‑1 (MQ=255)
 ttcttcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATa                                                                                                                                                    <  1:363783/150‑1 (MQ=255)
     tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcga                                                                                                                                                >  2:471534/1‑150 (MQ=255)
     tcCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGcga                                                                                                                                                >  1:704339/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:1201867/150‑1 (MQ=255)
            cgcgCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCgacgacgac                                                                                                                                         >  2:993245/1‑150 (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                                                                                                                         >  2:789339/1‑150 (MQ=255)
                            cTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTAcgcg                                                                                                                         <  1:1520240/150‑1 (MQ=255)
                            cTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGAT‑ACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACgcgc                                                                                                                        <  2:672213/150‑1 (MQ=255)
                              gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                       >  1:1117302/1‑150 (MQ=255)
                              gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                       <  2:515195/150‑1 (MQ=255)
                              gTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCa                                                                                                                       <  2:46438/150‑1 (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:1312753/1‑150 (MQ=255)
                                                cAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAgcggc                                                                                                     >  1:1084496/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                                                                                               <  1:416756/145‑1 (MQ=255)
                                                           cGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAac                                                                                               >  2:416756/1‑145 (MQ=255)
                                                            ggATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGAcgcg                                                                                         <  2:1496911/150‑1 (MQ=255)
                                                               tCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCa                                                                                      >  1:956419/1‑150 (MQ=255)
                                                                 cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                    >  1:688856/1‑150 (MQ=255)
                                                                 cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                    <  1:771147/150‑1 (MQ=255)
                                                                 cAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATg                                                                                    >  2:892266/1‑150 (MQ=255)
                                                                  aGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGc                                                                                   >  1:920374/1‑150 (MQ=255)
                                                                     cTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCAc                                                                                >  2:2875/1‑150 (MQ=255)
                                                                     cTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCAc                                                                                >  2:261430/1‑150 (MQ=255)
                                                                        tGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGt                                                                             <  2:1403899/150‑1 (MQ=255)
                                                                        tGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGt                                                                             >  1:1255140/1‑150 (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:682507/1‑150 (MQ=255)
                                                                                         ggAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCa                                                            >  2:1215292/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                                     <  2:1657284/150‑1 (MQ=255)
                                                                                                                cGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTgg                                     >  1:1299443/1‑150 (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                                               <  2:694567/112‑1 (MQ=255)
                                                                                                                                        gtcgtcgtcgCTATCGTGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTcttctt                                               >  1:694567/5‑116 (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)
                                                                                                                                                   |                                                                                                                                                     
CTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACGGTCAACGGCATTGGCGACTGGGTTGGATAACGCTCGTCACGATAGCGACGACGACGCTGACCACTTACGCGCAGGTGACGAGGCGAGCGGCGAGAACGACGCGGCATGCCACCGTTGTCACGGTTATCAGCATTGTTCTCTTCTTGCTGTTCTGGTGCAGTTTGCGCTACGACTGGCAGCGGGACTTTCA  >  NC_000913/1141775‑1142071

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