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

AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  NC_000913/4604201‑4604500
                                                                                                                                                      |                                                                                                                                                        
aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                                                                                          <  2:160429/150‑1 (MQ=255)
 gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                         <  2:556253/150‑1 (MQ=255)
            cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTTCCATgag                                                                                                                                                 >  1:1093106/1‑147 (MQ=255)
             aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGTGaa                                                                                                                                               >  1:622874/1‑148 (MQ=255)
              tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGaa                                                                                                                                               >  2:274729/1‑147 (MQ=255)
               gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                           <  2:174207/150‑1 (MQ=255)
                 gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCCTGAGAAGGGAgcg                                                                                                                                        >  1:177330/1‑151 (MQ=255)
                  tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTg                                                                                                                                            >  1:60638/1‑146 (MQ=255)
                        cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                                 <  2:700582/151‑1 (MQ=255)
                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGGGCGCCGAGAgg                                                                                                                                >  2:343443/1‑151 (MQ=255)
                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                                <  2:942322/151‑1 (MQ=255)
                            aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCTGTTTTTTTTTTTATTCCCTGAGAagagaggggaggggggg                                                                                                                               >  2:1201741/1‑138 (MQ=255)
                                 gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCATGAGGAGGGAGCGCAGAGAGGGGTTATCt                                                                                                                        >  1:192347/1‑151 (MQ=255)
                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                                  <  1:171656/151‑1 (MQ=255)
                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                                   <  2:608074/150‑1 (MQ=255)
                                          cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                                <  2:718163/150‑1 (MQ=255)
                                            gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                             <  1:49445/151‑1 (MQ=255)
                                               aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTATTTTGCCTGACTGAACt                                                                                                          >  2:518550/1‑151 (MQ=255)
                                               aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGGAGGGGGGGGTTTTTTTGCCTGGCTGaa                                                                                                            >  2:308768/1‑149 (MQ=255)
                                                  ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                                        <  1:1223345/150‑1 (MQ=255)
                                                     cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCCGGGGGGGATTATCTTTTCTGACTGAAACGGCg                                                                                                     >  1:1118448/1‑149 (MQ=255)
                                                      gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAtt                                                                                                  <  2:863644/151‑1 (MQ=255)
                                                       ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCCGAGAGGGGTTTTTTTGCCTGACTGAACTGGCGGtt                                                                                                  >  2:1221338/1‑150 (MQ=255)
                                                          aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAg                                                                                               <  1:974247/150‑1 (MQ=255)
                                                            cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                            <  2:785093/151‑1 (MQ=255)
                                                                tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                          <  2:396712/149‑1 (MQ=255)
                                                                tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcg                                                                                         <  1:560132/150‑1 (MQ=255)
                                                                 ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgcg                                                                                       <  1:384568/151‑1 (MQ=255)
                                                                 ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAATGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgag                                                                                       <  1:666274/151‑3 (MQ=255)
                                                                  taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGGGGGGTTTTTTTTCCTTTCTTAAATGGGGGTTTAGTTTcc                                                                                          >  1:298479/1‑147 (MQ=255)
                                                                   aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGTGAGCGCCGAGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTTAAGttcc                                                                                           >  2:653560/1‑144 (MQ=255)
                                                                    aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGGGAAGTGAGCGCAGAGAGGGATTATTTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTa                                                                                    >  2:852595/1‑150 (MQ=255)
                                                                         cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGGGCAGGGAGGGGGTTTTTTTCCTGAATGCAACGGGGATTAAAGTTTCGCGGAACCgg                                                                               >  1:600124/1‑150 (MQ=255)
                                                                              cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGGAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAAATGGCGGTTAAGTTTTCGCGTACCCGGCGttt                                                                          >  2:108348/1‑150 (MQ=255)
                                                                               tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGAATGTACTGGCGGTTAGGTTTTCGCGGAGCCCGCGttt                                                                          >  1:1241565/1‑149 (MQ=255)
                                                                                ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                                         <  2:471465/149‑1 (MQ=255)
                                                                                ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                        <  1:1028127/150‑1 (MQ=255)
                                                                                 taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTCCCATGAGAAGTGAGGGCAGGGAGGGGATTTCTTGTCTGTATGAACTGGGGATTAAAATTCCGCGGACCCGGCGTTTTGt                                                                       >  2:234682/1‑150 (MQ=255)
                                                                                  aCCTGTCAGGAGTCATTTCAGGATAGCGTGTCGCGGATCTCGTTTGCGGCGGTCAGTTTGTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATGTTGCCTGACTGAACTGGCTATTAAGTTTCCGCGTACCCGGCGTTGAGtt                                                                      <  1:114923/150‑1 (MQ=255)
                                                                                   ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGGGGGGGCAGGGAGGGGTTTTCTTGCCTTAACGAAATGGCGATTAAGGTTTCGCGTAACCCGCGCTTTGTTAt                                                                    >  1:813543/1‑151 (MQ=255)
                                                                                   ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTCCATGGGGAGTGAGCGCCGGGAGGGAGTTTCTTTCCCGACTGGACTGGCGAGTTAATTTTCGCGTAGCCCGCGGTTTGtt                                                                      >  1:404853/1‑149 (MQ=255)
                                                                                    cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                                   <  2:122882/151‑1 (MQ=255)
                                                                                          ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTCGCGTTCCCGGCGTTTTGTTATTGCgg                                                               >  2:784127/1‑149 (MQ=255)
                                                                                               cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                                         <  1:138563/150‑1 (MQ=255)
                                                                                               cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATATTTTCTTTTTTTTGCCTTTGGAATTAGCCCCGAGGGGGAGTTAATTTCCTGACTTAACTAGCGATTTAGTTTTCGAGTCGCCCGGCGtttgtttttttgggt                                                              >  1:1181636/1‑134 (MQ=255)
                                                                                                   tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                                     <  2:282861/150‑1 (MQ=255)
                                                                                                    cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                    <  1:5747/150‑1 (MQ=255)
                                                                                                     aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTTAGTTTTCGCGTAACCGGCGGTTAGTTATTTCGGATGGTGATAtt                                                    >  1:1264624/1‑149 (MQ=255)
                                                                                                       gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                <  2:109231/151‑1 (MQ=255)
                                                                                                       gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                <  2:899216/151‑1 (MQ=255)
                                                                                                        aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                               <  2:228766/151‑1 (MQ=255)
                                                                                                           ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGGGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                            <  1:319124/151‑1 (MQ=255)
                                                                                                           ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                            <  1:789124/151‑1 (MQ=255)
                                                                                                               gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                         <  1:393934/150‑1 (MQ=255)
                                                                                                                   cgGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTCCTGTCTGAACTGGGGGTTTAGTTTCCGCGGACCCCGCGTTTTGTTTTTGCGGGTGATGATATTTAAGGTCGGCTggt                                      >  1:127501/1‑147 (MQ=255)
                                                                                                                    ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGTGGAGAGAGAGCAGAGGGGGGATTTCCTGACCGACATAATTAGCGGTTTAGGTTCCCCGGGACCGGCGGTTAGTTATTGTGGCTGTTTAATTTGAAGccccggtggt                                      >  2:1155131/1‑139 (MQ=255)
                                                                                                                      aTCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGTCATGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGAACTGGGGAATAAGTTTCCGCGTACCCCGCGGTTTGTTATTGCGGGTGATGATATTGAAGCCCGGCTGCtt                                     >  2:1041547/1‑147 (MQ=255)
                                                                                                                          cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                               <  1:524250/149‑1 (MQ=255)
                                                                                                                           gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                            <  2:914267/151‑1 (MQ=255)
                                                                                                                            tttGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTTCTGTCTGAACTGGCGATTTAATTTTCGCGGACCCCGCGGTTTGTTTTTTCGGGTGGTGATATTGAAACTCGGCTGGTTGGGTTAtttt                           >  1:555924/1‑149 (MQ=255)
                                                                                                                                cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTGCCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                       >  2:1241592/1‑151 (MQ=255)
                                                                                                                                 ggcggTCATTTTTTCTTTTTCTGCCTTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                       >  1:544288/1‑150 (MQ=255)
                                                                                                                                 ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCa                      <  2:1251891/151‑1 (MQ=255)
                                                                                                                                  gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                     <  1:899087/151‑1 (MQ=255)
                                                                                                                                   cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                    >  1:419817/1‑151 (MQ=255)
                                                                                                                                          tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt              <  1:215322/150‑1 (MQ=255)
                                                                                                                                             ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc          >  2:960306/1‑151 (MQ=255)
                                                                                                                                                   tttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtg      <  2:1093172/150‑1 (MQ=255)
                                                                                                                                                    ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt     >  1:805419/1‑150 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    >  1:1069486/1‑150 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    <  1:657119/150‑1 (MQ=255)
                                                                                                                                                     tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa    <  1:472848/150‑1 (MQ=255)
                                                                                                                                                      ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAtt  >  1:825288/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATT  >  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: