Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,169 T→G 5.8% intergenic (‑227/+57) pheS ← / ← pheM phenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,1690TG 5.8% 278.2 / 7.8 86intergenic (‑227/+57)pheS/pheMphenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base G (3/2);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAG  >  NC_000913/1799024‑1799317
                                                                                                                                                  |                                                                                                                                                    
ggTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGa                                                                                                                                                  >  2:425055/1‑150 (MQ=255)
 gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                                 >  2:1047694/1‑150 (MQ=255)
 gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                                 >  2:532650/1‑150 (MQ=255)
    ccGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCt                                                                                                                                              >  2:1187456/1‑150 (MQ=255)
       ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                           <  1:469482/150‑1 (MQ=255)
       ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                           <  1:516386/150‑1 (MQ=255)
       ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                           <  1:941707/150‑1 (MQ=255)
       ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                           <  2:1359662/150‑1 (MQ=255)
         cTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCCAGTGGAGGCtttttt                                                                                                                                         <  2:1162423/150‑4 (MQ=38)
         cTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAg                                                                                                                                         <  2:194326/150‑1 (MQ=255)
          tgtgCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTATTCCTGTTGAGAGAGGCTTTTTATTTTTCAAAACGCGCATAGAAAAAAAGTCTCCCAGTGGAGGCttttttt                                                                                                                                        <  2:298964/150‑5 (MQ=25)
              cTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCt                                                                                                                                    >  2:830276/1‑150 (MQ=255)
                cGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGt                                                                                                                                  >  2:318842/1‑150 (MQ=255)
                     ttCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTcc                                                                                                                             <  1:914578/150‑1 (MQ=255)
                      tctcACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCg                                                                                                                            <  1:631086/150‑1 (MQ=255)
                        tcACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGtt                                                                                                                          <  1:1284934/150‑1 (MQ=255)
                             tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAACCGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGGttctct                                                                                                                     >  2:994809/1‑150 (MQ=255)
                             tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctct                                                                                                                     <  1:1112145/150‑1 (MQ=255)
                             tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctct                                                                                                                     <  1:364648/150‑1 (MQ=255)
                             tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACATGATGAGAAAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctct                                                                                                                     >  2:999720/1‑150 (MQ=255)
                              aCTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctctt                                                                                                                    <  2:934983/150‑1 (MQ=255)
                               cTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttc                                                                                                                   <  2:1342708/150‑1 (MQ=255)
                                 tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                                 >  2:781092/1‑150 (MQ=255)
                                 tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                                 >  2:824611/1‑150 (MQ=255)
                                 tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                                 <  1:1125480/150‑1 (MQ=255)
                                  aaaTACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCa                                                                                                                >  2:1147690/1‑150 (MQ=255)
                                       cGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACgcgc                                                                                                           >  1:869386/1‑150 (MQ=255)
                                         cTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTa                                                                                                         >  1:376653/1‑150 (MQ=255)
                                          tCCGGTTTCTCCGCGCTGGGCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTGAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAg                                                                                                        <  1:830276/150‑1 (MQ=255)
                                          tCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAg                                                                                                        <  1:1187456/150‑1 (MQ=255)
                                           ccGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGc                                                                                                       >  1:1001965/1‑150 (MQ=255)
                                               tttCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTATTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCGCTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctcc                                                                                                   >  2:145159/1‑150 (MQ=255)
                                                 tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                                 <  1:467899/150‑1 (MQ=255)
                                                 tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                                 >  1:438005/1‑150 (MQ=255)
                                                 tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                                 <  2:1401134/150‑1 (MQ=255)
                                                 tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                                 <  1:318842/150‑1 (MQ=255)
                                                  ctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTgg                                                                                                >  2:50967/1‑150 (MQ=255)
                                                  ctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTgg                                                                                                >  2:452032/1‑150 (MQ=255)
                                                      gcgcTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTc                                                                                            <  1:1322677/150‑1 (MQ=255)
                                                       cgcTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCa                                                                                           <  1:425055/150‑1 (MQ=255)
                                                          tGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGt                                                                                        >  2:839414/1‑150 (MQ=255)
                                                            gCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGc                                                                                      <  1:442116/150‑1 (MQ=255)
                                                            gCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGc                                                                                      >  1:482495/1‑150 (MQ=255)
                                                                 gtCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAg                                                                                 <  1:351004/150‑1 (MQ=255)
                                                                 gtCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAg                                                                                 >  2:367803/1‑150 (MQ=255)
                                                                         tGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaaaaaaa                                                                         >  1:395955/1‑150 (MQ=255)
                                                                          ggCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaag                                                                        <  2:46231/150‑1 (MQ=255)
                                                                           gCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaaga                                                                       <  2:991170/150‑1 (MQ=255)
                                                                            cTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaagaa                                                                      >  2:550991/1‑150 (MQ=255)
                                                                             tGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaagaag                                                                     >  1:163464/1‑150 (MQ=255)
                                                                              gCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGc                                                                    >  1:822608/1‑150 (MQ=255)
                                                                               cGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCg                                                                   >  1:1209560/1‑150 (MQ=255)
                                                                               cGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCg                                                                   >  2:253638/1‑150 (MQ=255)
                                                                                gACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCgg                                                                  >  2:608365/1‑150 (MQ=255)
                                                                                   aaTTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACAGCGCCT‑‑‑‑‑GAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATag                                                           >  2:705508/1‑150 (MQ=255)
                                                                                   aaTTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGaaa                                                               <  1:279784/150‑1 (MQ=255)
                                                                                      tACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATa                                                            <  1:50967/150‑1 (MQ=255)
                                                                                      tACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATa                                                            <  2:1405665/150‑1 (MQ=255)
                                                                                        cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                          >  2:649040/1‑150 (MQ=255)
                                                                                        cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                          <  1:532650/150‑1 (MQ=255)
                                                                                        cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                          >  2:592843/1‑150 (MQ=255)
                                                                                          ccTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGCCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagcag                                                        <  1:1020668/150‑1 (MQ=255)
                                                                                            tgatgaGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAgcagca                                                      >  1:1351570/1‑150 (MQ=255)
                                                                                                    cAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgc                                              <  1:1047694/150‑1 (MQ=255)
                                                                                                       gCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgcttg                                           >  2:870811/1‑150 (MQ=255)
                                                                                                       gCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgcttg                                           >  2:452995/1‑150 (MQ=255)
                                                                                                       gCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgcttg                                           <  1:1228754/150‑1 (MQ=255)
                                                                                                         tttttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCg                                         <  2:354944/150‑1 (MQ=255)
                                                                                                                 tttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCttt                                 >  1:1018204/1‑150 (MQ=255)
                                                                                                                    cAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTgg                              >  1:1393973/1‑150 (MQ=255)
                                                                                                                            cgcATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCtt                       >  2:79499/1‑150 (MQ=255)
                                                                                                                               aTACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCCCGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaa                    >  1:642739/1‑150 (MQ=255)
                                                                                                                                tACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaaa                   >  2:1238846/1‑150 (MQ=255)
                                                                                                                                  cAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAg                 <  1:824611/150‑1 (MQ=255)
                                                                                                                                  cAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAg                 >  2:228750/1‑150 (MQ=255)
                                                                                                                                   cctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGa                >  2:1318550/5‑150 (MQ=255)
                                                                                                                                    ctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAc               >  2:388758/4‑150 (MQ=255)
                                                                                                                                     tgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAcc              >  1:1254773/3‑150 (MQ=255)
                                                                                                                                          gCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt         >  2:155704/1‑150 (MQ=255)
                                                                                                                                          gCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt         >  1:1324087/1‑150 (MQ=255)
                                                                                                                                           ccTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAtt        >  1:822566/1‑150 (MQ=255)
                                                                                                                                            cTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg       <  2:1393973/150‑1 (MQ=255)
                                                                                                                                            cTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg       >  1:1210195/1‑150 (MQ=255)
                                                                                                                                              ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa     <  1:1402803/150‑1 (MQ=255)
                                                                                                                                              ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa     >  2:475584/1‑150 (MQ=255)
                                                                                                                                                 cTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAg  >  2:825901/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAG  >  NC_000913/1799024‑1799317

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