Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT56.9% 20.0 / 126.9 116S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (10/40);  new base T (66/0);  total (76/40)
Fisher's exact test for biased strand distribution p-value = 4.66e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-15

Only 100 of 125 total aligned reads displayed.

CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604196‑4604492
                                                                                                                                                    |                                                                                                                                                      
cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttttt                                                                                                                                                      >  1:1586791/1‑151 (MQ=255)
cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcgt                                                                                                                                                      <  1:616069/151‑3 (MQ=255)
  ccaacaaccctGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttttttt                                                                                                                                                    >  2:528631/11‑151 (MQ=255)
  cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                                    >  1:927674/1‑151 (MQ=255)
  cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                                    >  2:1478518/1‑151 (MQ=255)
              cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                        <  2:1172389/151‑1 (MQ=255)
               gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCTGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCGCtttttctatttttcccctctcaa                                                                                                                                       >  2:769741/1‑139 (MQ=255)
               gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGTGaa                                                                                                                                       >  1:1267187/1‑151 (MQ=255)
               gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGTCATGAGaa                                                                                                                                       >  2:491131/1‑151 (MQ=255)
                tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                      <  1:1407240/151‑1 (MQ=255)
                tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                      <  1:1220274/151‑1 (MQ=255)
                tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                      <  1:1057433/151‑1 (MQ=255)
                  aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                    <  2:591119/151‑1 (MQ=255)
                  aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                    <  1:1461358/151‑1 (MQ=255)
                    gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAg                                                                                                                                  >  1:290199/1‑151 (MQ=255)
                    gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAg                                                                                                                                  >  1:801561/1‑151 (MQ=255)
                    gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTCTTTTTTTCCATggatagtggg                                                                                                                                  >  2:869950/1‑142 (MQ=255)
                      gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAgcg                                                                                                                                >  2:559971/1‑151 (MQ=255)
                         gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGAGAAGGGAgcgccg                                                                                                                             >  1:592464/1‑149 (MQ=255)
                          taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                             <  2:1488623/150‑1 (MQ=255)
                            tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                                                                          <  2:1059727/151‑1 (MQ=255)
                              ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCTTGAGAAGTGAGCGCAGAGAgg                                                                                                                        >  2:912937/1‑151 (MQ=255)
                               gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGgggg                                                                                                                       >  1:30005/1‑151 (MQ=255)
                               gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAgggggagccggggggg                                                                                                                       >  1:1089806/1‑136 (MQ=255)
                                tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                       <  2:824776/150‑1 (MQ=255)
                                 aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                     <  2:1127943/151‑1 (MQ=255)
                                 aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                     <  2:20755/151‑1 (MQ=255)
                                 aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                     <  1:1603316/151‑1 (MQ=255)
                                   tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                   <  1:441046/151‑1 (MQ=255)
                                       agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCCTGAGAAGTGGGGGCAGAGGGGGGTTATCtt                                                                                                               >  1:923158/1‑151 (MQ=255)
                                            aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTTTGAAGGGGGCGCCGGGAGGGGttttttttcccg                                                                                                          >  1:575261/1‑141 (MQ=255)
                                             ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATTTTGCctga                                                                                                         >  2:107693/1‑151 (MQ=255)
                                              tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                        <  1:1514572/151‑1 (MQ=255)
                                                 gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTTCCTGAGAAGTGGGGGCCGAGGGGGGTTTTCTTGCctgactga                                                                                                     >  2:654617/1‑151 (MQ=255)
                                                    aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTCGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCCGAGCGGGATTATCTTGCctgactgga                                                                                                    >  2:657276/1‑147 (MQ=255)
                                                     gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATTTTTCCTGACTGAAATg                                                                                                 >  1:1594806/1‑151 (MQ=255)
                                                     gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGTCTGACTGAACtt                                                                                                 >  2:231551/1‑150 (MQ=255)
                                                     gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGGGAGGGGTTTTTTTTCcttactgga                                                                                                    >  2:114168/1‑146 (MQ=255)
                                                      tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                                  <  2:558375/149‑1 (MQ=255)
                                                           gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGGGGGGGGTTTTCTTTTCTGGCTTAAATGGCGGt‑t                                                                                          >  1:295842/1‑151 (MQ=255)
                                                              tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAg                                                                                       <  1:381623/151‑1 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCATGAGCAGTGAGCGCGGAGAGGGGTTATCTTTCCTGACTGTACTGGCGCG‑TACG‑ttt                                                                                   >  2:1137071/1‑151 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATCAGAAGGGGGGGCCGGGAGGGGTTTTCTTTCCTTACTTGACCGGGGGT‑TAAA‑ttt                                                                                   >  1:702594/1‑151 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAAATGGCGAT‑TAAG‑gtt                                                                                   >  1:1152425/1‑151 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTCCCCTGGGAAGTGAGCGCCGAGAGGGATTATCTTGCCTTGATGAACTGGGGGT‑TAAG‑ttt                                                                                   >  1:356909/1‑151 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGGGAAGGGGGGGCCGGGAGGGGTTATTTTGCCCGACTGAACCGGGGGT‑TAAA‑ttt                                                                                   >  1:6020/1‑151 (MQ=255)
                                                                 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCCGAGGGGGGTTATTTTGTCTGACTGAACTGGCGGT‑TAAG‑ttt                                                                                   >  2:657293/1‑151 (MQ=255)
                                                                   cTTATAAGCTTCATACCTATACTAGCCAGGAGTCATTTCAGGCTGCCATGTCGCGGGTCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTCCCCTGGGAAACGAAGCCAGAGCGGGGTTTTCTCGCCTGACCGAACCGGCGGC‑TCAG‑TTTc                                                                                  >  2:488691/1‑150 (MQ=255)
                                                                   cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGT‑TAAG‑tttt                                                                                  >  2:74556/1‑149 (MQ=255)
                                                                   cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGGGCGCCGAGAGGGGTTATCTTTCCTTGCTTCACTTGGGAT‑TAAG‑GTTc                                                                                  >  1:750019/1‑150 (MQ=255)
                                                                   cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTGACTGTACTGGCGGT‑TAAG‑tttt                                                                                  >  2:1570255/1‑149 (MQ=255)
                                                                          gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAGGGGTTATTTTGCCTGACTGAACCGGCGGT‑TTAG‑TTTCCGCGTAc                                                                           >  1:1229907/1‑150 (MQ=255)
                                                                          gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTTCTGACTGAACTGGCGAT‑TTAG‑TTTTCGCGTAc                                                                           >  2:1447393/1‑150 (MQ=255)
                                                                          gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTTTTTTTCTGTCATGTGCAGAGAGCGCAGCGAGGGGAGAACTTTCCTGAATGAAATGACGATGAAAG‑TTTctattaa                                                                            >  2:1073682/1‑144 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGGTGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATTAGAAGTGCGCGCAGAGGGGGGTTTTCTTGCCTGACTGAACTGGCGGT‑TTAG‑GTTCCGCGtgac                                                                          >  2:1440289/1‑147 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGAGAAGAGAGGGCGGAGAGGGATTTTCTTTCCTGGCTTTACTTGTGAT‑TAAG‑ATTTCCCGGAcc                                                                          >  2:1169416/1‑150 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTTTCTGTACAGGCGAT‑TTAG‑TTTTCGCGTAcc                                                                          >  2:119107/1‑150 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                     <  2:794579/109‑1 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTGCCCTGAGAAGGGGGGGCAGAGGGGGGTTATCTTGCCTGACTGAAATGGCGAA‑TTAG‑TTTTCCCGTAcc                                                                          >  1:1066328/1‑150 (MQ=255)
                                                                               aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCAGAGAGGGGTTATTTTTGCCGGCTGTAACGGCGAT‑TAAG‑TTTTCGCGGACCCGGCg                                                                     >  2:1278121/1‑151 (MQ=255)
                                                                               aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCg                                                                     <  1:925974/151‑1 (MQ=255)
                                                                                tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCTGTTTTTTATTTTATCCCATTAGAAGGGAGCGCAGAGAGGGATTATTTTTTCTTACTGTAATGGGCGT‑TAAA‑TTTTCGCCGACCCCGCGt                                                                    >  1:257452/1‑151 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCCTTAGAAGTGAGCGCCGCGAGGGGGTTTCTTGGCTTACTTAACTAGCGAT‑TAAA‑TTTTCGCGGACACGGCggt                                                                   >  2:1085894/1‑149 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGGGAGCGCCGAGAGGGGTTTTCTTGCCTGACTTAAATGGCGAT‑TAAG‑TTTTCCCGGACCCCGCGtt                                                                   >  1:661086/1‑151 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGGA‑TAAG‑GTTTCGCGGAACCGGCGtt                                                                   >  2:1186066/1‑151 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGGGGGGGGTTATCTTTCCTGTCTGAACTGGCGGT‑TAAG‑TTTTCGCGGGCCCGGCggt                                                                   >  2:1143640/1‑149 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTCTGCGGCGGCCCCTTTTTCTTTTTCTTCCACGAGGAAGGAGCGCCGGGAGGGGTTATCTTGCCTCACTGAAATGGCGAC‑TAAG‑TTTCCCCGTACCCCGNGnt                                                                   >  1:333180/1‑151 (MQ=255)
                                                                                   cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAATGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTTCGCGTACCCGGCGgttt                                                                 >  2:157131/1‑150 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTTTCTTTTCCGTCCGTACTGGCGAT‑TAAG‑TTTCcgcgccccccccgctttgt                                                               >  1:1318844/1‑136 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGTGAAGGGAGCGCAGAGAGGGGTTATCTTTCCCGGCTGAACTGGCGGT‑TTAG‑TTTTCGCGGACCCCGCGGTTTGt                                                               >  1:452959/1‑151 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGAT‑TAAG‑GTTCCGCGGACCCGGCGGTTTGt                                                               >  2:174574/1‑151 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGt                                                               <  1:1389218/151‑1 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCCGAGGGGGGTTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑GTTCCGCGTACCCCGCGGTTTGt                                                               >  2:890065/1‑151 (MQ=255)
                                                                                      taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGtt                                                              <  1:658063/151‑1 (MQ=255)
                                                                                          aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCCGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTTCGCGTACCCGGCGCTTTGTttttt                                                          >  2:553541/1‑150 (MQ=255)
                                                                                          aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTTCCAGAGGAGGGGGGGGCGAGAGGGGTTATTTTGCCCTAATTGATTGGGGAT‑TAAGATTTCCCCGAACCC‑GCGTTTTGTTTTTg                                                          >  2:1408253/1‑151 (MQ=255)
                                                                                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTg                                                          <  1:174610/149‑1 (MQ=255)
                                                                                                    cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCCGAGGGGGGTTTTCTTTCCTGACTGAACTGGGGAT‑TAAG‑TTTTCGCGGACCCGGCGGTTTGTTTTTGCggggggtga                                                >  1:1284911/1‑144 (MQ=255)
                                                                                                    cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCCGAGAGGGGTTTTCTTTCCTGACTGAACTGGGGAT‑TAAG‑TTTTCCCGGAACCGGCGGTTTGTTATTTGGgatggtgg                                                >  2:1389689/1‑150 (MQ=255)
                                                                                                    cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                                <  2:389582/151‑1 (MQ=255)
                                                                                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                               <  2:1052434/151‑1 (MQ=255)
                                                                                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                               <  2:454241/151‑1 (MQ=255)
                                                                                                      tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGAGAAGGGGGCGCAGGGGGGGGTTTTCTTTCCTGACTTAACTGGCGGT‑TAAG‑TTTTCCCGTACCCCGCGTTTTTTTTTTTCGGGTGATGata                                              >  1:128630/1‑151 (MQ=255)
                                                                                                      tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATTTTTCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCCGCGTTTAGTTATTGCGGATGATGata                                              >  2:228208/1‑151 (MQ=255)
                                                                                                        tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGTGAGGGCAGGGAGGGATTTTCTTGCCTGACTGTACTGGCGGT‑TTAG‑TTTTCGCCTACCCCGCGCTTCGTTATTTCGGGTGGTGGTAtt                                            >  2:1581498/1‑151 (MQ=255)
                                                                                                           ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGGCCCGAGAAGTGAGCGCAGGGAGGGGTTATCTTTCCTGACTGGACTGGCGGT‑TAAG‑TTTTCGCGGACGCGGCGTTTTGGTATTGCGGATTATGGTAATGaa                                         >  2:355093/1‑151 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTGCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  1:1137600/150‑1 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                       <  2:546717/151‑1 (MQ=255)
                                                                                                                  aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                  <  1:890539/151‑1 (MQ=255)
                                                                                                                    gTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGGGAAGGGGGCGCAGAGAGGGGGTATCTTGCCTGACTAAACTGGCGGT‑TAAG‑GTTTCGCGTACCCCGCGTTTAGTTATTGCGGTGGTTGTTATTGTAGCTCgg                                   >  2:1333378/1‑148 (MQ=255)
                                                                                                                        cgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                             <  1:1143213/150‑1 (MQ=255)
                                                                                                                        cgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTg                            <  1:528826/151‑1 (MQ=255)
                                                                                                                            tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc                        <  2:372934/151‑1 (MQ=255)
                                                                                                                            tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc                        <  2:308610/151‑1 (MQ=255)
                                                                                                                                 tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  1:43006/151‑1 (MQ=255)
                                                                                                                                         ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa           >  1:375603/1‑151 (MQ=255)
                                                                                                                                             aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc       <  1:660588/151‑1 (MQ=255)
                                                                                                                                                  ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTTGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  2:89049/151‑1 (MQ=255)
                                                                                                                                                  ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:842416/151‑1 (MQ=255)
                                                                                                                                                  ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGAGGATATTGAAGCTCGGGTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  2:1584868/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAT‑TAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604196‑4604492

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