Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913826,9540AT39.7% 104.8 / 29.1 69I100N (ATT→AAT) ybhSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/5);  new base T (0/27);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 2.96e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTG  >  NC_000913/826814‑827097
                                                                                                                                               |                                                                                                                                                
gccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTcc                                                                                                                                            <  1:476073/147‑1 (MQ=255)
         cTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtttt                                                                                                                              >  1:823946/1‑151 (MQ=255)
              aCATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCg                                                                                                                          >  1:116422/1‑150 (MQ=255)
               cATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGAt                                                                                                                        >  1:554903/1‑151 (MQ=255)
                aTACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATc                                                                                                                       >  2:1091641/1‑151 (MQ=255)
                 tcccccTGTAAAAGGTTAGCGGTTTTCGGCTCACTGCGGTGGGTGATCCCCTGAATCGGTGCGGTGGCGTGGGCGCGCTCCATCTTTCCCGCAAAACCCCCCGGAATACCCACCAGACCGCGATTTTTCCCCGCCTGCATTTTGGCGATc                                                                                                                       <  1:847302/148‑1 (MQ=255)
                 tcccccTGTAAAAAGTATGCGGTTTCCGCCTCACCGCGGTCGTTGACCCCCTAAATCGGTGGGGTGGGGTGGGCGCGCTCCTTTTGTCCCGCAAAATCCCCCGGAATACCCCCCAGCCCGAAATTTTCCCCCGCCTGCATTTTGGCGAt                                                                                                                        <  1:897525/147‑1 (MQ=255)
                 tACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATc                                                                                                                       >  2:994168/1‑150 (MQ=255)
                      cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTcc                                                                                                                 >  1:1231397/1‑151 (MQ=255)
                          aaaaaGTTGGGGGTATTCGGCTCCCTCCGGTGGGTGATCCCCTGAATGGGTGGGGTGGCTTGGGCGCCCTCCTTCTGTCCCGAAAAACCCCCGGAAATACCCACCAGCCCGGGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCt                                                                                                                <  1:12679/146‑1 (MQ=255)
                          aaaaaGTTGGCGGTGTTCGGCTCCCTGGCGGTGATCATCCCCTGAATCGGTGCGGTGGTGTGGGCCCGCTCCATCGGTCCCCCAAAACCCCCCGCAATACCCTCCAGGCCGCAAGTTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTg                                                                                                               <  1:619449/147‑1 (MQ=255)
                                ttAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTAt                                                                                                       >  1:549841/1‑151 (MQ=255)
                                     ggTATTCGCCTCTGCCCCGTCGTTGACCCCCTAAACCGGTGGGGGGGCTTGGGCGCCCTCCTTCTTTTCCACAAACCCCCCCAAAATACCCCCCGGCCCGCTATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGc                                                                                                    <  1:628266/149‑1 (MQ=255)
                                       tATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAt                                                                                                >  2:317902/1‑151 (MQ=255)
                                           cggcccccgcccGGCGGTGACCCCCTGAATCGGGGGGGTGGCGTGGGCGCGCCCCTTCTGTCCCAAAAAACCCACGGAAATACCCACCAGACCGCGATTTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtg                                                                                            <  1:776991/141‑1 (MQ=255)
                                              cTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCg                                                                                         >  1:954466/1‑151 (MQ=255)
                                                      cGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:1122734/1‑150 (MQ=255)
                                                      cGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  2:1174745/1‑150 (MQ=255)
                                                       gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:256360/1‑149 (MQ=255)
                                                       gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                                  >  1:1316187/1‑149 (MQ=255)
                                                         cGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAggg                                                                              >  1:410720/1‑151 (MQ=255)
                                                            tGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCATCGATGTAGGGCg                                                                            >  2:798762/1‑150 (MQ=255)
                                                             gATCCCCTGAATGGGGGCGGTGGCTTGGGCCCGCTCTTTTGTTCCCCCAAAATCCCCGGAAATAAACCCCGGGCAGCAAATTTGCCCCGCCTGCATTTTGGCGGTCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGa                                                                           <  2:688105/150‑1 (MQ=255)
                                                             gATCCCCGGAATCGGGGCGGTGGCGTTGGCGCGCTCCATCTGTCCCGAAAAACCCCCCGAAATAACCCCCAGCCCGCAAATTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGaa                                                                          <  1:420603/151‑1 (MQ=255)
                                                               tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAcc                                                                        >  1:1211475/1‑151 (MQ=255)
                                                                ccccTGATTCGGTGGGTGGGCGTGGCCGCGCTCAATCTGTTCCGCAAAACCCCACGGAAAACCCACAAGCCCGGGATTTTCCCCCGCCTGCATTTTGGCGAGCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGc                                                                             <  2:5938/143‑1 (MQ=255)
                                                                 aCCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCgg                                                                      >  1:840818/1‑151 (MQ=255)
                                                                  cctaacccggggcggggGGGTGGGCGCCCTCCATCTGTCCCGAAAAACCCCCCGAAATACCCCCCAGACCGAGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGaa                                                                          <  2:1231333/139‑1 (MQ=255)
                                                                    tGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTc                                                                    >  2:1168568/1‑150 (MQ=255)
                                                                    tGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTc                                                                    >  2:998534/1‑150 (MQ=255)
                                                                      actggttggtggggCGTTGGCGCGCCCCATCTTTCCCGCAAAATCCACCGGAATACCCCCCAGACCGCTATTTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                         <  1:89998/131‑1 (MQ=255)
                                                                      aaggggtggggggGCGTGGGCCCGCTCCATCTGTCCCAAAAACCCCCCCGGAATACCCCCCAGCCCGAGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGt                                                                     <  2:1101848/143‑1 (MQ=255)
                                                                             ggggtgggtggggCGCGCTCCTTTTGTCCCGCAAATCCCCCCGGATAAACCACCGCCCCCTGATTTCCCACCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAgggtg                                                                            <  2:250315/122‑3 (MQ=255)
                                                                             gcggtgGCGTTGGCGCCCTCCATCTTTCCCGAAAAACCCCCCGAAATACCCCCCAGCCCGCGATTTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGtgtg                                                            <  1:395551/149‑1 (MQ=255)
                                                                               ggtggtGTTGGCGCCCTCCATCTGTTCCAAAAAACCCACGGAAATACCCACCAGCCCGGGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTgg                                                           <  1:914265/148‑1 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAGAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTg                                                        >  2:1135618/1‑151 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGt                                                         >  2:936813/1‑150 (MQ=255)
                                                                               ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTg                                                        >  2:1216255/1‑151 (MQ=255)
                                                                                     gTTGGCGCGCTCTATCTGTACCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATc                                                   >  2:329466/1‑150 (MQ=255)
                                                                                        ggcccccTCCTTCTGTCCCAAAAAACCCCCCGAAATACCCCCCGGCCCGGAAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATc                                                   <  1:360621/141‑1 (MQ=255)
                                                                                        ggCGCGCCCCACCTTTCCCGCAAACGCCCCCGAAATAACCACCCGGCCGCTATTTTTCGCCGCCTGCATTTTGGCGATCAGTTCCTGGCGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCa                                                 <  2:453793/149‑1 (MQ=255)
                                                                                        gcctccccccctttGTCCCGAAAACCCCCCCGAAAAACCCCCCAGACCCCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAg                                                <  2:903098/137‑1 (MQ=255)
                                                                                            cgcTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGcgccg                                           >  2:538180/1‑151 (MQ=255)
                                                                                             gcTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGcgccgc                                          >  1:1114343/1‑151 (MQ=255)
                                                                                                ccATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTc                                       >  2:808812/1‑151 (MQ=255)
                                                                                                 cttCTGTCCCGAAAAACCCCCGGAAATACCCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCtt                                        <  1:856325/147‑1 (MQ=255)
                                                                                                 cATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTc                                       >  1:74747/1‑150 (MQ=255)
                                                                                                  aTCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                     >  2:1272684/1‑151 (MQ=255)
                                                                                                     tctttgcGCAAAGATCCCCCGTAATAACCACAAGCCGGCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAgc                                               <  1:284511/134‑1 (MQ=255)
                                                                                                              aaaaTCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTcc                          >  1:919098/1‑151 (MQ=255)
                                                                                                              aaaaTCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTc                           >  1:1238418/1‑150 (MQ=255)
                                                                                                               aaaTCCCCCGGGAATACCCACCACCCCGAGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                                <  2:1075575/145‑1 (MQ=255)
                                                                                                                     ccccagaaaaaccccccaGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTaa                     <  2:325540/136‑1 (MQ=255)
                                                                                                                     cACCGAAATAACCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTaaa                    <  2:116418/151‑1 (MQ=255)
                                                                                                                       ccGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAAt                  >  2:967460/1‑151 (MQ=255)
                                                                                                                       ccGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAAt                  >  1:908156/1‑151 (MQ=255)
                                                                                                                          gcaaTACCCCCCAGACCGCGATTTTTCTCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATcc                <  2:1008556/149‑1 (MQ=255)
                                                                                                                            aaTACCCCCCAGCCCGCATTTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATc                 <  2:948477/148‑1 (MQ=255)
                                                                                                                              tAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAc            >  2:956493/1‑151 (MQ=255)
                                                                                                                                 cccccaGCCCGGGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGa             <  1:371404/144‑1 (MQ=255)
                                                                                                                                   acccGACCGCGAATTTTCCCCGCAGGCATTTTTGCGATCAGTTCCTGACGGTTATCGCTGTTAGTGGTGACGATGTAGGTTGAACCGGTCATGGGGTGGGTGATATCCACCGCCGCTTCGCTACGCTGTTCCAGCAAACTCCCGAACCGc        >  2:1107989/1‑150 (MQ=255)
                                                                                                                                   accaGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGc        >  1:333840/1‑150 (MQ=255)
                                                                                                                                      aGACCGGGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa       <  1:1237019/148‑1 (MQ=255)
                                                                                                                                       gACCGCGATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAg      <  1:1179052/148‑1 (MQ=255)
                                                                                                                                       gACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCtt   >  2:1061086/1‑151 (MQ=255)
                                                                                                                                         ccGCTATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTg  <  2:1272626/150‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa       <  2:763210/145‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAg      <  2:1041259/146‑1 (MQ=255)
                                                                                                                                         ccGCAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCtt   <  1:70541/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCAGCTTG  >  NC_000913/826814‑827097

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