Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA42.9% 102.9 / 37.3 77Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (35/9);  new base A (0/33);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 7.81e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436371‑436666
                                                                                                                                                   |                                                                                                                                                      
tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                                   >  2:717890/1‑151 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTCCCTGTATCGGGCAAAt                                                                                                                                                    <  2:502994/149‑1 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaaa                                                                                                                                                     <  3:1444360/148‑1 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                                    <  3:856010/149‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                               <  3:1224442/151‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                               >  3:831247/1‑151 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                               <  2:473199/151‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                               >  2:480967/1‑151 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                              <  2:424841/151‑1 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                              >  2:1110611/1‑151 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGCGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                              >  3:1254050/1‑151 (MQ=255)
       gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                            <  3:1456943/151‑1 (MQ=255)
           aaCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGatat                                                                                                                                        >  3:1537401/1‑151 (MQ=255)
              cGGAACAGGCGCGGCGCGGGGCGCTCTCGGGCGGTAAAGATTACGAGTTGTGTTTCACTGGGCCGGAACGGACCCGTGGCGCGGTGGATGTGGCTCCGGGCCCCCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGatat                                                                                                                                        <  3:1015512/148‑1 (MQ=255)
               ggAACGGGCGCGGGGTGGGGCGCTCTCGGGGGGTGAAAATTAGGAGTTTTTTTTCACTGTGCGGAAACTGACCGGTGGCGCGTTGGATGGGGCTCTGGGACCCCTGGGCGTCCCGTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGaa                                                                                                                                    <  2:549488/151‑1 (MQ=255)
                    aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                               >  3:460837/1‑151 (MQ=255)
                            cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                        >  3:1549586/1‑150 (MQ=255)
                                gggCGCTCTCTGGGGGTGAAGATTACGAGTTGTTTTTCCCTGGGCCGGAACTGACCCGTGGCGCGTTGGTTGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                     <  2:1520721/149‑1 (MQ=255)
                                gggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACAACTGGGCGTACCGTTTTACTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                   >  3:768023/1‑151 (MQ=255)
                                gggCCCTCTCGGGGGGTGAAGATTACGGGTTGTGTTTCTCTGTGCCGGAACTGACCCGGGGCGCGCTGGTGGGGCCTCTCGGCCCCGGGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  2:1583330/148‑1 (MQ=255)
                                        tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           >  2:642132/1‑151 (MQ=255)
                                         ctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           >  2:1095188/1‑150 (MQ=255)
                                          tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                          >  3:884263/1‑150 (MQ=255)
                                           gggggTGAAGATTACGGTTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCGGGTTGTGGCTCTCGGACCCCGGGGCGTACCTTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           <  2:1079614/145‑1 (MQ=255)
                                           ggCGGTGAAGATTACGAGTTGTGTTTCTCTGTGCGGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTGGGACCCCTGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                              <  2:1412501/145‑1 (MQ=255)
                                              ggTGAAGATTACGGGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                       <  2:546879/149‑1 (MQ=255)
                                              ggTGAAGATTACGGGTTGTGTTTCTCTGTGCCGGAACTGACCGGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                      <  2:379750/150‑1 (MQ=255)
                                               gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGt                                                                                                    >  3:1565731/1‑151 (MQ=255)
                                                tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   >  2:77252/1‑151 (MQ=255)
                                                    gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               >  3:498555/1‑151 (MQ=255)
                                                      ttacgttttttttttCCCTGGGCGGAAACTGACCGGGGGCGCGTTGGATGGGGCTCTCGGCCCCCGGGGCGTCCTTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  3:480932/135‑1 (MQ=255)
                                                      ttACGGGTTGTTTTTCCCTGTGCGGGAACTGACCGGGGGCGCGTTGGTTGGGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTACGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               <  3:370582/149‑1 (MQ=255)
                                                      ttACGAGTTGTTTTTCCCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               <  2:346036/149‑1 (MQ=255)
                                                      ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCGGGTTGTGGCTCTCGGACCCGTGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              <  2:902351/150‑1 (MQ=255)
                                                          gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                         >  3:1122505/1‑151 (MQ=255)
                                                            gTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                       >  2:189617/1‑151 (MQ=255)
                                                             ttGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGTTGGATGGGGCTCTGGGACCCGTGGGCGTCCCGTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  2:799027/149‑1 (MQ=255)
                                                             tggtgtTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTGGGACGCTTGGGCGTCCCGTTTCCTTGTATCGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  2:554077/147‑1 (MQ=255)
                                                                  tttCGCTGTGCCGGAACTGACCCGTGGCGCGTTGGTTGGGGCTCTCGGACGCCGGGGCGTCCCGTTTCCGTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  <  3:996910/150‑1 (MQ=255)
                                                                    tcccTGTGCCGGAAGGTACGCGGGGCGCGTTGGATGTGGCTCCGGGGCATCGGGGGGTCCCTTTTACCTTTTGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                                <  2:1474663/147‑1 (MQ=255)
                                                                         gggcgggAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCTTGGGCGTCCCTTTTACCTGTTTCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa                                                                            <  2:203826/144‑1 (MQ=255)
                                                                             cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      >  2:1237969/1‑151 (MQ=255)
                                                                              ggAACTGAACCGTGGCGCGTTGGATGTGGCTCGCGGACCCGGGGGCGTACCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                        <  2:1295073/148‑1 (MQ=255)
                                                                                aaCTGACCGGGGGCGCGGTGGTTGTGGCTCTCGGACCCGGGGGCTTACCTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   <  3:374706/151‑1 (MQ=255)
                                                                                aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   >  3:1512825/1‑151 (MQ=255)
                                                                                  cTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                                  >  2:1313268/1‑150 (MQ=255)
                                                                                   tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                                >  2:1621653/1‑151 (MQ=255)
                                                                                   tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                                >  2:1072248/1‑151 (MQ=255)
                                                                                     aaCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                              >  2:168148/1‑151 (MQ=255)
                                                                                       cggTGGCGCGCTGGTTGTGGCTCTCGGCCCCTTGGGCGTACCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                               <  2:1304083/146‑1 (MQ=255)
                                                                                       ccGGGGCGCGCGGGTTGTGGCTCGGGAACACGGGGGCTTCCTGTTTGCTGGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   <  3:387871/144‑1 (MQ=255)
                                                                                          tGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCTGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                          >  2:650222/1‑150 (MQ=255)
                                                                                          gggCGCGGTGGTTGGGTCTCCCGGACCCGGGGGGTTCCTTTTTCCTTGTATGGGAAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                          <  3:457042/149‑1 (MQ=255)
                                                                                                gtgggtgggggCTCTCGAACCCTGGGGGTTACTTTTTCCTTGTATGGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGa                                                      <  2:652137/139‑1 (MQ=255)
                                                                                                  tGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  >  3:1076530/1‑150 (MQ=255)
                                                                                                   ggttggggCTCGCGGCCCCGTGGGCGTCCTTTTTCCCTATTTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                                   <  2:1319787/142‑1 (MQ=255)
                                                                                                   ggATGTGGCTCTCGGCCCCCTGGGCTTACCTTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                 <  2:830080/150‑1 (MQ=255)
                                                                                                      tgtgGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGa                                              >  3:1063384/1‑150 (MQ=255)
                                                                                                         ggCTCTCGGACCCTGGGGCGTACTTTTTCCTTTTTTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                             <  2:1461866/148‑1 (MQ=255)
                                                                                                         ggCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                            >  2:626205/1‑149 (MQ=255)
                                                                                                            tctcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCATGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       >  2:949162/1‑151 (MQ=255)
                                                                                                             cgggcccccgtGGGCGTCCGGTTTCCTTTTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                          <  2:275531/137‑1 (MQ=255)
                                                                                                             cgcGGACACGGGGGCGTCCCTTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  2:1453104/148‑1 (MQ=255)
                                                                                                                ggcccccggggCGTCCCTTTTCCTTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                     <  2:895007/125‑1 (MQ=255)
                                                                                                                ggACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                                   >  2:94613/1‑151 (MQ=255)
                                                                                                                  ggcccccggggcttccGTTTTCCTTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                            <  2:1589510/128‑1 (MQ=255)
                                                                                                                   cacCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  3:140935/1‑149 (MQ=255)
                                                                                                                     ccTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  2:861147/1‑147 (MQ=255)
                                                                                                                       ggggCGTCCCTTTTCCTTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             <  2:192942/149‑1 (MQ=255)
                                                                                                                        gggCGTCCGTTTTCCTTGTAGGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa                            <  2:377937/150‑1 (MQ=255)
                                                                                                                               aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  2:1643569/1‑151 (MQ=255)
                                                                                                                                   tttCCCTGTATCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  <  2:360854/150‑1 (MQ=255)
                                                                                                                                   tttACTTGTAGGGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                    <  3:1043164/148‑1 (MQ=255)
                                                                                                                                          tGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt            >  3:157927/1‑150 (MQ=255)
                                                                                                                                               cGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc      >  3:723656/1‑151 (MQ=255)
                                                                                                                                                gggAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg       <  2:1426017/149‑1 (MQ=255)
                                                                                                                                                 ggCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTg    >  3:1593165/1‑151 (MQ=255)
                                                                                                                                                   cAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  2:1535313/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436371‑436666

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