Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9940GT25.0% 229.1 / 12.0 163intergenic (+392/+631)yhaC/garKpentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (20/91);  new base T (1/36);  total (21/142)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

AAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTT  >  NC_000913/3269846‑3270134
                                                                                                                                                    |                                                                                                                                             
catctcgccgTACCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGTACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaa                                                                                                                                             <  1:538904/146‑1 (MQ=11)
aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                                             <  3:1466320/151‑1 (MQ=255)
aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                                             <  3:1373066/151‑1 (MQ=255)
aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                                             <  6:973595/151‑1 (MQ=255)
 aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                                            <  2:535648/151‑1 (MQ=255)
  aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                                           <  3:664843/151‑1 (MQ=255)
   cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                                          <  1:1487688/151‑1 (MQ=255)
   cccGCCGAATCGGTTTATTACGTACATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt                                                                                                                                          <  5:542988/151‑1 (MQ=255)
     cgcctGAACGGTTTATTACTTACATCAGGTAAAACTGACCGATAAGGCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACACGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                                                        <  1:754559/151‑1 (MQ=255)
     cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTTTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                                                                        <  2:436455/151‑1 (MQ=255)
     cgccgGAACGGTTTATTACGTAAATCAGGTAAAAATGACCCATAAGCAGCTTTCTATTGGGTATACTGTAGTGTACAGTCATTAAACATTAAGACCCTAAAAAAAAAAAAAAAACCCGAAGCGGGTTTTTAACAAAAACAGGTGaatatcc                                                                                                                                        >  5:1189204/1‑146 (MQ=255)
       ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGGCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                                      <  7:363816/151‑1 (MQ=255)
        cgtagcGGTTTATTACGTAAATCATGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCg                                                                                                                                     <  1:1022047/146‑1 (MQ=11)
        cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                                                                     <  6:496950/151‑1 (MQ=255)
        cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                                                                     <  3:609516/151‑1 (MQ=255)
        cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAAACCCCCGAAGCGGGGTTTTACGTAAAAAAGCTgcgaagtgccc                                                                                                                                     >  2:1139345/1‑141 (MQ=14)
          gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                                                   <  2:277404/151‑1 (MQ=255)
          gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                                                   <  5:356869/151‑1 (MQ=255)
          gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGAAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGaaagatcccaa                                                                                                                                    >  3:1412945/1‑142 (MQ=14)
            aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                                                 <  5:616378/151‑1 (MQ=255)
             cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                                                <  2:1308898/151‑1 (MQ=255)
             cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                                                <  5:372752/151‑1 (MQ=255)
             cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAg                                                                                                                                <  2:531028/151‑1 (MQ=255)
             cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTT‑GGTATAGTGTCGTGGACAGTAAATAATATTTATGCCCCGCCAAAAAAAAAAAACCCCCCAAGCCGGTTTTTAAGTAAAACAAGTGAAAAAGACCAAAAAGc                                                                                                                               >  7:572040/1‑151 (MQ=9)
              ggTTTATTACGTACATCAGTTAAAACTTACCGATAAGCGGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                                               <  6:1243478/151‑1 (MQ=255)
              ggTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                                               <  5:432225/151‑1 (MQ=255)
              ggTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                                               <  3:21383/151‑1 (MQ=255)
              ggTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                                               <  3:1344472/151‑1 (MQ=255)
                 ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGccccgttcggtgg                                                                                                                            >  2:937200/1‑141 (MQ=14)
                  tattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGCTTCATCTTTCTGGCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTTAATCAGGTGAAACGGACCgcaaacggggg                                                                                                                           >  2:851382/1‑141 (MQ=11)
                   attaCGTACAACAAGATACACTTAGCGATAAGCGGCATTCATTTGTCTATCGTGTCGCGGACAGTCACTCATATTAATGCCCCTCCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctt                                                                                                                          <  8:1293879/151‑1 (MQ=2)
                    ttaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGTGTACAGGGTCGTGGGCCCTCATGTCTCATATCGACGCACCAAAACGAAAAACCCCCCCCCGCGGGTTGTTTCATAAAACACGTGCAACGGGCCgggtcgggggggt                                                                                                                         >  8:1268395/1‑139 (MQ=2)
                      acttACATCAGGTAACACTGACCGATAAGCAGCTTTCTTTTGGGTATAGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc                                                                                                                        <  2:1091044/147‑1 (MQ=255)
                       cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGGGGACAGGCAATCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                                                      <  8:1075826/151‑1 (MQ=255)
                        gTAGATCAGGTAAAACTGTCCGATAAGCCGCTTTCTTTTGTTTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                                                     <  4:711689/151‑1 (MQ=255)
                         tACATAAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                                                     <  7:1387606/150‑1 (MQ=38)
                         gaCATCAGGTAAAACTGACCGAGCAGCCGCTGTCTTTTGCGTATCGTGTCGGGGCCAGGCCTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                                    <  8:1147311/150‑1 (MQ=14)
                          aCATCATGTATAAATGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGCGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                                    <  2:351894/150‑1 (MQ=11)
                          aCATCATGTAAAACTGACCGATAAGCGGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                                    <  8:118256/150‑1 (MQ=255)
                          aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGGGGACAGTCATTCATCTATCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                                    <  8:603084/150‑1 (MQ=255)
                          aCATCACTTAAAACTGTCCGATCAGCGGCATTCATTTGCCTGTAGTGACGTGGACAGTCTTTCCTCTATCTGGCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                                                   <  8:190321/151‑1 (MQ=11)
                           cATCAGGTAAAACTGTCCGTTAAGCCGCTTTCTTTTTGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                                  <  7:1480796/151‑1 (MQ=255)
                            atcatttAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTCTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                                                 <  1:216769/145‑1 (MQ=255)
                              cacttAATACTGACCGATAAGCCGCTTTCATTTGCCTATTGTGTCGTCGAAAGTGTTTCTTCTTTCTGCCACTCCAAAAGTTAAAACCCGCCGAATCGGGTTATTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                                               <  1:185171/147‑1 (MQ=11)
                                   aaaaCTGAGCGATAAGCAGCTTTCATTTGTTTATAGTGTCGTGGACAGTCATTCTTCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtg                                                                                                          <  4:378562/151‑1 (MQ=11)
                                   aaaaCTGACCGATAAGCAGCTTTCTTTTGGCTATAGTGTCGTGGACAGTGATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt                                                                                                           <  5:491791/150‑1 (MQ=255)
                                    aaaCTGACCGATAAGCCGCTGTCTTTTGGGTATAGTGTCGTGGACAGGCATTCATCTATCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                                                         <  2:671154/151‑1 (MQ=255)
                                     aaCTGACCGATAAGCCGCTTTCTTTTGGGTATCTTGTCGTGGAAAGTCATTCATATTTCTGCCACTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc                                                                                                        <  8:305902/151‑1 (MQ=11)
                                       attaCCGATAATCTGCTTTCTTTTTCCTATAGTGTCGTGGACAGTCATTCTTCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                                      <  1:455879/148‑1 (MQ=11)
                                          aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg                                                                                                    <  2:493180/150‑1 (MQ=35)
                                            cgagcatcagCTTTCTTTATCGTATCGTCTCGCGGAAAGGGAGTCATCTTTCTGCCCCTCCAAAATGACAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                                                 <  6:162455/142‑1 (MQ=11)
                                              atattctgttttctttTGTGTTTAGTGTCGGCGAAAGTGTTACTTCTATCTGGCACTCCAAAAGTAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGt                                                                                               <  2:1191709/142‑1 (MQ=11)
                                               gattcaggtttctttATGTTATAGTCTCGTGGACTGGGTGTAGTATTTCTGCCACTCAAAAAGCATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                                              <  8:1041063/143‑1 (MQ=2)
                                                   ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTGATCTTTCTGCCCATACAAAACAAAAAACCCGCCGAATACGGGTTTTAAGTAAAACAGGTGAAAATTAAAGATAAGACGGGTTATTTTTGGTAGAatgttttttcaacgccagc                                                                                          >  2:184505/1‑136 (MQ=2)
                                                     tctgtctttTGGGTATAGTCTCGAGGACAGTAATTCATATTTCTGCCACTCCAAAAGTACAAACCCGCCGAAGCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGGGGACAGtcattca                                                                                         <  8:773757/146‑1 (MQ=11)
                                                     gCTTTGTTTTGGGTATTGTGTCTTGTACAGTCTCTCATCTTTCTGGCACTCCAAAATTAAAAACCCGCCGTACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcat                                                                                        <  6:740644/151‑1 (MQ=11)
                                                         tctttTGACTAGAGTCTCGGCGACAGTCACTCCTCTATCTGGCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCtt                                                                                     <  6:611090/150‑1 (MQ=11)
                                                          catttttttatcttgtCGTGGACTGTGTTTGATCTTTCTGCCACTCAAAAATTAAAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                                                    <  3:556635/137‑1 (MQ=11)
                                                                       tgacggcgACAGTCAGTCATATATCTCGCACTCAAAAAGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                                                        <  4:1179581/142‑1 (MQ=11)
                                                                         tagtctaatttcattcatGTATCTGTACCTATATAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                                     <  4:140792/140‑1 (MQ=11)
                                                                              taatgtgttttttatttCTGGCACTCAAAAATTTTAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaa                                                                <  4:761046/136‑1 (MQ=1)
                                                                                aaggctttgatCTTTCTGCCCCTCAAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                             <  2:108214/142‑1 (MQ=255)
                                                                                  gtgattgttCTATCTGCCACTCAAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                                           <  4:494442/143‑1 (MQ=11)
                                                                                  gagattcatATTTCTGGCACTCAAAAATAACAAACCCGCCGAATCGGTTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                                           <  6:698701/148‑1 (MQ=11)
                                                                                     attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTTATATTTTTTTTACGTAAAACAGGTGAAACTGACTTATAAGTAGCTTTTTTTTGTGTATATTGTCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                                              <  3:248678/145‑1 (MQ=11)
                                                                                      gtctaatttatcccacTCCAAAATTAAAAAGCGGCCGTAGCGGTTTTTTACGTAAAACAGGTCAAAATGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                                        <  2:896072/135‑1 (MQ=11)
                                                                                      gtcctatatctcCCCCTCCACAAGCACAAACCCGCCGAACCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                                        <  8:454163/142‑1 (MQ=11)
                                                                                       tggtCTTTCTCGCCCTCCAAAAGGACAAACCCGCCGTACCGGGTTTTTACGTAAAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                       <  8:344498/147‑1 (MQ=11)
                                                                                       atttattaatggcacTCAAAAATTATAAACCCGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                       <  4:1257326/136‑1 (MQ=11)
                                                                                        aaactaacTGCCACTCAAAAAGTACAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                     <  2:759433/144‑1 (MQ=11)
                                                                                          acTTACTGCCACTCAAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑Cg                                                  <  4:492150/150‑1 (MQ=11)
                                                                                           atttCTCCCACTCAAAAATTTTAAACCCGCCTTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTTGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑Cg                                                  <  4:819707/149‑1 (MQ=11)
                                                                                           atatctcCCACTCAACAAGAAAAAACCCGCCGAACCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑Cgg                                                 <  6:83696/148‑1 (MQ=11)
                                                                                                  ccactcaacaattataaaCCCGCCGAAGCGGTTTATTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTa                                          <  6:74476/136‑1 (MQ=11)
                                                                                                      tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAGAGAGACCGATAAGCCGATTTCTTTTGGGTATAGTGTCGTGGACTGGGATTCCGCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCGG‑GGGGGGTTTACGtg                                      >  2:434187/1‑150 (MQ=11)
                                                                                                       ccAAAAGAACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGGCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGGCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTaa                                     <  6:286775/151‑1 (MQ=255)
                                                                                                        cAAAAGGACAAACCCGCCGCACCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTaaa                                    <  5:1325595/151‑1 (MQ=255)
                                                                                                         aaaattttaaaCCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAAc                                   <  4:819462/143‑1 (MQ=11)
                                                                                                           aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGCCAGTCATTGATCTTTGTGCGCCTGCAAAACATAAAACAAGACGAAG‑CGGCttgtatccatagacc                                 >  2:140481/1‑138 (MQ=11)
                                                                                                            agatcaaaCCCGCCGTAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAg                                <  7:1250212/146‑1 (MQ=255)
                                                                                                             gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTAATCTTTATTCCCCTCCAAAAACAAAAAACAACACCAG‑CGGGGTTTTAAATAAAACAgg                               >  4:674698/1‑151 (MQ=11)
                                                                                                                taaaCCCGCCGTATCTGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTg                             <  7:700910/149‑1 (MQ=11)
                                                                                                                gaaaCCCGCCGCACCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGa                            <  8:304612/150‑1 (MQ=255)
                                                                                                                aaaaCCCGCCGCATCTTTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGa                            <  3:304974/151‑1 (MQ=11)
                                                                                                                 aatccgcgcgtaccGGGTTATTACGTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaa                           <  8:1405657/138‑1 (MQ=11)
                                                                                                                 aaaCCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGaa                           <  5:199821/151‑1 (MQ=25)
                                                                                                                    ggcgccgATGCGATCTTTTACGTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACt                        <  2:1414598/149‑1 (MQ=11)
                                                                                                                    cccGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGAGCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACt                        <  3:1188847/151‑1 (MQ=255)
                                                                                                                     ccgccgTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATNGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTg                       <  1:318812/151‑1 (MQ=255)
                                                                                                                     ccgccgAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTg                       <  6:1086554/151‑1 (MQ=25)
                                                                                                                       gccgTAGCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAc                     <  7:903371/151‑1 (MQ=255)
                                                                                                                       gccgTAGCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAc                     <  3:844912/151‑1 (MQ=255)
                                                                                                                        ccgcaccGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                    <  6:273233/145‑1 (MQ=255)
                                                                                                                        ccgAAGCGGGTTTTTACGTAAAAGAGGTGAAACTGCCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                    <  6:569942/151‑1 (MQ=25)
                                                                                                                          taagatgttttttACGTAAATCAGGTTAAACTGACCGATAAGCTGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGAcc                    <  4:938162/141‑1 (MQ=11)
                                                                                                                             acGGGTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATaa               <  2:590359/150‑1 (MQ=25)
                                                                                                                              cGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGATTTCATTTGGGTATCGTGTCGCCGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAg              <  5:146720/151‑1 (MQ=11)
                                                                                                                                  cttccactcacaaCACGTAAAACTGACCGATAAGCCGCTTTCTTTTTTTTTTATTTTATTGGACATTCATTGATATTTCTGCAATTCCAAAATCAATTACCCTACTAAT‑CGGTTTTTTTATTTAATCAGGTTAAACTGACCGATAAGcc            <  4:1031102/138‑1 (MQ=2)
                                                                                                                                          tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAGAGTCATTCATTTTTTCGCGCCGGCAAACCCAAAAACCCGCAGAAGAGGGGTTTTTACTTAAA‑CACGGGCAACTTGCCCATAAGCCGGGTTGtttt  >  6:1277486/1‑151 (MQ=2)
                                                                                                                                           aaaaCAGCTGAAACTGACATACAAGCCCCTATCTCTTGTCTTTAGTGGAGACGACTGGGATTCATCTATATCCCACTCAAAAATCACAAACCCGCCGCAG‑CGGGTTTTTACGTAAACCAGGTCAAACTGACCGATAAGCCGCTTTCtttt  <  2:279745/150‑1 (MQ=9)
                                                                                                                                                    |                                                                                                                                             
AAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAG‑CGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTT  >  NC_000913/3269846‑3270134

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑