Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,212 C→T 5.2% intergenic (‑60/+38) 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,2120CT 5.2% 367.2 / 4.9 116intergenic (‑60/+38)mutM/rpmGformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base C (59/50);  new base T (4/2);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

Only 100 of 122 total aligned reads displayed.

TGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAA  >  NC_000913/3811065‑3811360
                                                                                                                                                    |                                                                                                                                                    
tGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCccg                                                                                                                                                    >  1:1133074/1‑150 (MQ=255)
 gcgcACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCccgc                                                                                                                                                   <  1:437511/150‑1 (MQ=255)
  cgcACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTTGTTTCAACTTAGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCccgcc                                                                                                                                                  >  1:1250895/1‑150 (MQ=255)
      ccacTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAg                                                                                                                                              <  2:1059673/150‑1 (MQ=255)
       cacTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGACCAGATGCAGTCAATTT‑GGGGGCAAGCAGATAACAAAAAACCACGCGGGAGc                                                                                                                                             >  1:165349/1‑150 (MQ=255)
         cTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGa                                                                                                                                           <  1:1213592/150‑1 (MQ=255)
           gCATGAAGAATGGTTGCACCAACGAGATGAGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAg                                                                                                                                          >  1:1375170/1‑150 (MQ=255)
              tGAAGTATGGATGCACCAACTCGATGCGGTTCTATGCGGCGGCGGCTGGTTTCAACTTCGGGTAATTCTGGCATAGCATCTCAAGGAATGAAGAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCTCGCTGCGGCGGGGttt                                                                                                                                      <  2:1253215/150‑1 (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)
                             ccAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg                                                                                                                       >  2:210029/1‑150 (MQ=255)
                              cAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc                                                                                                                      >  1:865313/1‑150 (MQ=255)
                              cAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc                                                                                                                      >  1:892511/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)
                                              tATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATATCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCATATAACAAAAAACACCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATATGAttt                                                                                                      >  1:1459800/1‑150 (MQ=255)
                                               aTGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGAtttt                                                                                                     >  2:1364357/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)
                                                                   aaCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTAGATTTTCGCTTCTTTGTAGATCACGt                                                                                 >  1:297956/1‑150 (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                                                                            <  1:1339308/150‑1 (MQ=255)
                                                                        cGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTgg                                                                            <  2:539332/150‑1 (MQ=255)
                                                                          ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCTAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                          <  2:798521/150‑1 (MQ=255)
                                                                          ggTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCg                                                                          >  2:1389916/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                                                                   <  1:1701912/150‑1 (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)
                                                                                     cATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGTGAGAATTATTTGATTTTCGCTTCTTTATAGATCACg                                                                                  <  2:660944/131‑1 (MQ=255)
                                                                                     cATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGTGAGAATTATTTGATTTTCGCTTCTTTATAGATCACg                                                                                  >  1:660944/1‑131 (MQ=255)
                                                                                     cATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGAt                                                               >  1:1176838/1‑150 (MQ=255)
                                                                                      aTAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                              <  2:128247/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:1020269/1‑150 (MQ=255)
                                                                                                        aaCAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAg                                            >  2:700697/1‑150 (MQ=255)
                                                                                                        aaCAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATAGAATTTTTTCAGttaaag                                            >  2:822290/1‑146 (MQ=255)
                                                                                                         ccAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGCACAACTGGATCGAATTTTTTAAGTTCCAGt                                           >  2:320313/2‑150 (MQ=255)
                                                                                                          cAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCTAATTTTTTCAGTTCCAGtt                                          >  2:1323638/1‑150 (MQ=255)
                                                                                                          cAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtt                                          >  1:1338007/1‑150 (MQ=255)
                                                                                                          cAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtt                                          >  2:675180/1‑150 (MQ=255)
                                                                                                            gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                        >  1:1093232/1‑150 (MQ=255)
                                                                                                            gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                        >  2:449847/1‑150 (MQ=255)
                                                                                                            gATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGtttt                                        >  2:172640/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:191972/150‑1 (MQ=255)
                                                                                                               gCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTcc                                     >  1:906654/1‑150 (MQ=255)
                                                                                                                     aaTAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTa                               <  1:939347/150‑1 (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)
                                                                                                                                      gATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCg               <  1:700697/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)
                                                                                                                                           cAAAAAACCTCGCTCCGGCGGGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg          >  2:1634617/1‑150 (MQ=255)
                                                                                                                                           cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg          >  1:1540200/1‑150 (MQ=255)
                                                                                                                                           cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg          <  1:818752/150‑1 (MQ=255)
                                                                                                                                           cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTg          >  1:27142/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:822290/150‑1 (MQ=255)
                                                                                                                                              aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGta       <  1:336978/150‑1 (MQ=255)
                                                                                                                                                   cccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCATCTTTGTAGATCACGTGCTGGCGAACAAATGGATCGAATTTTTTCAGTTCGAGTATTTCAGGCTTAGTCCGTTTGTTCTTCTTAGTGGTATAGaa  >  1:973359/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATAT‑GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAA  >  NC_000913/3811065‑3811360

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