Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT40.3% 76.9 / 31.1 63S64F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (5/32);  new base T (25/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 9.67e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  NC_000913/4604201‑4604500
                                                                                                                                                      |                                                                                                                                                        
aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttc‑tttttc                                                                                                                                                          <  2:160429/150‑1 (MQ=255)
 gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttc‑tttttct                                                                                                                                                         <  2:556253/150‑1 (MQ=255)
            cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTT‑TTTTTTTTCCATgag                                                                                                                                                 >  1:1093106/1‑147 (MQ=255)
             aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCCTGTGaa                                                                                                                                               >  1:622874/1‑148 (MQ=255)
              tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCCTGAGaa                                                                                                                                               >  2:274729/1‑147 (MQ=255)
               gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGa                                                                                                                                           <  2:174207/150‑1 (MQ=255)
                 gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTT‑TTTTTTTGCCCTGAGAAGGGAgcg                                                                                                                                        >  1:177330/1‑151 (MQ=255)
                  tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGCCATGAGAAGTg                                                                                                                                            >  1:60638/1‑146 (MQ=255)
                        cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                                 <  2:700582/151‑1 (MQ=255)
                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGCCCTGAGAAGGGGGCGCCGAGAgg                                                                                                                                >  2:343443/1‑151 (MQ=255)
                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                                <  2:942322/151‑1 (MQ=255)
                            aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCTGTTTTTT‑TTTTTATTCCCTGAGAagagaggggaggggggg                                                                                                                               >  2:1201741/1‑138 (MQ=255)
                                 gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGTCATGAGGAGGGAGCGCAGAGAGGGGTTATCt                                                                                                                        >  1:192347/1‑151 (MQ=255)
                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                                  <  1:171656/151‑1 (MQ=255)
                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                                   <  2:608074/150‑1 (MQ=255)
                                          cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                                <  2:718163/150‑1 (MQ=255)
                                            gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                             <  1:49445/151‑1 (MQ=255)
                                               aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTATTTTGCCTGACTGAACt                                                                                                          >  2:518550/1‑151 (MQ=255)
                                               aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGCCATGAGAAGGGAGCGGAGGGGGGGGTTTTTTTGCCTGGCTGaa                                                                                                            >  2:308768/1‑149 (MQ=255)
                                                  ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                                        <  1:1223345/150‑1 (MQ=255)
                                                     cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTT‑TTTTTTTGCCCTGAGGAGTGAGCGCCGGGGGGGATTATCTTTTCTGACTGAAACGGC‑g                                                                                                     >  1:1118448/1‑149 (MQ=255)
                                                      gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GAtt                                                                                                  <  2:863644/151‑1 (MQ=255)
                                                       ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTT‑TTTTTTTTCCATGAGAAGTGGGCGCCGAGAGGGGTTTTTTTGCCTGACTGAACTGGC‑GGtt                                                                                                  >  2:1221338/1‑150 (MQ=255)
                                                          aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAg                                                                                               <  1:974247/150‑1 (MQ=255)
                                                            cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGttt                                                                                            <  2:785093/151‑1 (MQ=255)
                                                                tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTcc                                                                                          <  2:396712/149‑1 (MQ=255)
                                                                tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCcg                                                                                         <  1:560132/150‑1 (MQ=255)
                                                                 ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCcgcg                                                                                       <  1:384568/151‑1 (MQ=255)
                                                                 ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAATGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCcgag                                                                                       <  1:666274/151‑3 (MQ=255)
                                                                  taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCCTGAGAAGGGAGCGCAGGGGGGGGTTTTTTTTCCTTTCTTAAATGGG‑GGTTTAGTTTcc                                                                                          >  1:298479/1‑147 (MQ=255)
                                                                   aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGCCATGAGGAGTGAGCGCCGAGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTTAAGttcc                                                                                           >  2:653560/1‑144 (MQ=255)
                                                                    aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTCTGCCCTGGGAAGTGAGCGCAGAGAGGGATTATTTTTCCTGACTGAACTGGC‑GGTTAAGTTTTCGCGT‑a                                                                                    >  2:852595/1‑150 (MQ=255)
                                                                         cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCATGGGAAGTGAGGGCAGGGAGGGGGTTTTTTTCCTGAATGCAACGGG‑GATTAAAGTTTCGCGG‑AACCgg                                                                               >  1:600124/1‑150 (MQ=255)
                                                                              cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTTTTCCCTGAGGAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAAATGGC‑GGTTAAGTTTTCGCGT‑ACCCGGCGttt                                                                          >  2:108348/1‑150 (MQ=255)
                                                                               tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCCTTAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGAATGTACTGGC‑GGTTAGGTTTTCGCGG‑AGCCCGCGttt                                                                          >  1:1241565/1‑149 (MQ=255)
                                                                                ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTa                                                                         <  2:471465/149‑1 (MQ=255)
                                                                                ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAg                                                                        <  1:1028127/150‑1 (MQ=255)
                                                                                 taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTT‑TTTTTTTCCCATGAGAAGTGAGGGCAGGGAGGGGATTTCTTGTCTGTATGAACTGGG‑GATTAAAATTCCGCGG‑ACCCGGCGTTTTGt                                                                       >  2:234682/1‑150 (MQ=255)
                                                                                  aCCTGTCAGGAGTCATTTCAGGATAGCGTGTCGCGGATCTCGTTTGCGGCGGTCAGTTTGTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATGTTGCCTGACTGAACTGGC‑TATTAAGTTTCCGCGT‑ACCCGGCGTTGAGtt                                                                      <  1:114923/150‑1 (MQ=255)
                                                                                   ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTGTCCTGAGAAGGGGGGGCAGGGAGGGGTTTTCTTGCCTTAACGAAATGGC‑GATTAAGGTTTCGCGT‑AACCCGCGCTTTGTTAt                                                                    >  1:813543/1‑151 (MQ=255)
                                                                                   ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTT‑TTTTTTTTCCATGGGGAGTGAGCGCCGGGAGGGAGTTTCTTTCCCGACTGGACTGGC‑GAGTTAATTTTCGCGT‑AGCCCGCGGTTTGtt                                                                      >  1:404853/1‑149 (MQ=255)
                                                                                    cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTAtt                                                                   <  2:122882/151‑1 (MQ=255)
                                                                                          ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTT‑TTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGC‑GATTAAGGTTTCGCGT‑TCCCGGCGTTTTGTTATTGCgg                                                               >  2:784127/1‑149 (MQ=255)
                                                                                               cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGgatgatg                                                         <  1:138563/150‑1 (MQ=255)
                                                                                               cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATATTTTCTTTTTTTTGCCTTTGGAA‑TTAGCCCCGAGGGGGAGTTAATTTCCTGACTTAACTAGC‑GATTTAGTTTTCGAGTCGCCCGGCGtttgtttttttgggt                                                              >  1:1181636/1‑134 (MQ=255)
                                                                                                   tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGatat                                                     <  2:282861/150‑1 (MQ=255)
                                                                                                    cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                    <  1:5747/150‑1 (MQ=255)
                                                                                                     aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTT‑TTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGAACTGGC‑GGTTTAGTTTTCGCGT‑AACCGGCGGTTAGTTATTTCGGATGGTGATAtt                                                    >  1:1264624/1‑149 (MQ=255)
                                                                                                       gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                <  2:109231/151‑1 (MQ=255)
                                                                                                       gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                <  2:899216/151‑1 (MQ=255)
                                                                                                        aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                               <  2:228766/151‑1 (MQ=255)
                                                                                                           ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGGGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                            <  1:319124/151‑1 (MQ=255)
                                                                                                           ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                            <  1:789124/151‑1 (MQ=255)
                                                                                                               gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                         <  1:393934/150‑1 (MQ=255)
                                                                                                                   cgGATCTCGTTTGCGGCGGTCCTTTTTTT‑TTTTTTTGCCCTGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTCCTGTCTGAACTGGG‑GGTTTAGTTTCCGCGG‑ACCCCGCGTTTTGTTTTTGCGGGTGATGATATTTAAGGTCGGCTggt                                      >  1:127501/1‑147 (MQ=255)
                                                                                                                    ggATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCCTGTGGAGAGAGAGCAGAGGGGGGATTTCCTGACCGACATAATTAGC‑GGTTTAGGTTCCCCGG‑GACCGGCGGTTAGTTATTGTGGCTGTTTAATTTGAAGccccggtggt                                      >  2:1155131/1‑139 (MQ=255)
                                                                                                                      aTCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTTTGTCATGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGAACTGGG‑GAATAAGTTTCCGCGT‑ACCCCGCGGTTTGTTATTGCGGGTGATGATATTGAAGCCCGGCTGCtt                                     >  2:1041547/1‑147 (MQ=255)
                                                                                                                          cGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                               <  1:524250/149‑1 (MQ=255)
                                                                                                                           gTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                            <  2:914267/151‑1 (MQ=255)
                                                                                                                            tttGCGGCGGTCATTTTTTT‑TTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTTCTGTCTGAACTGGC‑GATTTAATTTTCGCGG‑ACCCCGCGGTTTGTTTTTTCGGGTGGTGATATTGAAACTCGGCTGGTTGGGTTAtttt                           >  1:555924/1‑149 (MQ=255)
                                                                                                                                cggcggTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑GCCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                       >  2:1241592/1‑151 (MQ=255)
                                                                                                                                 ggcggTCATTTTTTC‑TTTTTCTGCCTTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                       >  1:544288/1‑150 (MQ=255)
                                                                                                                                 ggcggTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCa                      <  2:1251891/151‑1 (MQ=255)
                                                                                                                                  gcggTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                     <  1:899087/151‑1 (MQ=255)
                                                                                                                                   cggTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                    >  1:419817/1‑151 (MQ=255)
                                                                                                                                          tttttc‑tttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt              <  1:215322/150‑1 (MQ=255)
                                                                                                                                             ttc‑tttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc          >  2:960306/1‑151 (MQ=255)
                                                                                                                                                   tttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtg      <  2:1093172/150‑1 (MQ=255)
                                                                                                                                                    ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt     >  1:805419/1‑150 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    >  1:1069486/1‑150 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    <  1:657119/150‑1 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    <  1:472848/150‑1 (MQ=255)
                                                                                                                                                      ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAtt  >  1:825288/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTC‑TTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  NC_000913/4604201‑4604500

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