Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 10.1% intergenic (+131/+14) rraB → / ← yjgM protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0200AT10.1% 141.3 / 5.0 49intergenic (+131/+14)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/17);  new base T (2/3);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTGAATGCCGATCTGATCGATGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGGCGCTATCCCACCG  >  NC_000913/4478725‑4479316
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
tggagttcagacgtgtgctcttccgatcttGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                                                                                                                                                                                                                                                                                     <  2:546313/273‑1 (MQ=255)
      gCCGATCTGATCGTTGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                                                                                                                                                                                                                               <  1:188992/301‑2 (MQ=255)
      gCCGATCTGATCGATGTCCAGGTTGTACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACGTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                                                                                                                                                                                                                               <  2:276901/301‑2 (MQ=255)
                             tGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggagatcggaagagcgtcgtgtagggaaaga                                                                                                                                                                                                                                                                        >  1:545640/1‑272 (MQ=255)
                                     tgattaCGCTGTCAGATAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGGATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCACACTTcac                                                                                                                                                                                                                                                                <  1:214644/301‑1 (MQ=255)
                                                      aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                                                                                                                                                                               >  2:298471/1‑301 (MQ=255)
                                                                    gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTAGCGCAACATCCACACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                                                                                                                                                                                                 >  2:164453/1‑301 (MQ=255)
                                                                                   cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGGCGATATGCTCaaa                                                                                                                                                                                                                  >  1:287877/1‑301 (MQ=255)
                                                                                       aagttggAAGATCCTAAGGGTGAAGATGGGGAGGTGGAGTGTTTGGCGGATGAGGAGGAGGACGAGGTCGGCCACTATGAAAGTGTGGCGAACAGGGCCAAGGTAGTTCCGGCAACTGGGGAAATTTGGAGCCCGGTTAAGGCATTCAGGGAACATCCGGCATAACAAAAACGTAGGTGGTTAAGGTGCGAAGGCCCGCGAGTGGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGccc                                                                                                                                                                                                              <  1:105547/295‑1 (MQ=255)
                                                                                             gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACCTCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCCCAAAGCCCAAATGc                                                                                                                                                                                                        >  1:497307/1‑301 (MQ=255)
                                                                                                gATCCTTACGGCGAAGATGGTGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCAACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTGGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                                                                                                                                                                                                     <  1:6794/301‑1 (MQ=255)
                                                                                                         ggCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGATGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCAAAGCCCAAAAGCACAAAAAACGCa                                                                                                                                                                                            >  1:393613/1‑301 (MQ=255)
                                                                                                           cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCAAAGCCCAAATGCTCATAAAGCGCaaa                                                                                                                                                                                          >  2:565901/1‑300 (MQ=255)
                                                                                                           cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAACGCATAGACGAAATGCTCACAGCCCAAATGCCAAAAAAGCGCAAt                                                                                                                                                                                          >  1:122049/1‑301 (MQ=255)
                                                                                                                                         ttGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGATGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCa                                                                                                                                                             <  2:394127/300‑1 (MQ=255)
                                                                                                                                           gAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGACCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGa                                                                                                                                                           >  2:553738/1‑301 (MQ=255)
                                                                                                                                                            tCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAAATAGCAGCGTTTGAAAcc                                                                                                                                          >  1:224256/1‑301 (MQ=255)
                                                                                                                                                                    aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCAtctctc                                                                                                                                  >  2:110533/1‑301 (MQ=255)
                                                                                                                                                                      atAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTcgc                                                                                                                                <  1:110409/300‑1 (MQ=255)
                                                                                                                                                                       tAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTcgcg                                                                                                                               >  2:475210/1‑301 (MQ=255)
                                                                                                                                                                        ctttcccttcacgacgctcttccgatctCCTTCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCgcgc                                                                                                                              <  1:409621/274‑1 (MQ=255)
                                                                                                                                                                         aGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGcc                                                                                                                             <  2:224499/301‑1 (MQ=255)
                                                                                                                                                                                     tcagtcgtgtgctcttccgatctcTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGcc                                                                                                                 <  2:55078/278‑1 (MQ=255)
                                                                                                                                                                                         cccAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTATCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCCCCATCGCCCttt                                                                                                             >  1:635895/1‑300 (MQ=255)
                                                                                                                                                                                                aGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAg                                                                                                      <  2:497954/301‑1 (MQ=255)
                                                                                                                                                                                                 gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAga                                                                                                     >  2:211512/1‑300 (MQ=255)
                                                                                                                                                                                                    ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCgcgcagatcggaagagcacacgtctgaactcc                                                                                                  >  2:410184/1‑273 (MQ=255)
                                                                                                                                                                                                            cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGATCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCACCATCGCCagaacagaagagcctcgtgtagg                                                                                          >  1:55010/1‑279 (MQ=255)
                                                                                                                                                                                                               ttGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGccc                                                                                       <  2:122173/300‑1 (MQ=255)
                                                                                                                                                                                                                                   aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                                                   >  2:604881/1‑301 (MQ=255)
                                                                                                                                                                                                                                   aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTACTCCATCGCCCATAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgac                                                                   >  1:583606/1‑300 (MQ=255)
                                                                                                                                                                                                                                   aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATACTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                                                   >  2:630918/1‑301 (MQ=255)
                                                                                                                                                                                                                                     gcatgcatGTCACATCCGGCATAAACAAAACGTACTTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaaga                                                                 <  1:565187/301‑1 (MQ=255)
                                                                                                                                                                                                                                               cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGAAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTAGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCAAGCCCTTTGCCACGGACAGCCGACAGAAAATACAACt                                                       >  1:160486/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  aTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACAACtttt                                                    >  1:204064/1‑301 (MQ=255)
                                                                                                                                                                                                                                                    ccgtttgctcttccgatctcgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGc                                                  <  2:600709/282‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTAAACAAAACGGCCGTTTTTAATTTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATACTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTc                                             <  1:630041/301‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGTAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCAGACAGAAAAAACAACttttcaagtacacaa                                         >  1:441645/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                    ttccgtttgttttaaaaaaaTGCCCGCGTGTTCAACTCGCGGTCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCg                                  <  1:113329/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGAAAGCCGGAAGAAAATAAAACTTTTGCAGATCGGaagaacgtcgtg                               >  1:599948/1‑291 (MQ=255)
                                                                                                                                                                                                                                                                         ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGactca                             >  2:553951/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             ctctaaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGagatcggaa                         >  2:113453/10‑293 (MQ=255)
                                                                                                                                                                                                                                                                                tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCgg                      <  2:160637/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATACCCGGAAGAACACCCAACTTTTGCAGTTCGAACATAcccacaccactccccttt                    >  1:11543/5‑286 (MQ=255)
                                                                                                                                                                                                                                                                                   tCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGCCGTTTCCCGCTAGCAGCGTTTGAAACCCATCTCTCGAGCCTGCTCCATCGCCATTAAGGCCAGATTTTTTGCCAGCCCTTTGCCGCGGATCGCCGGAAGAAAAGACACCTATTGCCGTTCACAAATAACCTACTCACCCCCGGttt                   >  1:484433/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                    cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATAACAGACTCACTCCCGGTTaa                  >  2:332467/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGgcg             <  2:204307/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    gAACACGCGTGTATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGTAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGGCGCTATCCCAccg  <  2:442209/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    gAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATAACCGACTCACTCCCGGTTAATGGCGCTAACCCAccc  >  1:277721/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
TTGAATGCCGATCTGATCGATGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGGCGCTATCCCACCG  >  NC_000913/4478725‑4479316

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: