Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,167 A→C 5.9% intergenic (‑225/+59) 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,1670AC 5.9% 260.2 / 6.0 85intergenic (‑225/+59)pheS/pheMphenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (44/35);  new base C (3/2);  total (48/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 = 9.76e-01

CATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAA  >  NC_000913/1799018‑1799314
                                                                                                                                                      |                                                                                                                                                   
cATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCa                                                                                                                                                     >  1:1405665/1‑150 (MQ=255)
      ggTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGa                                                                                                                                               >  2:425055/1‑150 (MQ=255)
       gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                              >  2:532650/1‑150 (MQ=255)
       gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                              >  2:1047694/1‑150 (MQ=255)
          ccGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCt                                                                                                                                           >  2:1187456/1‑150 (MQ=255)
             ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                        <  2:1359662/150‑1 (MQ=255)
             ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                        <  1:941707/150‑1 (MQ=255)
             ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                        <  1:516386/150‑1 (MQ=255)
             ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                        <  1:469482/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:824611/1‑150 (MQ=255)
                                       tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                              >  2:781092/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:318842/150‑1 (MQ=255)
                                                       tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                              <  2:1401134/150‑1 (MQ=255)
                                                       tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                              >  1:438005/1‑150 (MQ=255)
                                                       tctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTg                                                                                              <  1:467899/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                                                                >  2:253638/1‑150 (MQ=255)
                                                                                     cGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCg                                                                >  1:1209560/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:592843/1‑150 (MQ=255)
                                                                                              cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                       >  2:649040/1‑150 (MQ=255)
                                                                                              cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                       <  1:532650/150‑1 (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:452995/1‑150 (MQ=255)
                                                                                                             gCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgcttg                                        >  2:870811/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      >  1:1324087/1‑150 (MQ=255)
                                                                                                                                                gCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt      >  2:155704/1‑150 (MQ=255)
                                                                                                                                                 ccTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAtt     >  1:822566/1‑150 (MQ=255)
                                                                                                                                                  cTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg    >  1:1210195/1‑150 (MQ=255)
                                                                                                                                                  cTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg    <  2:1393973/150‑1 (MQ=255)
                                                                                                                                                    ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa  <  1:1402803/150‑1 (MQ=255)
                                                                                                                                                    ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa  >  2:475584/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAA  >  NC_000913/1799018‑1799314

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: