Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC43.5% 111.8 / 12.4 62T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/27);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                           |                                                                                                                                                      
gtttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACATGTCGCCGCGCGTACCGGATACGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       <  1:282733/149‑1 (MQ=255)
 cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                      >  1:552235/1‑151 (MQ=255)
  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  2:816017/1‑151 (MQ=255)
   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  1:453431/1‑151 (MQ=255)
   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  1:446461/1‑151 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                  <  1:929970/151‑1 (MQ=255)
      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                 <  2:382406/151‑1 (MQ=255)
               ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                        <  1:323576/151‑1 (MQ=255)
                 tCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                      >  2:957034/1‑151 (MQ=255)
                     cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                   >  2:1109170/1‑150 (MQ=255)
                         tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                              >  1:3333/1‑151 (MQ=255)
                            ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           >  1:810612/1‑151 (MQ=255)
                                   ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    >  2:1048165/1‑151 (MQ=255)
                                        ccGCAAAAGTGGTTAGAGTTCCCCGGCTCCTACGGTAAACGGCCGCCGCGCGTCCGGAACGCTTGTTGTATCCCAACCCTTAAACCGTCTACGAGGCAAATATGCAATCCCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 <  1:581595/150‑1 (MQ=255)
                                        ccGAAAAATGGTGAGTGTTCCCCGGCGCGTACGGGAACAGTTCCCCGCGCGTACCGAACCCGTGGTGTATCTCAACCCTTTAAACCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                <  2:883595/150‑1 (MQ=255)
                                         cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               >  2:293826/1‑150 (MQ=255)
                                            aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           >  1:768043/1‑151 (MQ=255)
                                            aaaaGGGTGAGTGTTCCCCTGCGCGTACGGAAACGGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCTTTAACCCGTTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                             <  2:494454/149‑1 (MQ=255)
                                             aaaGGGGGAGTGATCCCCGGCGCGTACTGAACCGGGTCCCCGCGCGTACCGGAACGGGGGTGTATCTCGACCCGTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  2:1067830/151‑1 (MQ=255)
                                                       gtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  2:1012051/1‑147 (MQ=255)
                                                           tccccTGCGCGTACTGGAACGGGTCGCCGCGCGTCCGGGACCCGGGGTGTATCTCGACCCTTAAACCCGCTACAAGGATCAAATAATAAACCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  2:840604/148‑1 (MQ=255)
                                                             cccgccgcgTACGGGACCGGGCCCCCGCGCGTACCGGACCGGGGGTTTACCTCGAACCGTTAAACCGTTATCAGGACAAAATGATAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                             <  1:145853/144‑1 (MQ=255)
                                                              ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          >  2:354114/1‑151 (MQ=255)
                                                                 gcgcGTACTAAAACGGGCCCCCGCGCGCCCCGAAACGGGGGTTTATCCCAACCTGTTAACCCTCTAACGGGATAAAATACACAACCCCCTCGGCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  1:116017/144‑1 (MQ=255)
                                                                     gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   >  1:249019/1‑151 (MQ=255)
                                                                      tacggaaagaGGCCGCCGCGCGTCCCGAACCGGGGTTTTATCTAGAACCTTTAACCCCTTACAGGGATAAAATAATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                      <  1:15691/138‑1 (MQ=255)
                                                                          ggacaaGGTCCCGGCGCGTACGGAACCCGGGTTTTATTACGAACCTTAAAACCCTTATCAGAATAAAATATTCACCCCCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                        <  2:106141/136‑1 (MQ=255)
                                                                            accgggccgccgcGCGTCCGGAACCGGGGTTTTATCTCGACCCTTAAACCCGCTACAAGAATAAAATGATAAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:336094/144‑1 (MQ=255)
                                                                             acgggTCGCCGCGCGACCGGAACCCGGGGTGTACCCCAACCCTTAAACCCGTTATCAGAATAATATAAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                               <  1:259881/144‑1 (MQ=255)
                                                                                ggTCCCCGCGGGTACGGGAACGGGGGTGTATCTCGACCCGTTAACCCGCTACAAGAATCATATGAACAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                             <  1:35613/146‑1 (MQ=255)
                                                                                    gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  1:200251/1‑151 (MQ=255)
                                                                                    gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  2:1040992/1‑151 (MQ=255)
                                                                                     ccgcGGGTCCCGAACCGGTGGTTTTTCCCAACCTGTAAAACCGCTATAGGGATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   <  1:406997/151‑1 (MQ=255)
                                                                                     ccgcGCGTCCCGGACCCGTGGTTTATCTCGACCCTTAAACCCGCTATCAGGATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                     <  2:90448/149‑1 (MQ=255)
                                                                                         gcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:39326/1‑150 (MQ=255)
                                                                                          cgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:219000/1‑149 (MQ=255)
                                                                                          cgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  2:565268/1‑149 (MQ=255)
                                                                                             aCCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            >  2:418623/1‑150 (MQ=255)
                                                                                                 gAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       >  1:1019727/1‑151 (MQ=255)
                                                                                                   acccGTGGTGTCTCTCGACCCGTTAACCCGCTATCAGAATCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:1073713/144‑1 (MQ=255)
                                                                                                        gtggttTATCTCGACCCTTTAACCCGCAAAAAGGATAATATTATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     <  2:108698/148‑1 (MQ=255)
                                                                                                        gtggtgTATCCCGAACCGTAAACCCGCTATCAGGATCATATAACAAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   <  2:1019809/150‑1 (MQ=255)
                                                                                                          ggtgTATCTCGACCCGTAAACCCGCTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:458652/151‑1 (MQ=255)
                                                                                                            ggtttacccagaCCCTTAAACCCGCTATCAGGATAATTTGAACACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    <  2:977726/136‑1 (MQ=255)
                                                                                                               tATCTCAACCCTTAAACCCGCTATAAGAATCATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  2:21516/149‑1 (MQ=255)
                                                                                                                  cccaaaCCGTTAACCCGCTATCAGGATAATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          <  2:623568/146‑1 (MQ=255)
                                                                                                                      aaCCGTAAACCCGCTATCAGGATCATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                       <  2:442079/149‑1 (MQ=255)
                                                                                                                        ccGTAAACCCGTTATCAGGATAATATGACCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   <  2:15110/151‑1 (MQ=255)
                                                                                                                            tAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAActg                               >  1:669903/1‑149 (MQ=255)
                                                                                                                                 ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  1:1050509/1‑151 (MQ=255)
                                                                                                                                 ccGCTATCAGAATAAAAGAATAACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:560523/149‑1 (MQ=255)
                                                                                                                                   gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         >  2:410898/1‑151 (MQ=255)
                                                                                                                                     cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        >  2:1089737/1‑151 (MQ=255)
                                                                                                                                       aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      >  2:1044512/1‑151 (MQ=255)
                                                                                                                                       aTCAGGATAATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        <  2:329945/149‑1 (MQ=255)
                                                                                                                                           ggacaatataaaaaaCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                  <  1:1090332/138‑1 (MQ=255)
                                                                                                                                           ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                  >  1:827976/1‑151 (MQ=255)
                                                                                                                                              tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt               >  1:1132444/1‑151 (MQ=255)
                                                                                                                                                 tataaccAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 <  2:453459/140‑1 (MQ=255)
                                                                                                                                                  atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:719503/1‑151 (MQ=255)
                                                                                                                                                  atGACCACCCCCCCCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  2:1018537/142‑1 (MQ=255)
                                                                                                                                                   tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:1062395/1‑150 (MQ=255)
                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  1:844966/1‑151 (MQ=255)
                                                                                                                                                         acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc    >  2:142369/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  2:65516/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655

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