Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,788 T→C 5.4% 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 5.4% 808.4 / 23.2 255intergenic (+20/+33)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base T (122/119);  new base C (9/5);  total (131/124)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 263 total aligned reads displayed.

AACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCAC  >  NC_000913/1306640‑1306934
                                                                                                                                                    |                                                                                                                                                  
aaCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGt                                                                                                                                                    >  1:1378599/1‑149 (MQ=255)
 aCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGtt                                                                                                                                                   <  3:1976475/149‑1 (MQ=255)
 aCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGtt                                                                                                                                                   <  3:541098/149‑1 (MQ=255)
   ccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGt                                                                                                                                                 <  4:1143916/149‑1 (MQ=255)
    cGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                                                                                <  1:1907287/149‑1 (MQ=255)
     gCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCa                                                                                                                                               <  4:1079834/149‑1 (MQ=255)
     gCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCa                                                                                                                                               >  1:1524408/1‑149 (MQ=255)
       cTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                                                                             <  1:397620/149‑1 (MQ=255)
        ttCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                                                                            >  4:2425870/1‑149 (MQ=255)
         tCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                                                                           <  4:2186508/149‑1 (MQ=255)
            ggTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                                                                        <  2:2203571/149‑1 (MQ=255)
             gTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtttt                                                                                                                                       <  2:1197624/149‑1 (MQ=255)
              ttGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttttt                                                                                                                                      >  4:1517551/1‑149 (MQ=255)
              ttGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttatt                                                                                                                                      <  3:1907533/149‑1 (MQ=255)
                gtgtTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                                                                    <  1:1956625/149‑1 (MQ=255)
                 tgtTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattgtt                                                                                                                                   <  4:675665/149‑4 (MQ=255)
                 tgtTTTCCGTACCCGTTGCCCGATTGCCGGTACAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                                                                   <  2:2356524/149‑1 (MQ=255)
                  gtTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                                                                                  <  3:701400/149‑1 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                 >  1:1793923/1‑149 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                 >  1:1196775/1‑149 (MQ=255)
                        cGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaaa                                                                                                                            <  2:1263759/149‑1 (MQ=255)
                         gTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAg                                                                                                                           >  2:2449968/1‑149 (MQ=255)
                         gTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAg                                                                                                                           >  1:2372096/1‑149 (MQ=255)
                         gTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAg                                                                                                                           >  3:1056211/1‑149 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                       >  3:538736/1‑149 (MQ=255)
                                ttGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattat                                                                                                                    <  4:1054452/149‑1 (MQ=255)
                                  gCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                                                                    <  2:1458383/131‑1 (MQ=255)
                                     cGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATtatttatt                                                                                                               <  2:632712/149‑1 (MQ=255)
                                     cGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATtatttatt                                                                                                               <  4:215910/149‑1 (MQ=255)
                                                 cAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatga                                                                                                   >  3:353262/1‑149 (MQ=255)
                                                 cAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatga                                                                                                   >  2:2411976/1‑149 (MQ=255)
                                                   gagTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGata                                                                                                 <  3:579683/149‑1 (MQ=255)
                                                         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                                                                           >  1:2009064/1‑149 (MQ=255)
                                                         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                                                                           >  1:1267518/1‑149 (MQ=255)
                                                         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                                                                           >  4:760992/1‑149 (MQ=255)
                                                           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCAGCTTTAAAAATAATAAGGGATGACAGTTGTCAGCCCTTTTTCACTCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCagga                                                                                         >  2:1721677/1‑146 (MQ=255)
                                                            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGt                                                                                        >  1:2077160/1‑149 (MQ=255)
                                                                  cGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                                                                  <  3:1006988/149‑1 (MQ=255)
                                                                      cTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                                                                              >  3:1513041/1‑149 (MQ=255)
                                                                       tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttt                                                                             >  3:168591/1‑149 (MQ=255)
                                                                       tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttt                                                                             >  2:842075/1‑149 (MQ=255)
                                                                       tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttt                                                                             <  4:901661/149‑1 (MQ=255)
                                                                       tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttt                                                                             <  1:2149079/149‑1 (MQ=255)
                                                                         ttCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                                                                           >  1:475320/1‑149 (MQ=255)
                                                                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                                                                      >  1:2294995/1‑149 (MQ=255)
                                                                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                                                                      >  3:2157129/1‑149 (MQ=255)
                                                                                 gggAGTTTCAGACACTCCGTTTCTTGCCTGAGAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                                                                  <  2:1917891/134‑1 (MQ=255)
                                                                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa                                                               <  3:760992/149‑1 (MQ=255)
                                                                                      tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg                                                              >  2:2014516/1‑149 (MQ=255)
                                                                                            acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg                                                        <  4:1586904/149‑1 (MQ=255)
                                                                                            acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg                                                        <  1:2122083/149‑1 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                       <  4:799230/149‑1 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                       <  1:1592393/149‑1 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                       >  1:291637/1‑149 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                       >  2:2169605/1‑149 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                       <  4:1287099/149‑1 (MQ=255)
                                                                                              acTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTTAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg                                                      >  2:1673820/1‑149 (MQ=255)
                                                                                              acTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCATATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg                                                      >  3:1219471/1‑149 (MQ=255)
                                                                                              acTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg                                                      >  1:725813/1‑149 (MQ=255)
                                                                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa                                                    >  4:330830/1‑149 (MQ=255)
                                                                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa                                                    >  3:320315/1‑149 (MQ=255)
                                                                                                 ccGTTTCTTGCCTGAAAGTAGATGCGCTTTAAAAACAATAAGGCCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAATATTTATTCCCGCCAGATGATATGGCAAAGTCTGTGATCTTTTTCGCGGCATAACAGTATGCGGGCaa                                                   >  2:177981/1‑149 (MQ=255)
                                                                                                   gTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                                                                                <  3:2392806/102‑1 (MQ=255)
                                                                                                    tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaaaa                                                >  4:1916677/1‑149 (MQ=255)
                                                                                                       cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACAt                                             <  2:123135/149‑1 (MQ=255)
                                                                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg                                                        <  3:1790569/137‑1 (MQ=255)
                                                                                                         tGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                                                                       >  3:2457237/1‑121 (MQ=255)
                                                                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGtt                                         <  1:986935/149‑1 (MQ=255)
                                                                                                            cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                                                                  <  2:1122760/107‑1 (MQ=255)
                                                                                                            cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTg                                        >  2:1830294/1‑149 (MQ=255)
                                                                                                            cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCATTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTg                                        <  3:862842/149‑1 (MQ=255)
                                                                                                              gAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGAt                                      <  3:2077637/149‑1 (MQ=255)
                                                                                                              gAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGAt                                      >  3:619708/1‑149 (MQ=255)
                                                                                                                aaGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGAttt                                    >  3:901820/1‑149 (MQ=255)
                                                                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAACGCGATTATTCATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCa                                  >  1:2512042/1‑149 (MQ=255)
                                                                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCa                                  <  3:2277692/149‑1 (MQ=255)
                                                                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCa                                  >  3:1836653/1‑149 (MQ=255)
                                                                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCa                               <  4:2271787/149‑1 (MQ=255)
                                                                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCAt                              >  1:806529/1‑149 (MQ=255)
                                                                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCAt                              >  1:1814697/1‑149 (MQ=255)
                                                                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACCTCGTTGATTTCACCATCtt                           <  3:1458147/149‑1 (MQ=255)
                                                                                                                          gCTTTAAAAACAATAAGGGCTAACAGTTGTCTGCCCTTTTACACGCTAAAAGGGATTATTTATTCGCCCCAGATGATATGGCAAAGTTTGTGATCTTGTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTTATTTCACCATCTTc                          <  4:1648145/149‑1 (MQ=255)
                                                                                                                           cTTTAAAAACAAGAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCAGCTTCa                         <  2:542137/149‑1 (MQ=255)
                                                                                                                            tttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCAc                        >  3:448721/1‑149 (MQ=255)
                                                                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCATATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCAc                                 <  4:49723/136‑1 (MQ=255)
                                                                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTc                    >  2:414245/1‑149 (MQ=255)
                                                                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTcc                   >  1:1044337/1‑149 (MQ=255)
                                                                                                                                 aaaCAATAAGGGCTGACAGTCGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTcc                   >  2:1340332/1‑149 (MQ=255)
                                                                                                                                  aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa                                                               <  1:1105777/104‑1 (MQ=255)
                                                                                                                                  aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCg                  >  3:219168/1‑149 (MQ=255)
                                                                                                                                  aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCg                  >  1:1127649/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCa               >  2:291309/4‑149 (MQ=255)
                                                                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACa                                              <  3:125874/114‑1 (MQ=255)
                                                                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACa                                              >  4:125874/1‑114 (MQ=255)
                                                                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGAcc           >  1:1405926/1‑149 (MQ=255)
                                                                                                                                          gggCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCa          >  4:1810900/1‑149 (MQ=255)
                                                                                                                                           ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGt                                          >  2:1542871/1‑116 (MQ=255)
                                                                                                                                           ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCaa         >  4:516826/1‑149 (MQ=255)
                                                                                                                                           ggCGGACAGTTGTCAGCCCTTTTTCACGCTAAGAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCaa         >  4:250031/1‑149 (MQ=255)
                                                                                                                                                  aGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATcaccac  >  1:1638528/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCAC  >  NC_000913/1306640‑1306934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: