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

TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436371‑436666
                                                                                                                                                   |                                                                                                                                                      
tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATg                                                                                                                                                   >  2:717890/1‑151 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTCCC‑TGTATCGGGCAAAt                                                                                                                                                    <  2:502994/149‑1 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCaaa                                                                                                                                                     <  3:1444360/148‑1 (MQ=255)
 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAAt                                                                                                                                                    <  3:856010/149‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccg                                                                                                                                               <  3:1224442/151‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccg                                                                                                                                               >  3:831247/1‑151 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccg                                                                                                                                               <  2:473199/151‑1 (MQ=255)
    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccg                                                                                                                                               >  2:480967/1‑151 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgc                                                                                                                                              <  2:424841/151‑1 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgc                                                                                                                                              >  2:1110611/1‑151 (MQ=255)
     aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGCGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgc                                                                                                                                              >  3:1254050/1‑151 (MQ=255)
       gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGAccgccg                                                                                                                                            <  3:1456943/151‑1 (MQ=255)
           aaCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGatat                                                                                                                                        >  3:1537401/1‑151 (MQ=255)
              cGGAACAGGCGCGGCGCGGGGCGCTCTCGGGCGGTAAAGATTACGAGTTGTGTTTCACTGGGCCGGAACGGACCCGTGGCGCGGTGGATGTGGCTCCGGGCCCCCGGGGCG‑TCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGatat                                                                                                                                        <  3:1015512/148‑1 (MQ=255)
               ggAACGGGCGCGGGGTGGGGCGCTCTCGGGGGGTGAAAATTAGGAGTTTTTTTTCACTGTGCGGAAACTGACCGGTGGCGCGTTGGATGGGGCTCTGGGACCCCTGGGCG‑TCCCGTTTCCC‑TGTTTGGGGAAAATGACCGCCGATATCGaa                                                                                                                                    <  2:549488/151‑1 (MQ=255)
                    aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                               >  3:460837/1‑151 (MQ=255)
                            cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                        >  3:1549586/1‑150 (MQ=255)
                                gggCGCTCTCTGGGGGTGAAGATTACGAGTTGTTTTTCCCTGGGCCGGAACTGACCCGTGGCGCGTTGGTTGTGGCTCTCGGACCCCTGGGCG‑TCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                     <  2:1520721/149‑1 (MQ=255)
                                gggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACAACTGGGCG‑TACCGTTTTAC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                   >  3:768023/1‑151 (MQ=255)
                                gggCCCTCTCGGGGGGTGAAGATTACGGGTTGTGTTTCTCTGTGCCGGAACTGACCCGGGGCGCGCTGGTGGGGCCTCTCGGCCCCGGGGGCG‑TCCCTTTTCCC‑TGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                      <  2:1583330/148‑1 (MQ=255)
                                        tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           >  2:642132/1‑151 (MQ=255)
                                         ctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           >  2:1095188/1‑150 (MQ=255)
                                          tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                          >  3:884263/1‑150 (MQ=255)
                                           gggggTGAAGATTACGGTTTGTGTTTCACTGTGCCGGAACTGACCCGTGGCGCGCGGGTTGTGGCTCTCGGACCCCGGGGCG‑TACCTTTTCCC‑TGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           <  2:1079614/145‑1 (MQ=255)
                                           ggCGGTGAAGATTACGAGTTGTGTTTCTCTGTGCGGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTGGGACCCCTGGGCG‑TCCCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                              <  2:1412501/145‑1 (MQ=255)
                                              ggTGAAGATTACGGGTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCG‑TCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                                       <  2:546879/149‑1 (MQ=255)
                                              ggTGAAGATTACGGGTTGTGTTTCTCTGTGCCGGAACTGACCGGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCG‑TCCCTTTTCCC‑TGTATCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                      <  2:379750/150‑1 (MQ=255)
                                               gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGt                                                                                                    >  3:1565731/1‑151 (MQ=255)
                                                tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   >  2:77252/1‑151 (MQ=255)
                                                    gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               >  3:498555/1‑151 (MQ=255)
                                                      ttacgttttttttttCCCTGGGCGGAAACTGACCGGGGGCGCGTTGGATGGGGCTCTCGGCCCCCGGGGCG‑TCCTTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  3:480932/135‑1 (MQ=255)
                                                      ttACGGGTTGTTTTTCCCTGTGCGGGAACTGACCGGGGGCGCGTTGGTTGGGGCTCTCGGACCCCGGGGCG‑TCCCTTTTCCC‑TGTACGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               <  3:370582/149‑1 (MQ=255)
                                                      ttACGAGTTGTTTTTCCCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCG‑TCCCTTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                               <  2:346036/149‑1 (MQ=255)
                                                      ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCGGGTTGTGGCTCTCGGACCCGTGGGCG‑TCCCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              <  2:902351/150‑1 (MQ=255)
                                                          gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                         >  3:1122505/1‑151 (MQ=255)
                                                            gTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                       >  2:189617/1‑151 (MQ=255)
                                                             ttGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGTTGGATGGGGCTCTGGGACCCGTGGGCG‑TCCCGTTTCCC‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  2:799027/149‑1 (MQ=255)
                                                             tggtgtTTCACTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTGGGACGCTTGGGCG‑TCCCGTTTCCT‑TGTATCGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  2:554077/147‑1 (MQ=255)
                                                                  tttCGCTGTGCCGGAACTGACCCGTGGCGCGTTGGTTGGGGCTCTCGGACGCCGGGGCG‑TCCCGTTTCCG‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                                  <  3:996910/150‑1 (MQ=255)
                                                                    tcccTGTGCCGGAAGGTACGCGGGGCGCGTTGGATGTGGCTCCGGGGCATCGGGGGG‑TCCCTTTTACC‑TTTTGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                                <  2:1474663/147‑1 (MQ=255)
                                                                         gggcgggAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCTTGGGCG‑TCCCTTTTACC‑TGTTTCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa                                                                            <  2:203826/144‑1 (MQ=255)
                                                                             cGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                      >  2:1237969/1‑151 (MQ=255)
                                                                              ggAACTGAACCGTGGCGCGTTGGATGTGGCTCGCGGACCCGGGGGCG‑TACCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                        <  2:1295073/148‑1 (MQ=255)
                                                                                aaCTGACCGGGGGCGCGGTGGTTGTGGCTCTCGGACCCGGGGGCT‑TACCTTTTCCT‑TGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   <  3:374706/151‑1 (MQ=255)
                                                                                aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   >  3:1512825/1‑151 (MQ=255)
                                                                                  cTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                                  >  2:1313268/1‑150 (MQ=255)
                                                                                   tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                                >  2:1621653/1‑151 (MQ=255)
                                                                                   tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                                >  2:1072248/1‑151 (MQ=255)
                                                                                     aaCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                              >  2:168148/1‑151 (MQ=255)
                                                                                       cggTGGCGCGCTGGTTGTGGCTCTCGGCCCCTTGGGCG‑TACCGTTTCCC‑TGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                               <  2:1304083/146‑1 (MQ=255)
                                                                                       ccGGGGCGCGCGGGTTGTGGCTCGGGAACACGGGGGCT‑TCCTGTTTGCT‑GGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   <  3:387871/144‑1 (MQ=255)
                                                                                          tGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCTGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                          >  2:650222/1‑150 (MQ=255)
                                                                                          gggCGCGGTGGTTGGGTCTCCCGGAC‑CCGGGGGG‑TTCCTTTTTCCTTGTATGGGAAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                          <  3:457042/149‑1 (MQ=255)
                                                                                                gtgggtgggggCTCTCGAAC‑CCTGGGGGTTACTTTTTCCT‑TGTATGGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGa                                                      <  2:652137/139‑1 (MQ=255)
                                                                                                  tGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                                  >  3:1076530/1‑150 (MQ=255)
                                                                                                   ggttggggCTCGCGGCCCCGTGGGCG‑TCCTTTTTCCC‑TATTTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                                   <  2:1319787/142‑1 (MQ=255)
                                                                                                   ggATGTGGCTCTCGGCCCCCTGGGCT‑TACCTTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                                 <  2:830080/150‑1 (MQ=255)
                                                                                                      tgtgGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGa                                              >  3:1063384/1‑150 (MQ=255)
                                                                                                         ggCTCTCGGACCCTGGGGCG‑TACTTTTTCCT‑TTTTTGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                             <  2:1461866/148‑1 (MQ=255)
                                                                                                         ggCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                            >  2:626205/1‑149 (MQ=255)
                                                                                                            tctcGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCATGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       >  2:949162/1‑151 (MQ=255)
                                                                                                             cgggcccccgtGGGCG‑TCCGGTTTCCT‑TTTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                          <  2:275531/137‑1 (MQ=255)
                                                                                                             cgcGGACACGGGGGCG‑TCCCTTTTCCT‑TGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                       <  2:1453104/148‑1 (MQ=255)
                                                                                                                ggcccccggggCG‑TCCCTTTTCCT‑TGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt                                                     <  2:895007/125‑1 (MQ=255)
                                                                                                                ggACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                                   >  2:94613/1‑151 (MQ=255)
                                                                                                                  ggcccccggggcttccGTTTTCCTTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                            <  2:1589510/128‑1 (MQ=255)
                                                                                                                   cacCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  3:140935/1‑149 (MQ=255)
                                                                                                                     ccTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                                  >  2:861147/1‑147 (MQ=255)
                                                                                                                       ggggCG‑TCCCTTTTCCT‑TTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                             <  2:192942/149‑1 (MQ=255)
                                                                                                                        gggCG‑T‑CCGTTTTCCTTGTAGGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa                            <  2:377937/150‑1 (MQ=255)
                                                                                                                               aCCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                     >  2:1643569/1‑151 (MQ=255)
                                                                                                                                   tttCCC‑TGTATCGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg                  <  2:360854/150‑1 (MQ=255)
                                                                                                                                   tttACT‑TGTAGGGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                    <  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)
                                                                                                                                                   |                                                                                                                                                      
TCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCG‑TACCGTTTACC‑TGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  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: ‑