Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,788 T→C 6.1% intergenic (+20/+33) oppF → / ← yciU oligopeptide ABC transporter ATPase/UPF0263 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,7880TC 6.1% 200.3 / 2.4 65intergenic (+20/+33)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (32/33)
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.45e-01

TACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGAACAGCAGAAC  >  NC_000913/1306520‑1307048
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
tACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGACCTGGAGACGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAATCCGCTTTCCGGTTGTGTTTTCCGTACCTGTTGCCCGTTTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGCGGGGATTCTCAGACACTACGTTTCTTGCCCGAAGGTCGATCAGCTTTAAAAACAATAAGGGCTGACTGTTgccca                                                                                                                                                                                                                                                                >  4:213701/1‑271 (MQ=255)
       aTCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGCGTGCGCCAAAACTCGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGAGCCGCTTTAGAAACAATAAGGCCTGACAGTTATCAGCCCtttt                                                                                                                                                                                                                                                         >  2:363526/1‑275 (MQ=255)
                    ttACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCATCCTAAAAg                                                                                                                                                                                                                                             >  2:211118/1‑274 (MQ=255)
                      aaacCAGGGCATTGATGTCGTCAGTCCCCATACCTGAACCGGAACTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAAGTACCGTCGCCGATCAACCCCCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                                                                                                                                                                                          <  4:78926/273‑1 (MQ=255)
                               cATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATtatttatt                                                                                                                                                                                                                                 <  2:431300/275‑1 (MQ=255)
                                       tCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGCGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                                                                                                                                                                                            >  1:72519/1‑240 (MQ=255)
                                                          aTCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCACTTTTTCACGCTAAAAGCGCTTATTTATTCCCGCCAGATGATTTGGCAACGTTtgt                                                                                                                                                                                                      >  1:43615/1‑275 (MQ=255)
                                                               gATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCAGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                                                                                                                                                                                  >  2:93715/1‑274 (MQ=255)
                                                                    ggAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGTTTGTGATCTTTTTc                                                                                                                                                                                            >  2:68365/1‑275 (MQ=255)
                                                                     gagaAGAACAAAACCATCCAGTTACTGGAAGTGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTATGTTTTCCGTACCTGTTGCCCGATTGCCGGTCCAGAGTACGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcg                                                                                                                                                                                           <  2:213704/275‑1 (MQ=255)
                                                                         agaaCAAACCCAGCCAGTTACTGAAAGGGGAACTACCATCGCCGTTCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCGGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                                                                                                                                                                                       <  2:28918/275‑1 (MQ=255)
                                                                           aaCAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAACGTTTGTGATCTTTTTCGCGGCAt                                                                                                                                                                                     >  2:105932/1‑275 (MQ=255)
                                                                                     tCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGCCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc                                                                                                                                                                           >  2:399251/1‑275 (MQ=255)
                                                                                     tCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACCCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTACCGACAGATGATATGGCACAGTTTGTGATCTTTTTCGCGTCATAAACGAAtcg                                                                                                                                                                           >  4:233945/1‑273 (MQ=255)
                                                                                         gTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcg                                                                                                                                                                                           <  2:219449/255‑1 (MQ=255)
                                                                                            aCTGGAAGGGGAATTACCGTCACCGATCAACCCGCCTGCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGAGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaa                                                                                                                                                                    <  2:312734/275‑1 (MQ=255)
                                                                                                    gggAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGCTTATTTATTCCCGc                                                                                                                                                                                                                            >  2:535511/1‑211 (MQ=255)
                                                                                                                         aCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGtt                                                                                                                                                                                                                                                                     <  6:301086/149‑1 (MQ=255)
                                                                                                                          cccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                                                                                                                                                                                                    <  5:188239/149‑1 (MQ=255)
                                                                                                                                  ccGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                                                                                                                                                                                            >  7:168299/1‑149 (MQ=255)
                                                                                                                                    ggTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                                                                                                                                                                                          <  6:1197/149‑1 (MQ=255)
                                                                                                                                       tgtgTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTc                                                                                                                                                                                                                                                       <  7:315939/149‑1 (MQ=255)
                                                                                                                                           ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                                                                                                                                   >  8:187656/1‑149 (MQ=255)
                                                                                                                                             ttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTa                                                                                                                                                                                                                                                 >  7:116587/1‑149 (MQ=255)
                                                                                                                                                cGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCACACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTCTCAGTTGTCAGCCCTTTTTCACGCTaaaa                                                                                                                                                                                                                                              <  5:320662/149‑1 (MQ=255)
                                                                                                                                                     ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                                                                                                                                         >  6:464544/1‑149 (MQ=255)
                                                                                                                                                                  gCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtta                                                                                                                                                                                                                                                          >  5:401242/1‑118 (MQ=255)
                                                                                                                                                                  gCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtta                                                                                                                                                                                                                                                          <  6:401242/119‑2 (MQ=255)
                                                                                                                                                                    cGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGTGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                                                                                                                                                                                                  >  5:389527/1‑109 (MQ=255)
                                                                                                                                                                    cGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                                                                                                                                                                                                  <  6:389527/109‑1 (MQ=255)
                                                                                                                                                                       tacAGAGTTGGACAAAACAGGTCCTGTTGTGGAGGGGGGTTTTAAGCCATACGTTTGTTGTCTGAAATGCGATCCGATTTAAAAAGAATCAGGGCTGACAGTTGTCAGCAGTTTTTCACACGAAAAGCGTTTATTTAGTCCCGCCAGATAATATGGCAAAGGTAGTGCTCTTTTTCGCGGCATAACAGAATGCGGGCAAAAGCATCGTTGATTTCACCATCTTCACTATCCGCCAGACCAATCACCACTACGGCAAAAAAGTCCGGCTTCAAAt                                                                                          <  4:399246/272‑1 (MQ=255)
                                                                                                                                                                       tCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                                                                                                                                                                                       >  6:206975/1‑149 (MQ=255)
                                                                                                                                                                            agTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                                                                                                                                                                                                  <  5:206975/149‑1 (MQ=255)
                                                                                                                                                                                cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCaa                                                                                 >  4:149957/1‑274 (MQ=255)
                                                                                                                                                                                cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCaa                                                                                 >  2:181235/1‑274 (MQ=255)
                                                                                                                                                                                cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAAc                                                                                <  2:233949/275‑1 (MQ=255)
                                                                                                                                                                                cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGACAGATGATATGGCACAGTTTGTGATCTTTTTCGCGGCATAACAGATTGCAGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAAACCAATCACCACTTCGGCAAACAAGGCCGGTTTCAAATCGAAATCAAc                                                                                >  2:450798/1‑275 (MQ=255)
                                                                                                                                                                                      aaCACGTCGGGTTCTGAGGGTGAGGTTCAGACACTCCGTTTCTTGCCTGAAAGTGTCTCCACGTTAACAACAATAAGGGCTGACAGTTGTCAGCCCTGTTTAACTCTAAAGTGGATTATTGATTCCCGCCAGTTGATATTGCAAAGTTTGTGATCTGTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCGTCTTCACTGTCCTCCAGACCAATCACCGCTTCGGCAAAATAGTCCGGGTTCAAATCGAAATCAACATGtt                                                                           <  4:105931/274‑1 (MQ=255)
                                                                                                                                                                                          cGTCCTGTTCCGGTGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGGCCTTTTTCACGCTAAAAGCGTTTATTTATTCCCGCCCGATGATATGGCAAAGTTTGTGCTCTTTTTCGTGGCATAACAGAATGCGGGCAACAACATCTTCAGTTACCACCTCTTCACTTTCAGCAAGCCCAACACCCACCTGGTCAAACATGTCAGGGTTCACACGGAAAACAACATGTAccctc                                                                      >  4:142570/1‑272 (MQ=255)
                                                                                                                                                                                                 ttCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                                                                                                                                                                                             <  5:464544/149‑1 (MQ=255)
                                                                                                                                                                                                     ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTCTTTCACGCTACAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTTGCGGCATAACAGACTGCGGGCAACAACATCTTTGATTTCACCATCTTCACTGTCCGCCAGTCCAACCACCACTTCGGCCAAAAAGTACGGGTTCAAATCCAAACCAACATTTTCCTGCCatgtttgcg                                                            >  4:313388/1‑266 (MQ=255)
                                                                                                                                                                                                      gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATAGTTTCCCGCCAGATGATATGGCCAAGTTTGTGATCGTTTTCGCCGCATAACAGAATGCGGCCAAAAACATCGATGCTTTCACTATCTTCACTGTCCGCCAGACCACTCACAACTTCGCCAAACAAGTCCGGGTCTAAACCAGAATCAACCTTCTCCTGCCAGTCTTccgt                                                           >  4:30415/1‑273 (MQ=255)
                                                                                                                                                                                                      gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTGAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAg                                                                                                                                >  2:447229/1‑205 (MQ=255)
                                                                                                                                                                                                         gggAGGTTCAGACAACGCGCTTGTGGCATGAAAGGCGATCAGCTTTAAAACCAGAAAGGTCTGACGGTGGTAAGCCCTTTTTCCCGCTAAAAGCGATGATTGATCCCCGCCAGAGGGTAGGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGGTGATTGCACCATCTTCGCTGTCCGCCAGACCAATGACCACTGCGGCAAAAAAGTCCGGGTTCAAATCGAAATCACCATGTTCCTGCCAGTCTTCCGCCGGa                                                       <  4:181233/275‑1 (MQ=255)
                                                                                                                                                                                                          ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTGAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAACGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCa                                                                                                                            <  4:447224/205‑1 (MQ=255)
                                                                                                                                                                                                           gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTNCGCGGCATAACAGAATGCGGGCAAAAACATCGGTGATTTCACCATCTTCACTGTCACCCAGACCCATCACCACCNCCGCACAAAAGTCCGGGTTCAAATAGAACANNAAAAGTACCTGCCCGACATCCGCCGGa                                                       >  3:62316/1‑273 (MQ=255)
                                                                                                                                                                                                                     cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                                                                                                                                         <  8:306116/149‑1 (MQ=255)
                                                                                                                                                                                                                        tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa                                                                                                                                                                      >  7:362283/1‑149 (MQ=255)
                                                                                                                                                                                                                        tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCcgcg                                        >  2:293662/1‑275 (MQ=255)
                                                                                                                                                                                                                                ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATc                                                                                                                                                              <  8:116587/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGct                              >  2:228254/1‑275 (MQ=255)
                                                                                                                                                                                                                                  gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCAAACTCACCATGCTCCTGCCAGTCTTCCGCCGGAACGACTAACTCCGCCCCGCCGCGct                              >  2:404525/1‑275 (MQ=255)
                                                                                                                                                                                                                                    cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                                                                                                                                                                                          <  5:53494/69‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                                                                                                                                                                                          >  6:53494/1‑69 (MQ=255)
                                                                                                                                                                                                                                    cTGAAAGGCGATCCGCTTTAAAAATGATAAGGGCTGACAGTTGTCATCCCGTTTTCACTATAAAAGGGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTGTTTCGCGGCATAACAGAATGCGAGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAGTCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGctct                            <  4:68364/275‑1 (MQ=255)
                                                                                                                                                                                                                                         gatCGATCCACTTTAAAAACAATAAGGGCTGTCAGTTGTCAGCCCTTTTTCCCGCTAAAAGCGATTATTTATTGCCGCCAGAGTCTATGGCAGAGTTTGTGATGATTTTCGCGGCATAACAGAGTGCGGGCAAAAACATCGTTGATTTCACCAACTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAGAAGTCCGGGTACAAATCGAAATCAACATGTGCCTGCCAGTCTTCGGCCGGGTCGAATAACTCCGCGCCGCCGCGCTCTTCaaa                       <  4:93714/273‑1 (MQ=255)
                                                                                                                                                                                                                                           tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACt                     <  2:149958/275‑1 (MQ=255)
                                                                                                                                                                                                                                                ccGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCt                                                                                                                                              <  8:168299/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGGGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGGTGATTTCACCAGCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGaa           <  2:142571/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGaa           <  4:293658/274‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaaaaaGTGCTGACATGTGTCAGACCTTTTTACCGCTACAAGCGATTATTTATTCCCGCCATATGATATGGCAAAGTTTGTGATCTTTTTCGCTGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGAAcagcag     <  4:228250/270‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGtt                                                                                                                                                           >  8:434517/3‑122 (MQ=255)
                                                                                                                                                                                                                                                              aaaaGGGCTGACAACTGACAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGAGCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGtt                                                                                                                                                           <  7:434517/120‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTAAGGGCTGTCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGAACAGCAGAAc  <  1:62316/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCGTCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGAACAGCAGAAc  <  2:30415/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATc                                                                                                                         <  7:187656/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
TACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGATCGAATAACTCCGCGCCGCCGCGCTCTTCAAACTGAAGATTGAACAGCAGAAC  >  NC_000913/1306520‑1307048

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