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

CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  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)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGGGCAGAGGGGGATTAt                                                                                                                      >  2:718924/1‑149 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCGTTTTTTTTTTTTTTGCCCCGCTGAAAGTGCGCAGAGAGGGGTTAt                                                                                                                      >  1:842715/1‑149 (MQ=255)
                                     agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                     <  2:903883/149‑1 (MQ=255)
                                      gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGTCATGAGGAGTGAGGGCAGAGAGGGATTATTTTg                                                                                                                  >  1:245213/1‑151 (MQ=255)
                                         tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAGGGAGGGGAGGGGGGGAGTTTCtttc                                                                                                                 >  2:567654/1‑147 (MQ=255)
                                                 aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                                                                                      <  1:718664/151‑1 (MQ=255)
                                                        cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGGGGGGGGTTATCTTGCCTGACTGAACTGGCggt                                                                                             >  2:1002121/1‑149 (MQ=255)
                                                         gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCCTGAGAAGTGAGAGCAGAGAGGGGTTTTCTTTCCTGTCTGAACTGGCgggt                                                                                            >  1:21135/1‑148 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCCGAGAGGGGTTTTTTTTCCTTACTGAAATGGCGGTTAAGgtt                                                                                     >  2:660347/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGGGAAGTGTGCGCTGAGAGGGGTTTTCTTGCCTTACTGAACTTGCGGTTAAGttt                                                                                     >  2:381787/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGGGGAGGGAGCGCCGGGGGGGGTTTTTTTGTCTTATTGGACTTGGGGTTAAGgtt                                                                                     >  2:633017/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                     >  1:1053012/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGGGTATCTTGCCTTACTTAAATGGCGGTTTAGttt                                                                                     >  2:406893/1‑151 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGGAGTGGGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGttt                                                                                     >  1:803220/1‑151 (MQ=255)
                                                                ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTTACTGAACTGGCGATTTAGtttt                                                                                    >  1:269131/1‑150 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgc                                                                                                                                      >  1:394209/1‑105 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCCGGGAGGGGTTTTCTTTCCTGGCTGAACTGGCGGTTAAAtttt                                                                                    >  2:127545/1‑149 (MQ=255)
                                                                 cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                      <  2:394345/105‑1 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTAGAAGGGGGGGCAGGGGGGGGTTATTTTTCCTTGCTGTACTTGGGATTTAGGTTTcg                                                                                  >  2:889720/1‑150 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAAGGAGGGCAGAGGGGGATTTTCTTGCCTGATTGAAACGGGGATTAAGTTTTcg                                                                                  >  2:687619/1‑150 (MQ=255)
                                                                        gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCTATTTTTTTTTTTTCTTCCCTGCGGAGTGAGGGCAGAGGGGGGTTATCTTTGCTCAATTAACTGGGGATTTAGTTTAcgccgac                                                                             >  1:52552/1‑146 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAGGGATTATTTTGCCTGACTGAACTGGCGATTAATGTTCcgccggcc                                                                            >  1:558378/1‑145 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGATATCTTTCCTGCCTGAACTGGCGATTAAGTTTCCCCGTAcc                                                                          >  1:368349/1‑150 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTTACTTAACTGGCGGTTAAGGTTTCGCGGAcc                                                                          >  1:801274/1‑150 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                          <  2:106811/109‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAc                                                                           <  1:887735/149‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAAGGAGCGCAGAGAggggt                                                                                                                          >  1:106796/1‑107 (MQ=255)
                                                                           ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGAAATGAGGGCAGAGAGGGGTTTTTTTTCCTTACTGAACTGGCGGTTAAGGTTTCGCGGAccccg                                                                       >  1:994264/1‑149 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTTTTTTTTTCCATCAGAAGTGGGGGCAGAGAGGGGTTTTTTTGCCTGGCTGAAAATGCGATTTAATTTTCTTCGTCCCGGCGt                                                                    >  2:999470/1‑151 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGt                                                                    >  1:990648/1‑151 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGAGAAGGGAGCGCCGAGAGGGATTATCTTGCCTTTCTGAACTGGCGGTTTAGGTTTCGCGGACCCGGCgg                                                                    >  2:126908/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTATCTTGTCTGACTGAACTGGCGATTAAATTTCCGCGTCACCCgcgg                                                                     >  1:286285/1‑148 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGGGGGCAGGGAGGGATTTTTTTGCCTGAATTAACAGGCGGTTAAGATTCCGCGGAACCGGCgg                                                                    >  2:389529/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGGGGTTAAATTTTCGCGGACCCGGCgg                                                                    >  2:119829/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt                                                                    <  1:278403/151‑1 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGATTTTTTTTTTTTCTGGCCTTAGAAGTGAGCGCCGGGAGGGGTTATATTTCCTCTATGTACATGCGACTAAGTTTCCGCGTCCCCGGCGt                                                                    >  1:602177/1‑151 (MQ=255)
                                                                               aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGTATTTTTTTTTTTTCTGCCATGAGGAGGGAGCGCCGAGGGGGATTTTTTTTTCTTACTGTACTTGGGATTTAGTTTTCGCGGACCCCGGGtt                                                                   >  2:178054/1‑151 (MQ=255)
                                                                                ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGTGGGGGCAGAGAGGGGTTATTTTGCCTGACTGAAACGGGGGTTTAAGTTCCGCGTAACCCGGGGtt                                                                   >  2:472765/1‑151 (MQ=255)
                                                                                ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCCTGAGAAGTGAGCGGAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGAGTAAGGTTCCGCGGACCCGGCGttt                                                                  >  1:215174/1‑151 (MQ=255)
                                                                                 cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGGGGGCCGGGAGGGGTTATTTTTGCCTGCTGAACTGGGGGTTTAAGTTTCGCGGACCCgggggttt                                                                 >  2:439531/1‑145 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                               <  1:1009451/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                               <  1:323617/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                               <  1:386632/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                               <  1:579561/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCCTTTGAAGGGAGGGCAGAGGGGGGGTTTTTTGTCTGTCTGAACTGGCGATTAAGTTTTCGCGGACCCCGGGGTTTat                                                                >  1:327848/1‑150 (MQ=255)
                                                                                      ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCGCGTTTGCGGCGGTCTTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtctg                                                            <  2:189764/151‑4 (MQ=255)
                                                                                        aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGGGGGCAGGGGGGGGTTTTTTTTCCTGACTGAACTTGCGGTTAAGGTTTCGCGGACCCGGCGTTTTGTTTTTg                                                          >  2:1004828/1‑151 (MQ=255)
                                                                                          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                          <  2:542221/149‑1 (MQ=255)
                                                                                          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                        <  2:796289/151‑1 (MQ=255)
                                                                                             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCAAGGGAAGTGGGCGGCGAGGGGGTTTTTCTTTTCCTACTTAACTTGCGGGTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCgggt                                                     >  2:391675/1‑149 (MQ=255)
                                                                                                    tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGTCTGTCTGAACCGGCGGTTTAGGTTTCCGGTACCCCGCGTTTTTGTATTTCGGGTGATGata                                              >  1:907148/1‑151 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                             <  2:327947/151‑1 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                             <  1:524077/151‑1 (MQ=255)
                                                                                                          ggtGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAACTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                           <  2:143996/146‑1 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGGAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCCCGGACCCGGCGTTTTGTTATTTCGGATGATGGTATTGAAg                                        >  1:729742/1‑150 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGGGCAGAGGGGGGTTTTTTTGCCTGTCTGAAATGGCGATTAAGTTTCCGCGTACCCCGCGGTTTGTTTTTGCGGGTGGTGATATTGAAg                                        >  2:795473/1‑150 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  1:732073/150‑1 (MQ=255)
                                                                                                           aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  2:505336/150‑1 (MQ=255)
                                                                                                            tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                      <  1:112581/151‑1 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgga                                  <  2:392116/151‑2 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCCTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                  <  1:112804/151‑1 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAGGGGTTATCTTTCCTGACTGAAATGGCGATTAAGTTTTCTCGGAACCCGCGGTTTGTTTTTTCGGATGATGATATTGAAGc                                       >  1:524215/1‑144 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                       <  2:524446/144‑1 (MQ=255)
                                                                                                                        gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCAAGAGAAGGGAGCGCAGAGGGGGGTTTTCTTGGCTTTCTGAACTGGCGATTAAGTTTCCGCGTAACCCGCGCTTAGTTATTTTCGGTGATGTAATTGAAACGCGGCCGATTTGt                          >  1:286289/1‑151 (MQ=255)
                                                                                                                               tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  1:472576/151‑1 (MQ=255)
                                                                                                                               tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  2:573818/151‑1 (MQ=255)
                                                                                                                                   cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc               <  2:893905/151‑1 (MQ=255)
                                                                                                                                          cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGGGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTTTTCATCACCg        >  2:458516/1‑151 (MQ=255)
                                                                                                                                                ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:779496/151‑1 (MQ=255)
                                                                                                                                                ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:761906/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  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: