Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,223 A→G 6.1% intergenic (‑71/+27) mutM ← / ← rpmG formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,811,2230AG 6.1% 347.6 / 10.2 114intergenic (‑71/+27)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base G (5/2);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 117 total aligned reads displayed.

GCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACC  >  NC_000913/3811076‑3811372
                                                                                                                                                   |                                                                                                                                                     
gCATGAAGAATGGTTGCACCAACGAGCTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGAAGCAGACAATATGCGGGCAAGCAGATAACAAAAAACCCCGCCGGACCGAgg                                                                                                                                                     >  2:1164704/1‑150 (MQ=255)
gCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCTGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAgg                                                                                                                                                     <  1:678939/150‑1 (MQ=255)
gCATGAACACTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATCTGGGGGCAAGCAGATAACAAAAACCCCCGCCGGAGCGAgg                                                                                                                                                     <  2:653503/150‑1 (MQ=255)
   tGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCTAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                                                                                                                                                  >  2:832760/1‑150 (MQ=255)
     aagaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAAATTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATTGGGGCAAGCAGATAACAAAAAACACCTCCGGAGCGAGGttttt                                                                                                                                                >  2:392510/1‑150 (MQ=255)
      agaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                                                                                                                                               >  1:1321745/1‑150 (MQ=255)
       gaaTGGTTGCACCAACGAGATGCGGTGCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                              <  2:1207740/150‑1 (MQ=255)
          tGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTa                                                                                                                                           >  2:947346/1‑150 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGCACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                          >  1:1169759/1‑150 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCATGCAGATAAGAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                          <  2:706958/150‑1 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTaa                                                                                                                                          >  2:1246035/1‑149 (MQ=255)
           ggTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                          <  2:1097771/150‑1 (MQ=255)
           ggTTGCACCAACGAGATGCGGTGCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATTGCATCTCCAGGAATGAACAGATGCAGTCACTATGGCGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                          <  2:467417/150‑1 (MQ=255)
            gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTatc                                                                                                                                         <  1:1286445/150‑3 (MQ=39)
            gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACa                                                                                                                                         >  2:669368/1‑150 (MQ=255)
              tGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGCTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATc                                                                                                                                       >  1:1066981/1‑150 (MQ=255)
               gCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                                      <  2:562007/150‑1 (MQ=255)
               gCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCa                                                                                                                                      <  1:938074/150‑1 (MQ=255)
                cACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTATTTTGTTACATCaa                                                                                                                                     >  1:238385/1‑150 (MQ=255)
                  ccATCGAGATGCGGTTCTATGCCGCGGCGGCTGGATTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATCATCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                                   <  1:976332/150‑1 (MQ=255)
                  ccAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                                   <  1:386979/150‑1 (MQ=255)
                    aaCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATCGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCACAGCg                                                                                                                                 >  2:825989/1‑150 (MQ=255)
                     aCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCga                                                                                                                                >  2:1197707/1‑150 (MQ=255)
                      cGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTATTGTTACATCAAAGCgag                                                                                                                               >  2:1323236/1‑150 (MQ=255)
                          aTGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatt                                                                                                                           <  2:967588/150‑1 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                       >  2:362454/1‑150 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                       >  2:340431/1‑150 (MQ=255)
                              ggTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAAATCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTAttt                                                                                                                       <  1:110455/150‑1 (MQ=255)
                               gTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTg                                                                                                                      <  2:1014991/150‑1 (MQ=255)
                               gTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTg                                                                                                                      <  1:514634/150‑1 (MQ=255)
                                          cggcggCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTGCAGGAATGAACAGATGCAGTGAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttc                                                                                                           <  1:308742/150‑1 (MQ=255)
                                            gcggcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGcttctt                                                                                                         <  1:547689/150‑1 (MQ=255)
                                               gcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGt                                                                                                      <  1:273763/150‑1 (MQ=255)
                                                cTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGTGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAAATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTa                                                                                                     >  1:660148/1‑150 (MQ=255)
                                                 tGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAg                                                                                                    <  2:303589/150‑1 (MQ=255)
                                                  ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   >  2:745612/1‑150 (MQ=255)
                                                  ggTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGa                                                                                                   >  2:846117/1‑150 (MQ=255)
                                                    tttCAACTTCGGGTATTTCAGTCATAGCATCTCCATGAATGAACTGATGCAGTCTATATGTGGGCAAGCAGATAACAAAATTCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGATAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                                 <  1:37741/150‑1 (MQ=255)
                                                    tttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAATAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                                 <  1:458730/150‑1 (MQ=255)
                                                    tttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCATTCAATATTGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTATCGCTTCTTTGTAGATc                                                                                                 >  2:689898/1‑150 (MQ=255)
                                                    tttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACACCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATc                                                                                                 >  1:582050/1‑150 (MQ=255)
                                                       cAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACg                                                                                              <  2:662621/150‑1 (MQ=255)
                                                          cTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGATGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGc                                                                                           >  1:892895/1‑150 (MQ=255)
                                                           ttCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGTTGCAGTCAATATGGGGGCAAGCCGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCt                                                                                          <  2:403806/150‑1 (MQ=255)
                                                             cGGGTAATTCAGTCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACTAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                        <  1:871172/150‑1 (MQ=255)
                                                             cGGGTAATTCAGGCATAGCATCTCCAGGAATCAACAGATGCAGTCAATATGGGGGAAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                        <  1:1039878/150‑1 (MQ=255)
                                                                  aaTTCAGGCATAGCATCTCCAGGAATGATCAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCCCGTGCTGGCGaac                                                                                   <  1:494198/150‑1 (MQ=255)
                                                                            tAGCATCTCCAGGAATGAACATATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGATCTGGCGAACAACTTGATCg                                                                         >  2:570674/1‑150 (MQ=255)
                                                                             aGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGa                                                                        <  1:667625/150‑1 (MQ=255)
                                                                              gCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGaa                                                                       >  1:810446/1‑150 (MQ=255)
                                                                               cATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAt                                                                      <  1:271569/150‑1 (MQ=255)
                                                                                aTCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCTAACAACTTCATCGAAtt                                                                     >  2:1013172/1‑150 (MQ=255)
                                                                                  ctcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttt                                                                   <  1:24758/150‑1 (MQ=255)
                                                                                   tcCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAttttt                                                                  >  1:240834/1‑150 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGTGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                                <  1:1340738/150‑1 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                                <  1:969945/150‑1 (MQ=255)
                                                                                     cAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                                <  2:1049510/150‑1 (MQ=255)
                                                                                       ggAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAACACCCCGCCGGAGCGAGTTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAg                                                              <  1:1164704/150‑1 (MQ=255)
                                                                                          aTGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTc                                                           >  1:423636/1‑150 (MQ=255)
                                                                                           tGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTcc                                                          <  2:165808/150‑1 (MQ=255)
                                                                                                 gATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                                    >  2:1066573/1‑150 (MQ=255)
                                                                                                  aTGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGttttt                                                   >  1:1294884/1‑150 (MQ=255)
                                                                                                      aGTCAATATGGGGGCAAGCAGATAACAAAACACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCgg                                               <  2:721631/150‑1 (MQ=255)
                                                                                                      aGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCgg                                               <  2:150817/150‑1 (MQ=255)
                                                                                                        tCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCt                                             <  1:589972/150‑1 (MQ=255)
                                                                                                         cAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCtt                                            >  2:233598/1‑150 (MQ=255)
                                                                                                         cAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCtt                                            >  1:109174/1‑150 (MQ=255)
                                                                                                          aaTATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTa                                           >  2:951919/1‑150 (MQ=255)
                                                                                                            tatGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGt                                         <  1:550759/150‑1 (MQ=255)
                                                                                                               gggggCAAGCAGATAAGAAAAAACCGCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACg                                      <  2:1079310/150‑1 (MQ=255)
                                                                                                                 gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTGGCTTGTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGGTTTTGCGGCTGAGTACGtt                                    >  2:117716/1‑150 (MQ=255)
                                                                                                                 gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGtt                                    <  1:947346/150‑1 (MQ=255)
                                                                                                                    cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGt                                 <  2:982446/150‑1 (MQ=255)
                                                                                                                      aGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTGCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCGGGCTTAGTACGTTTGtta                               >  1:843601/1‑149 (MQ=255)
                                                                                                                       tgatgtAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttct                              >  1:113795/6‑150 (MQ=38)
                                                                                                                       gCAGATAACAAAAAACCCCGCCGGAGCGAGGTGTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTGGTAGATCACGTGCTGGCGAACTACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttct                              >  1:33462/1‑150 (MQ=255)
                                                                                                                        cAGTTAAGAAAATACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGATTTTTTTGATTTTCTCTTCTTTGTAGATCACGTGCTGGCGAACTACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                             <  1:145728/150‑1 (MQ=255)
                                                                                                                        cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                             <  1:1075044/150‑1 (MQ=255)
                                                                                                                         aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttcttc                            <  2:1321745/150‑1 (MQ=255)
                                                                                                                           gtAACAATAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          >  1:1246022/2‑150 (MQ=255)
                                                                                                                           gtAACAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          >  2:265405/2‑150 (MQ=39)
                                                                                                                           aTAACAAAAAACCCCGCCGGAGCGATGTTTTTTGTTACATCAAAGCGAGTATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTTGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          <  1:669368/150‑1 (MQ=255)
                                                                                                                           aTAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGATTTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGt                          >  1:498103/1‑150 (MQ=255)
                                                                                                                            tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTa                         >  2:473304/1‑150 (MQ=255)
                                                                                                                              aCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt                       >  1:934429/1‑150 (MQ=255)
                                                                                                                               cAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATATTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTACAGTACCAGATTTTCCGGCTTAGTACGTTTGTTCTTCGGAGTg                      >  1:536014/1‑150 (MQ=255)
                                                                                                                               cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGtt                      >  2:879007/1‑149 (MQ=255)
                                                                                                                                aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTAGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg                     <  2:157763/150‑1 (MQ=255)
                                                                                                                                 aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt                    <  1:303013/150‑1 (MQ=255)
                                                                                                                                 aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt                    <  1:269432/150‑1 (MQ=255)
                                                                                                                                      ccccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGa               >  1:886238/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACACCTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGGTTGTTCTTCGTAGTGGTATAGaa              >  1:287201/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGaa              >  2:800943/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGaa              >  1:823693/1‑150 (MQ=255)
                                                                                                                                       cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTAGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGat              >  2:478406/1‑149 (MQ=255)
                                                                                                                                            cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAc         <  1:392510/150‑1 (MQ=255)
                                                                                                                                            cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGAc         <  2:1199919/150‑1 (MQ=255)
                                                                                                                                                   aGGTTTTTTGTTACATCAAATGGAGAATTATTTGTTTTTCGCTTCTTTGTAGATCACCTGCTGGCGAACAACTGGATCGAATTTTATCAGTTCCAGTTTTTCCGGCTTACTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc  <  1:229563/150‑1 (MQ=255)
                                                                                                                                                   aGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc  >  2:607877/1‑150 (MQ=255)
                                                                                                                                                   aGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc  <  2:912185/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACC  >  NC_000913/3811076‑3811372

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