Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,192,685 A→C 5.2% R363R (CGT→CGG mnmA ← tRNA‑specific 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,192,6850AC 5.2% 233.9 / 4.2 77R363R (CGT→CGGmnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base C (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

ACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCGGT  >  NC_000913/1192539‑1192831
                                                                                                                                                  |                                                                                                                                                  
acaATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgg                                                                                                                                                 <  1:76556/150‑2 (MQ=255)
  aaTGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgctg                                                                                                                                               >  1:384710/1‑150 (MQ=255)
      ccccTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGTCTGCCgctgccgg                                                                                                                                           <  1:1364306/150‑4 (MQ=255)
      ccccTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCa                                                                                                                                           >  2:695396/1‑150 (MQ=255)
       cccTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCaa                                                                                                                                          >  2:820944/1‑150 (MQ=255)
        ccTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGTCTGCCgctgccggtc                                                                                                                                         <  1:1268645/150‑6 (MQ=255)
         cTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCaata                                                                                                                                        <  2:187049/150‑1 (MQ=255)
             tGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATAc                                                                                                                                    <  1:790002/150‑1 (MQ=255)
                    gttgttGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCAcc                                                                                                                             <  2:1203998/150‑1 (MQ=255)
                           cACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGcac                                                                                                                      >  1:993669/1‑150 (MQ=255)
                             ccAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCacac                                                                                                                    <  1:695396/150‑1 (MQ=255)
                                 gCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCa                                                                                                                <  2:286657/150‑1 (MQ=255)
                                   gTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCAcc                                                                                                              <  2:993669/150‑1 (MQ=255)
                                     gCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGt                                                                                                            >  1:956549/1‑150 (MQ=255)
                                       cGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGTTAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTta                                                                                                          >  1:886380/1‑150 (MQ=255)
                                        gACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTtat                                                                                                         >  1:4553/1‑150 (MQ=255)
                                          cTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATag                                                                                                       >  2:237774/1‑150 (MQ=255)
                                           tGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATaga                                                                                                      <  1:676172/150‑1 (MQ=255)
                                            gACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaa                                                                                                     >  2:1139432/1‑150 (MQ=255)
                                             aCAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaag                                                                                                    <  2:876527/150‑1 (MQ=255)
                                              cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                   >  1:1268051/1‑150 (MQ=255)
                                              cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                   >  2:1157711/1‑150 (MQ=255)
                                              cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                   >  1:1047742/1‑150 (MQ=255)
                                                aaTACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACg                                                                                                 >  1:117131/1‑150 (MQ=255)
                                                  tACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGc                                                                                               <  1:1283315/150‑1 (MQ=255)
                                                  tACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGACACACTTCACCGTTATAGAAGACGGc                                                                                               <  2:437470/150‑1 (MQ=255)
                                                   aCCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCa                                                                                              >  1:1159211/1‑150 (MQ=255)
                                                   aCCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCa                                                                                              <  1:835852/150‑1 (MQ=255)
                                                    ccGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCTGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGTCAg                                                                                             <  1:443609/150‑1 (MQ=255)
                                                       gCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCAGCGAGGCACACTTCACCGTTATAGAAGACGGCAGACt                                                                                          <  2:1099481/150‑1 (MQ=255)
                                                       gCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACt                                                                                          >  1:444493/1‑150 (MQ=255)
                                                          aGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGc                                                                                       >  2:1086083/1‑150 (MQ=255)
                                                           gggCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGcc                                                                                      <  2:1268051/150‑1 (MQ=255)
                                                           gggCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGcc                                                                                      <  1:946747/150‑1 (MQ=255)
                                                                aGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCg                                                                                 <  2:1346617/150‑1 (MQ=255)
                                                                      aTGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcgg                                                                           <  1:687842/150‑1 (MQ=255)
                                                                        gTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcggcg                                                                         >  2:1338100/1‑150 (MQ=255)
                                                                          cATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACggcggc                                                                       >  2:833479/1‑150 (MQ=255)
                                                                            tAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCaa                                                                     >  2:1275186/1‑150 (MQ=255)
                                                                            tAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCaa                                                                     >  1:909038/1‑150 (MQ=255)
                                                                              gTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAAcc                                                                   <  1:820944/150‑1 (MQ=255)
                                                                                aaTTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGTTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCgg                                                                 >  1:230452/1‑150 (MQ=255)
                                                                                 aTTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGt                                                                >  2:786041/1‑150 (MQ=255)
                                                                                 aTTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACAGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCCGCGTCACGGCGGCAACCGGt                                                                >  2:796603/1‑150 (MQ=255)
                                                                                   tCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTc                                                              <  1:1139432/150‑1 (MQ=255)
                                                                                    cTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCa                                                             <  2:984980/150‑1 (MQ=255)
                                                                                      ttGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATc                                                           <  1:1338100/150‑1 (MQ=255)
                                                                                        gCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGa                                                         <  2:1047742/150‑1 (MQ=255)
                                                                                             gTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATc                                                    <  1:96315/150‑1 (MQ=255)
                                                                                             gTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATc                                                    >  2:507237/1‑150 (MQ=255)
                                                                                               tCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCAc                                                  >  1:1356773/1‑150 (MQ=255)
                                                                                                  cTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTc                                               <  1:237774/150‑1 (MQ=255)
                                                                                                  cTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTc                                               <  1:562645/150‑1 (MQ=255)
                                                                                                    gCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCaa                                             >  1:669310/1‑150 (MQ=255)
                                                                                                      ttCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATg                                           <  1:507237/150‑1 (MQ=255)
                                                                                                            gTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGatca                                     >  2:1360554/1‑150 (MQ=255)
                                                                                                               ttAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCg                                  <  1:876609/150‑1 (MQ=255)
                                                                                                                tAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGt                                 >  2:122251/1‑150 (MQ=255)
                                                                                                                tAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGt                                 >  2:1081370/1‑150 (MQ=255)
                                                                                                                 aaGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTc                                >  1:1235319/1‑150 (MQ=255)
                                                                                                                    tAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAg                             >  2:173074/1‑150 (MQ=255)
                                                                                                                    tAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAg                             >  1:92077/1‑150 (MQ=255)
                                                                                                                     aaaGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAgc                            <  2:4553/150‑1 (MQ=255)
                                                                                                                      aaGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAgcg                           <  1:122251/150‑1 (MQ=255)
                                                                                                                          taataaTCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCtt                       >  2:1274096/1‑150 (MQ=255)
                                                                                                                          taataaTCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCtt                       >  1:1089710/1‑150 (MQ=255)
                                                                                                                             taatCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGAc                    >  1:642519/1‑150 (MQ=255)
                                                                                                                                 ttattgagcagcgtctgCCGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTg                >  2:369452/8‑150 (MQ=255)
                                                                                                                                     ccGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAg            <  1:1157711/150‑1 (MQ=255)
                                                                                                                                      gagcagcgtctgCCGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAgg           >  2:399951/3‑150 (MQ=255)
                                                                                                                                      cggcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAgg           <  2:444493/150‑1 (MQ=255)
                                                                                                                                       ggcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAggg          >  1:377681/1‑150 (MQ=255)
                                                                                                                                       ggcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAggg          <  2:1089710/150‑1 (MQ=255)
                                                                                                                                        gcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGa         <  2:117131/150‑1 (MQ=255)
                                                                                                                                          agcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGata       <  2:909038/150‑1 (MQ=255)
                                                                                                                                             ggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCg    >  1:1403436/1‑150 (MQ=255)
                                                                                                                                               cagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCGGt  >  2:809540/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCGGT  >  NC_000913/1192539‑1192831

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