Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT48.9% 39.1 / 63.4 88S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (13/32);  new base T (43/0);  total (56/32)
Fisher's exact test for biased strand distribution p-value = 7.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604198‑4604492
                                                                                                                                                  |                                                                                                                                                           
cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                                         >  1:974762/1‑151 (MQ=255)
cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                                         >  2:830097/1‑151 (MQ=255)
         aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCCGCGGTCTGTTTTTTTTTTTTTACCaag                                                                                                                                                >  2:838107/1‑149 (MQ=255)
            cGTCATGATTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                             <  2:60715/151‑1 (MQ=255)
            cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttttttccctggga                                                                                                                                             >  2:67501/1‑140 (MQ=255)
             gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGaa                                                                                                                                            >  1:925501/1‑151 (MQ=255)
              tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                           <  1:975564/151‑1 (MQ=255)
                aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                         <  1:510973/151‑1 (MQ=255)
                  gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTCTGCCCTGCGAAGTgcg                                                                                                                                       >  1:53269/1‑149 (MQ=255)
                  gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAg                                                                                                                                       >  1:676935/1‑151 (MQ=255)
                  gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGAg                                                                                                                                       >  1:172489/1‑151 (MQ=255)
                    gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCATGAGAAGTGAgcg                                                                                                                                     >  1:308448/1‑151 (MQ=255)
                    gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAgcg                                                                                                                                     >  2:779766/1‑151 (MQ=255)
                       gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGAAGTGGGCGCag                                                                                                                                  >  1:42002/1‑151 (MQ=255)
                             gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCTTTTTTTTTTTTTTGCCATGGGAAATGGGCGGCGGGAggg                                                                                                                            >  1:605870/1‑151 (MQ=255)
                             gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                             <  2:580971/150‑1 (MQ=255)
                             gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGGGAAGGGAgggccggggggg                                                                                                                            >  2:762171/1‑143 (MQ=255)
                               aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGTGAggggt                                                                                                                          >  1:277902/1‑149 (MQ=255)
                                 tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                        <  2:67014/151‑1 (MQ=255)
                                 tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                        <  2:925828/151‑1 (MQ=255)
                                 tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                        <  2:159914/151‑1 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGGGAAGGGAGGGCAGAGAGGGGtttttt                                                                                                                     >  2:382220/1‑147 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGtttt                                                                                                                       >  2:310190/1‑147 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGGGCAGAGGGGGATTA‑t                                                                                                                      >  2:718924/1‑149 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCGTTTTTTTTTTTTTTGCCCCGCTGAAAGTGCGCAGAGAGGGGTTA‑t                                                                                                                      >  1:842715/1‑149 (MQ=255)
                                     agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑Tc                                                                                                                     <  2:903883/149‑1 (MQ=255)
                                      gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGTCATGAGGAGTGAGGGCAGAGAGGGATTA‑TTTTg                                                                                                                  >  1:245213/1‑151 (MQ=255)
                                         tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAGGGAGGGGAGGGGGGGAGTT‑TCtttc                                                                                                                 >  2:567654/1‑147 (MQ=255)
                                                 aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AAc                                                                                                      <  1:718664/151‑1 (MQ=255)
                                                        cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGGGGGGGGTTA‑TCTTGCCTGACTG‑AACT‑GGC‑ggt                                                                                             >  2:1002121/1‑149 (MQ=255)
                                                         gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCCTGAGAAGTGAGAGCAGAGAGGGGTTT‑TCTTTCCTGTCTG‑AACT‑GGC‑gggt                                                                                            >  1:21135/1‑148 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCCGAGAGGGGTTT‑TTTTTCCTTACTG‑AAAT‑GGC‑GG‑TTAAGgtt                                                                                     >  2:660347/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGGGAAGTGTGCGCTGAGAGGGGTTT‑TCTTGCCTTACTG‑AACT‑TGC‑GG‑TTAAGttt                                                                                     >  2:381787/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGGGGAGGGAGCGCCGGGGGGGGTTT‑TTTTGTCTTATTG‑GACT‑TGG‑GG‑TTAAGgtt                                                                                     >  2:633017/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGttt                                                                                     >  1:1053012/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGGGTA‑TCTTGCCTTACTT‑AAAT‑GGC‑GG‑TTTAGttt                                                                                     >  2:406893/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGGAGTGGGCGCAGAGAGGGGTTA‑TCTTTCCTGACTG‑AACT‑GGC‑GA‑TTTAGttt                                                                                     >  1:803220/1‑151 (MQ=255)
                                                                ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTA‑TCTTGCCTTACTG‑AACT‑GGC‑GA‑TTTAGtttt                                                                                    >  1:269131/1‑150 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgc                                                                                                                                      >  1:394209/1‑105 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCCGGGAGGGGTTT‑TCTTTCCTGGCTG‑AACT‑GGC‑GG‑TTAAAtttt                                                                                    >  2:127545/1‑149 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                      <  2:394345/105‑1 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTAGAAGGGGGGGCAGGGGGGGGTTA‑TTTTTCCTTGCTG‑TACT‑TGG‑GA‑TTTAGGTTTcg                                                                                  >  2:889720/1‑150 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAAGGAGGGCAGAGGGGGATTT‑TCTTGCCTGATTG‑AAAC‑GGG‑GA‑TTAAGTTTTcg                                                                                  >  2:687619/1‑150 (MQ=255)
                                                                        gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCTATTTTTTTTTTTTCTTCCCTGCGGAGTGAGGGCAGAGGGGGGTTA‑TCTTTGCTCAATT‑AACT‑GGG‑GA‑TTTAGTTTAcgccgac                                                                             >  1:52552/1‑146 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAGGGATTA‑TTTTGCCTGACTG‑AACT‑GGC‑GA‑TTAATGTTCcgccggcc                                                                            >  1:558378/1‑145 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGATA‑TCTTTCCTGCCTG‑AACT‑GGC‑GA‑TTAAGTTTCCCCGT‑A‑cc                                                                          >  1:368349/1‑150 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTA‑TTTTTCCTTACTT‑AACT‑GGC‑GG‑TTAAGGTTTCGCGG‑A‑cc                                                                          >  1:801274/1‑150 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                          <  2:106811/109‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑c                                                                           <  1:887735/149‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAAGGAGCGCAGAGAggggt                                                                                                                          >  1:106796/1‑107 (MQ=255)
                                                                           ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGAAATGAGGGCAGAGAGGGGTTT‑TTTTTCCTTACTG‑AACT‑GGC‑GG‑TTAAGGTTTCGCGG‑A‑ccccg                                                                       >  1:994264/1‑149 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTTTTTTTTTCCATCAGAAGTGGGGGCAGAGAGGGGTTT‑TTTTGCCTGGCTG‑AAAA‑TGC‑GATTTAATTTTCTTCGT‑‑‑CCCGGCGt                                                                    >  2:999470/1‑151 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTTCGCGT‑A‑CCCGGCGt                                                                    >  1:990648/1‑151 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGAGAAGGGAGCGCCGAGAGGGATTA‑TCTTGCCTTTCTG‑AACT‑GGC‑GG‑TTTAGGTTTCGCGG‑A‑CCCGGCgg                                                                    >  2:126908/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTA‑TCTTGTCTGACTG‑AACT‑GGC‑GA‑TTAAATTTCCGCGTCA‑CCCgcgg                                                                     >  1:286285/1‑148 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGGGGGCAGGGAGGGATTT‑TTTTGCCTGAATT‑AACA‑GGC‑GG‑TTAAGATTCCGCGG‑A‑ACCGGCgg                                                                    >  2:389529/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGGGGGGTTA‑TTTTGCCTGACTG‑AACT‑GGG‑GG‑TTAAATTTTCGCGG‑A‑CCCGGCgg                                                                    >  2:119829/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGt                                                                    <  1:278403/151‑1 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGATTTTTTTTTTTTCTGGCCTTAGAAGTGAGCGCCGGGAGGGGTTA‑TATTTCCTCTATG‑TACA‑TGC‑GA‑CTAAGTTTCCGCGT‑C‑CCCGGCGt                                                                    >  1:602177/1‑151 (MQ=255)
                                                                               aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGTATTTTTTTTTTTTCTGCCATGAGGAGGGAGCGCCGAGGGGGATTT‑TTTTTTCTTACTG‑TACT‑TGG‑GA‑TTTAGTTTTCGCGG‑A‑CCCCGGGtt                                                                   >  2:178054/1‑151 (MQ=255)
                                                                                ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGTGGGGGCAGAGAGGGGTTA‑TTTTGCCTGACTGAAACG‑GGG‑GT‑TTAAG‑TTCCGCGTAA‑CCCGGGGtt                                                                   >  2:472765/1‑151 (MQ=255)
                                                                                ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCCTGAGAAGTGAGCGGAGAGAGGGGTTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑GTAAGGTTCCGCGG‑A‑CCCGGCGttt                                                                  >  1:215174/1‑151 (MQ=255)
                                                                                 cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGGGGGCCGGGAGGGGTTATTTTTGCCTG‑CTG‑AACTGGGG‑GT‑TTAAGTTT‑CGCGG‑A‑CCCgggggttt                                                                 >  2:439531/1‑145 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGt                                                               <  1:1009451/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGt                                                               <  1:323617/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGt                                                               <  1:386632/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGt                                                               <  1:579561/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCCTTTGAAGGGAGGGCAGAGGGGGGGTT‑TTTTGTCTGTCTG‑AACT‑GGC‑GA‑TTAAGTTTTCGCGG‑ACCCCGGGGTTTat                                                                >  1:327848/1‑150 (MQ=255)
                                                                                      ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCGCGTTTGCGGCGGTCTTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGtctg                                                            <  2:189764/151‑4 (MQ=255)
                                                                                        aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGGGGGCAGGGGGGGGTTT‑TTTTTCCTGACTG‑AACT‑TGC‑GG‑TTAAGGTTTCGCGG‑A‑CCCGGCGTTTTGTTTTTg                                                          >  2:1004828/1‑151 (MQ=255)
                                                                                          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTg                                                          <  2:542221/149‑1 (MQ=255)
                                                                                          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCg                                                        <  2:796289/151‑1 (MQ=255)
                                                                                             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCAAGGGAAGTGGGCGGCGAGGGGGTTTT‑TCTTTTCCTACTT‑AACT‑TGC‑GG‑GTAAGTTTCCGCGT‑A‑CCCGGCGTTTTGTTTTTGCgggt                                                     >  2:391675/1‑149 (MQ=255)
                                                                                                    tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTA‑TTTTGTCTGTCTG‑AACC‑GGCGGT‑TTAGGTTTCCG‑GT‑A‑CCCCGCGTTTTTGTATTTCGGGTGATGata                                              >  1:907148/1‑151 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGatat                                             <  2:327947/151‑1 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGatat                                             <  1:524077/151‑1 (MQ=255)
                                                                                                          ggtGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑ACTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTg                                           <  2:143996/146‑1 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGGAGTGAGCGCCGAGAGGGGTTA‑TCTTTCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCCCGG‑A‑CCCGGCGTTTTGTTATTTCGGATGATGGTATTGAAg                                        >  1:729742/1‑150 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGGGCAGAGGGGGGTTT‑TTTTGCCTGTCTG‑AAAT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCCGCGGTTTGTTTTTGCGGGTGGTGATATTGAAg                                        >  2:795473/1‑150 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  1:732073/150‑1 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  2:505336/150‑1 (MQ=255)
                                                                                                            tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                      <  1:112581/151‑1 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgga                                  <  2:392116/151‑2 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCCTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                  <  1:112804/151‑1 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAGGGGTTA‑TCTTTCCTGACTG‑AAAT‑GGC‑GA‑TTAAGTTTTCTCGG‑A‑ACCCGCGGTTTGTTTTTTCGGATGATGATATTGAAGc                                       >  1:524215/1‑144 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                       <  2:524446/144‑1 (MQ=255)
                                                                                                                        gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCAAGAGAAGGGAGCGCAGAGGGGGGTTT‑TCTTGGCTTTCTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑ACCCGCGCTTAGTTATTTTCGGTGATGTAATTGAAACGCGGCCGATTTGt                          >  1:286289/1‑151 (MQ=255)
                                                                                                                               tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  1:472576/151‑1 (MQ=255)
                                                                                                                               tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  2:573818/151‑1 (MQ=255)
                                                                                                                                   cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc               <  2:893905/151‑1 (MQ=255)
                                                                                                                                          cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGGGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTTTTCATCACCg        >  2:458516/1‑151 (MQ=255)
                                                                                                                                                ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:779496/151‑1 (MQ=255)
                                                                                                                                                ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:761906/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTG‑AACT‑GGC‑GA‑TTAAGTTTCCGCGT‑A‑CCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604198‑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: ‑