Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,264,999 T→A intergenic (+29/‑13) hemA → / → prfA glutamyl tRNA reductase/peptide chain release factor RF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,264,9990TA98.5% 251.4 / NA 66intergenic (+29/‑13)hemA/prfAglutamyl tRNA reductase/peptide chain release factor RF‑1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (30/35);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAG  >  NC_000913/1264853‑1265148
                                                                                                                                                  |                                                                                                                                                     
ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCAt                                                                                                                                                   >  1:782713/1‑151 (MQ=255)
  cATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCAtt                                                                                                                                                  <  2:161796/150‑1 (MQ=255)
   aTGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTa                                                                                                                                                <  1:1203597/151‑1 (MQ=255)
       aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGcc                                                                                                                                            <  1:18412/151‑1 (MQ=255)
       aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGc                                                                                                                                             <  1:1132621/150‑1 (MQ=255)
          cTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTAt                                                                                                                                         >  2:1016655/1‑151 (MQ=255)
            gACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGa                                                                                                                                       >  1:644765/1‑151 (MQ=255)
            gACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGa                                                                                                                                       >  2:1334776/1‑151 (MQ=255)
             aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGa                                                                                                                                       <  2:318663/150‑1 (MQ=255)
              cTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAg                                                                                                                                     <  2:641249/151‑1 (MQ=255)
              cTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAg                                                                                                                                     >  2:391320/1‑151 (MQ=255)
               tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACCGGTCGCCCGTGACGGGGAGAACGAACGCCTGATTATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACACGGAGCATTTACGCCTATGAAGc                                                                                                                                    <  1:734481/151‑1 (MQ=255)
                aaCCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGcc                                                                                                                                   <  2:257506/151‑1 (MQ=255)
                  ccGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCtt                                                                                                                                 >  2:103724/1‑151 (MQ=255)
                    gCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCtt                                                                                                                                 >  2:1480872/1‑149 (MQ=255)
                    gCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCtt                                                                                                                                 >  1:1377607/1‑149 (MQ=255)
                      ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTAt                                                                                                                             <  2:715260/151‑1 (MQ=255)
                      ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTAt                                                                                                                             <  1:1370506/151‑1 (MQ=255)
                      ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTAt                                                                                                                             <  2:678961/151‑1 (MQ=255)
                         agccATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCg                                                                                                                           <  1:164779/148‑1 (MQ=255)
                         aTCCATTCGCCAACGAAATCACTTCAACAGGCCGCCCTTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGt                                                                                                                          >  1:822036/1‑151 (MQ=255)
                         aTCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCg                                                                                                                           <  2:1069390/150‑1 (MQ=255)
                         aTCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGt                                                                                                                          >  2:1219666/1‑151 (MQ=255)
                                  ccAACGAAATCGCTTCAACAGGCCGGCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTAGCGTTGCCAAAc                                                                                                                  <  1:1104301/150‑1 (MQ=255)
                                        aaaTCACTTCAACAGGGCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACATCCTCGGGCTGGAGTAGAAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGc                                                                                                           <  1:1413025/151‑1 (MQ=255)
                                          aTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGccc                                                                                                         <  2:708043/151‑1 (MQ=255)
                                              cTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCa                                                                                                     >  1:876342/1‑151 (MQ=255)
                                               ttCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCAt                                                                                                    <  2:1507150/151‑1 (MQ=255)
                                                             cGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaaga                                                                                      <  1:589093/151‑1 (MQ=255)
                                                             cGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaaga                                                                                      <  1:1484177/151‑1 (MQ=255)
                                                                 aCGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGtt                                                                                  <  1:894387/151‑1 (MQ=255)
                                                                       aTAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCg                                                                            <  2:738199/151‑1 (MQ=255)
                                                                       aTAACGAACGCCTGAATATTCTGCACGACAGCCTCGGGCTGGAGTGGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCg                                                                            <  1:1381723/151‑1 (MQ=255)
                                                                                 ccTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGt                                                                   <  1:857441/150‑1 (MQ=255)
                                                                                 ccTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGt                                                                   <  2:519500/150‑1 (MQ=255)
                                                                                   tGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGa                                                                 >  2:1308989/1‑150 (MQ=255)
                                                                                   tGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGa                                                                 >  2:704978/1‑150 (MQ=255)
                                                                                      atatTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGAcgcg                                                             >  1:764405/1‑151 (MQ=255)
                                                                                        atTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCa                                                           >  1:109154/1‑151 (MQ=255)
                                                                                           cTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAAc                                                        <  2:774240/151‑1 (MQ=255)
                                                                                            tGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTGCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACt                                                       <  1:224342/151‑1 (MQ=255)
                                                                                                      ccTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGAc                                              <  2:671038/150‑1 (MQ=255)
                                                                                                      ccTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGAc                                              >  1:51541/1‑150 (MQ=255)
                                                                                                           ggCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGaa                                        >  1:1312404/1‑151 (MQ=255)
                                                                                                                gagaAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGttt                                   >  1:800096/1‑151 (MQ=255)
                                                                                                                gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGttt                                   >  2:1177701/1‑151 (MQ=255)
                                                                                                                gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGttt                                   >  2:1117263/1‑151 (MQ=255)
                                                                                                                gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGttt                                   >  2:92306/1‑151 (MQ=255)
                                                                                                                gAGTAGCAGCACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCATGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGttt                                   >  1:1249914/1‑151 (MQ=255)
                                                                                                                     gCAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTcgcg                              >  1:143274/1‑151 (MQ=255)
                                                                                                                      cAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTcgcg                              <  1:1219090/150‑1 (MQ=255)
                                                                                                                      cAGTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTcgcg                              >  1:945878/1‑150 (MQ=255)
                                                                                                                        gTACATCATTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCAt                           >  2:558666/1‑151 (MQ=255)
                                                                                                                           catcatTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTAt                        <  2:1204466/151‑1 (MQ=255)
                                                                                                                           catcatTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTAt                        >  2:871613/1‑151 (MQ=255)
                                                                                                                             tcatTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCa                      <  2:951336/151‑1 (MQ=255)
                                                                                                                               atTTTCTTTTTTTACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCAcg                    >  1:714538/1‑151 (MQ=255)
                                                                                                                                     tttttttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAata              <  1:290374/151‑1 (MQ=255)
                                                                                                                                       tttttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATg            >  1:843101/1‑151 (MQ=255)
                                                                                                                                       tttttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATg            >  1:356266/1‑151 (MQ=255)
                                                                                                                                         tttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgcg          >  2:860613/1‑151 (MQ=255)
                                                                                                                                          ttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgcgc         <  2:1382236/151‑1 (MQ=255)
                                                                                                                                          ttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgcgc         <  2:730429/151‑1 (MQ=255)
                                                                                                                                          ttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgcgc         <  1:581368/151‑1 (MQ=255)
                                                                                                                                          ttACAGGGAGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgcgc         <  1:1036418/151‑1 (MQ=255)
                                                                                                                                                 gTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGCTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGCTTTCGCGCATTATCACGCGAATATGCGCAGTTAAg  >  1:856810/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAG  >  NC_000913/1264853‑1265148

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