Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,223 A→G 5.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 5.1% 359.4 / 6.6 120intergenic (‑71/+27)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base A (62/51);  new base G (4/2);  total (67/53)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 121 total aligned reads displayed.

CTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAC  >  NC_000913/3811074‑3811371
                                                                                                                                                      |                                                                                                                                                    
cTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGa                                                                                                                                                      <  1:1213592/150‑1 (MQ=255)
 tGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCATCGGCGGCTGGTTTCAACTTCGTGTAATTCAGGCATAGCATCTCCAGGATTGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAg                                                                                                                                                     <  2:685143/150‑1 (MQ=255)
  gCATGAAGAATGGTTGCACCAACGAGATGAGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAg                                                                                                                                                     >  1:1375170/1‑150 (MQ=255)
     tGAAGTATGGATGCACCAACTCGATGCGGTTCTATGCGGCGGCGGCTGGTTTCAACTTCGGGTAATTCTGGCATAGCATCTCAAGGAATGAAGAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCTCGCTGCGGCGGGGttt                                                                                                                                                 <  2:1253215/150‑1 (MQ=255)
     tGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                                                                                                                                                 >  1:1320219/1‑150 (MQ=255)
      gaagaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtttt                                                                                                                                                >  2:1198080/1‑150 (MQ=255)
         gaaTGGTTGCACCTACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                                                                                                                                             <  1:161942/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCTCGCTCCGGCGGGGTTTTTTg                                                                                                                                             <  1:1014891/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                                                                                                                                             <  2:1390358/150‑1 (MQ=255)
         gaaTGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                                                                                                                                             <  1:1414898/150‑1 (MQ=255)
         gaaTGGTGGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                                                                                                                                             <  1:585333/150‑1 (MQ=255)
             ggTTGCCCCAACGAGATGCGGGTCTATGCCGCGTCGGCTGGGTTCAACTTCGGTTAATTCATGCATAGCATCTCCAGGAATTAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                         <  2:1630617/150‑1 (MQ=255)
              gTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                                                                                                                                         <  2:1375170/150‑1 (MQ=255)
                   aCCAACGAGATGCCGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTAAGGCATAGCAACTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAACCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTCTTACATCaaa                                                                                                                                   >  2:649199/1‑150 (MQ=255)
                     cAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc                                                                                                                                 >  1:892511/1‑150 (MQ=255)
                     cAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc                                                                                                                                 >  1:865313/1‑150 (MQ=255)
                      aaCGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCg                                                                                                                                <  1:1015521/150‑1 (MQ=255)
                        cGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCgag                                                                                                                              <  2:527190/150‑1 (MQ=255)
                        cGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCgag                                                                                                                              <  2:520585/150‑1 (MQ=255)
                             tGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta                                                                                                                         <  1:818067/150‑1 (MQ=255)
                               cGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAattatt                                                                                                                       <  1:1504268/150‑1 (MQ=255)
                               cGGTTCTATGCCGCGGCCGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGTGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAattatt                                                                                                                       >  2:1142752/1‑150 (MQ=255)
                                     tATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGAttt                                                                                                                 >  2:77501/1‑150 (MQ=255)
                                      aTGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGAtttt                                                                                                                >  2:1364357/1‑150 (MQ=255)
                                        gCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTATTTTGTTACATCAAAGCGAGAATTATTTGATTTTCg                                                                                                              >  2:1640762/1‑150 (MQ=255)
                                          cgcgGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGct                                                                                                            <  2:817119/150‑1 (MQ=255)
                                                 gcTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGt                                                                                                     <  2:1557085/150‑1 (MQ=255)
                                                     gattCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGAt                                                                                                 >  2:862156/3‑150 (MQ=255)
                                                       ttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCa                                                                                               >  1:423173/1‑150 (MQ=255)
                                                       ttCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCa                                                                                               >  2:702328/1‑150 (MQ=255)
                                                         cAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGGCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACg                                                                                             <  2:88012/150‑1 (MQ=255)
                                                            cTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGc                                                                                          <  1:1198080/150‑1 (MQ=255)
                                                              tCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTg                                                                                        >  2:891558/1‑150 (MQ=255)
                                                               cGGGTAATTCAGGCATAGCATCTCCAGTAATGACCAGATGCAGTCAATAT‑GGGGGCAAGCAGATACCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATATGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                       <  2:325761/150‑1 (MQ=255)
                                                               cGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                       <  2:539332/150‑1 (MQ=255)
                                                               cGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                                       <  1:1339308/150‑1 (MQ=255)
                                                                 ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCTAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                                     <  2:798521/150‑1 (MQ=255)
                                                                 ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                                     >  2:1389916/1‑150 (MQ=255)
                                                                 ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                                     >  1:1548109/1‑150 (MQ=255)
                                                                  gTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGa                                                                                    >  2:818752/1‑150 (MQ=255)
                                                                   tAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGaa                                                                                   <  2:1323561/150‑1 (MQ=255)
                                                                    aaTTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGaac                                                                                  >  2:191972/1‑150 (MQ=255)
                                                                      ttCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGaacaa                                                                                >  2:1247864/1‑150 (MQ=255)
                                                                        cAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACt                                                                              >  2:245462/1‑150 (MQ=255)
                                                                         aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGTCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTACGCTTCTTTTTAGATCACGTCCTGCCGAACAACTg                                                                             >  2:336679/1‑150 (MQ=255)
                                                                         aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTg                                                                             >  2:594988/1‑150 (MQ=255)
                                                                         aGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTg                                                                             <  2:1320219/150‑1 (MQ=255)
                                                                         aGGAATAGCATATCCAGGAATGAACAGATTCAGACAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGTGTTTTTTTGTTACTTAAAAGCGAGAATTATTTGATTTTCGCTTGTTTGTAGATCACGAGCTGGCGAACACCTg                                                                             >  2:331638/1‑150 (MQ=255)
                                                                          ggCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTgg                                                                            >  1:63223/1‑150 (MQ=255)
                                                                           gCATAGCATCTAGAGGAATGAACAGATGCAGTCAATAT‑GGTGGCAATCAGATAACAAAAAACCCCGCCGGAGAGAGGTTTTTTGTTACATCAAAGCGATAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGa                                                                           <  1:84340/150‑1 (MQ=255)
                                                                            cATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGTGAGAATTATTTGATTTTCGCTTCTTTATAGATCACg                                                                                             >  1:660944/1‑131 (MQ=255)
                                                                            cATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGTGAGAATTATTTGATTTTCGCTTCTTTATAGATCACg                                                                                             <  2:660944/131‑1 (MQ=255)
                                                                             aTAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                                         <  2:128247/150‑1 (MQ=255)
                                                                              tAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCg                                                                        <  2:900208/150‑1 (MQ=255)
                                                                                 cATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAt                                                                     <  1:1308349/150‑1 (MQ=255)
                                                                                    ctcCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttt                                                                  >  1:1587923/1‑150 (MQ=255)
                                                                                       cAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTc                                                               <  2:1133074/150‑1 (MQ=255)
                                                                                          gAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGt                                                            >  1:453983/1‑150 (MQ=255)
                                                                                          gAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGt                                                            >  2:1708689/1‑150 (MQ=255)
                                                                                            aTGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTc                                                          <  2:578656/150‑1 (MQ=255)
                                                                                             tGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTcc                                                         <  2:423173/150‑1 (MQ=255)
                                                                                               aaCAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAg                                                       >  1:361228/1‑150 (MQ=255)
                                                                                               aaCAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAg                                                       >  1:1020269/1‑150 (MQ=255)
                                                                                                ccAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGCACAACTGGATCGAATTTTTTAAGTTCCAGt                                                      >  2:320313/2‑150 (MQ=255)
                                                                                                 cAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCTAATTTTTTCAGTTCCAGtt                                                     >  2:1323638/1‑150 (MQ=255)
                                                                                                 cAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtt                                                     >  2:675180/1‑150 (MQ=255)
                                                                                                   gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                                   >  2:172640/1‑150 (MQ=255)
                                                                                                   gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                                   >  2:449847/1‑150 (MQ=255)
                                                                                                   gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                                   >  1:1093232/1‑150 (MQ=255)
                                                                                                     tGCGGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTc                                                 >  2:1554447/1‑150 (MQ=255)
                                                                                                     tGCGGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTc                                                 >  2:1547668/1‑150 (MQ=255)
                                                                                                     tGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTc                                                 <  1:210029/150‑1 (MQ=255)
                                                                                                      gCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                                >  1:906654/1‑150 (MQ=255)
                                                                                                      gCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                                <  1:1364357/150‑1 (MQ=255)
                                                                                                      gCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                                <  1:191972/150‑1 (MQ=255)
                                                                                                         gTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGc                                             <  2:1027538/150‑1 (MQ=255)
                                                                                                         gTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGc                                             >  1:726892/1‑150 (MQ=255)
                                                                                                             atat‑GTGTGCAAGCAGATAACAAAAAACCCCGCCTGATCGAGTTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAg                                         <  2:1080131/150‑1 (MQ=255)
                                                                                                             atat‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAg                                         <  2:1176838/150‑1 (MQ=255)
                                                                                                                           cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttctt                            <  2:473078/150‑1 (MQ=255)
                                                                                                                            aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttcttc                           <  1:1327744/150‑1 (MQ=255)
                                                                                                                               tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTa                        >  1:665019/1‑150 (MQ=255)
                                                                                                                                aaCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAg                       >  1:756638/1‑150 (MQ=255)
                                                                                                                                aaCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAg                       >  2:378442/1‑150 (MQ=255)
                                                                                                                                 aCAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGCTCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGt                      >  1:1341787/1‑150 (MQ=255)
                                                                                                                                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                     <  1:818752/150‑1 (MQ=255)
                                                                                                                                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                     >  1:27142/1‑150 (MQ=255)
                                                                                                                                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                     >  1:1540200/1‑150 (MQ=255)
                                                                                                                                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGAGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACAGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg                     >  2:1207856/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTgg                    >  2:556341/1‑150 (MQ=255)
                                                                                                                                    aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt                   <  1:449847/150‑1 (MQ=255)
                                                                                                                                    aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTATTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGt                   >  2:1078641/1‑150 (MQ=255)
                                                                                                                                     aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGta                  <  1:336978/150‑1 (MQ=255)
                                                                                                                                      aaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGtat                 >  1:1713045/1‑150 (MQ=255)
                                                                                                                                         ccccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCGTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGa              <  1:1142752/150‑1 (MQ=255)
                                                                                                                                          cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCATCTTTGTAGATCACGTGCTGGCGAACAAATGGATCGAATTTTTTCAGTTCGAGTATTTCAGGCTTAGTCCGTTTGTTCTTCTTAGTGGTATAGaa             >  1:973359/1‑150 (MQ=255)
                                                                                                                                             gccgGTGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGTATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGTTATAGAAGTg          <  1:649199/150‑1 (MQ=255)
                                                                                                                                              ccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGa         >  2:1541289/1‑150 (MQ=255)
                                                                                                                                                    cGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTa   >  1:813485/1‑150 (MQ=255)
                                                                                                                                                     gAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAc  >  2:208388/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAC  >  NC_000913/3811074‑3811371

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