Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT43.1% 72.8 / 38.0 72S64F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (11/30);  new base T (31/0);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 2.13e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  NC_000913/4604202‑4604497
                                                                                                                                                    |                                                                                                                                                    
gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttttttttct                                                                                                                                                     >  2:1025406/1‑150 (MQ=255)
  aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                                                                                  <  2:392638/151‑1 (MQ=255)
    aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                                                                                  >  1:523433/1‑149 (MQ=255)
         gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCACTTTTTTTTTTTTTTCCAt                                                                                                                                                >  2:987024/1‑146 (MQ=255)
          tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTCTTTTTTTTCCATgag                                                                                                                                             >  1:223090/1‑148 (MQ=255)
          tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCAt                                                                                                                                                >  1:83885/1‑145 (MQ=255)
          tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATg                                                                                                                                               >  1:1041296/1‑146 (MQ=255)
           cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGaa                                                                                                                                           >  2:690592/1‑149 (MQ=255)
           cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGaa                                                                                                                                           >  1:761266/1‑149 (MQ=255)
           cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                          <  1:407531/150‑1 (MQ=255)
            aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTg                                                                                                                                        >  2:451111/1‑151 (MQ=255)
             tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGaa                                                                                                                                           >  2:515339/1‑147 (MQ=255)
             tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTg                                                                                                                                        >  2:515519/1‑150 (MQ=255)
                  agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGGCCTGAGAAGGGa                                                                                                                                       >  1:197130/1‑146 (MQ=255)
                  agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                   <  2:273410/150‑1 (MQ=255)
                  agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                   <  1:71419/150‑1 (MQ=255)
                       cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCAAGAGAAGTGAg                                                                                                                                      >  2:117012/1‑142 (MQ=255)
                          tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                           >  2:538305/1‑150 (MQ=255)
                             tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGAtt                                                                                                                        >  1:387248/1‑150 (MQ=255)
                             tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGATTTTTTTTTTTTTTGTCATTGGGAGTGAGGGCAGAGGGGGAtt                                                                                                                        >  2:356896/1‑150 (MQ=255)
                              gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                       <  1:231911/150‑1 (MQ=255)
                                 agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                   <  2:922775/151‑1 (MQ=255)
                                 agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                   <  2:231267/151‑1 (MQ=255)
                                    ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                                 <  1:704090/150‑1 (MQ=255)
                                      aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                               <  2:708836/150‑1 (MQ=255)
                                         cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                           <  2:219891/151‑1 (MQ=255)
                                          aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactg                                                                                                          <  1:971846/151‑1 (MQ=255)
                                           gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGGTATCTTTCctgac                                                                                                            >  2:88991/1‑148 (MQ=255)
                                            aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGaa                                                                                                        <  1:881476/151‑1 (MQ=255)
                                                 ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                                    <  2:574870/150‑1 (MQ=255)
                                                   gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGGTTTTTTTTCCTGACTGaaactg                                                                                                    >  1:398805/1‑144 (MQ=255)
                                                    cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGGGAAGGGGGGGCCGAGGGGGGTTTTCTTGACTGACTGAACTGGGGAt                                                                                               >  2:819319/1‑151 (MQ=255)
                                                     gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGGGGAGTGAGCGCAGAGAGGGATTTTTTTGCCTGACTGAACTgggg                                                                                                  >  2:794718/1‑146 (MQ=255)
                                                     gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAtt                                                                                              <  1:155229/151‑1 (MQ=255)
                                                        tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGGGAGGGGTTATTTTTCCTGCCTTAAATGGCGAtt                                                                                              >  2:343665/1‑148 (MQ=255)
                                                            ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTGCATGGGAAGTGGGCGCCGAGGGGGGTTTTTTTTGCTGTCTGAACTGGCGGTTAAGgtt                                                                                        >  2:458803/1‑150 (MQ=255)
                                                             cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCGGAGAGGGGTTTTATTGCCTGACTGAACTGGGGGTTAAGgtt                                                                                        >  1:628963/1‑149 (MQ=255)
                                                             cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTCTTTTTTTTCCCCGAGGAGAGAGGGCGGTGGGGGATTATCTTTCCTGTGTGAACTTGGGATTAAGttt                                                                                        >  1:1075289/1‑149 (MQ=255)
                                                              ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGGGAGCGCGGAGGGGGATTATTTTTGCTGACTGAACTGGCGGTTAAGTTTcc                                                                                      >  2:156950/1‑150 (MQ=255)
                                                               tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                      <  2:842234/149‑1 (MQ=255)
                                                                      tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTAccc                                                                              >  1:1056156/1‑150 (MQ=255)
                                                                      tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGGGGTTAAGTTTCCGCGTAccc                                                                              >  1:598765/1‑150 (MQ=255)
                                                                       ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                            <  1:764501/151‑1 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGGCATGATAAGGGAGCGCTGAGAGGGATTATCTTTGCTGACTGAAATGGCGCTTAAGATTCCGCGTGACCGGCGttt                                                                       >  1:408143/1‑151 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGgtt                                                                       >  2:272447/1‑151 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATGTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                       <  1:946888/151‑1 (MQ=255)
                                                                               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                    <  1:688149/151‑1 (MQ=255)
                                                                                  ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGGGGGGGGTTTTTTTGTCTGAATGAACTGGCGGTTTAGTTTTCGCGTACCCCGCGGTTTGtttt                                                                 >  2:94227/1‑149 (MQ=255)
                                                                                   cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                                <  1:606207/151‑1 (MQ=255)
                                                                                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                               <  1:884032/151‑1 (MQ=255)
                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                             <  1:718902/150‑1 (MQ=255)
                                                                                        aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGGGGGGGATTTTTTTTCCTGGCTGGACTGGCGGTTAAGTTTTCGCGTACCCCGCGGTTTGTTATTGCgg                                                            >  2:151992/1‑150 (MQ=255)
                                                                                        aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGGCTGACTGAAATGGCGGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCgg                                                            >  2:878684/1‑150 (MQ=255)
                                                                                         ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                          <  2:809987/151‑1 (MQ=255)
                                                                                          gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTTCgggt                                                          >  1:470181/1‑148 (MQ=255)
                                                                                            gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTTTTTTTCCTGACTGAAATGGCGATTAAGGTTTCGCGGACCCCGCGGTTAGTTAttttgg                                                            >  2:654987/1‑142 (MQ=255)
                                                                                             tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCCTGAGAAAGGAGAGCAGAGAGGGATTATCTTGCCTTACTGAAATGGGGATTTAGTTTTCCCGTAACCGGCGGTTTGATATTGCggg                                                           >  1:455413/1‑145 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTTCTGACTGAAATGGGGATTTAGTTTTCGCGTACCCGGGGGTTTGTTTTTGCgg                                                            >  2:840565/1‑143 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                    <  1:663776/151‑1 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                    <  2:230606/151‑1 (MQ=255)
                                                                                                tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGGTTTAAGTTTCGCGGACCCGGCGGTTTGTTATTGCGggtgatgat                                                    >  1:326521/1‑150 (MQ=255)
                                                                                                   cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                 <  1:622579/150‑1 (MQ=255)
                                                                                                    aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGTGAGGGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAGGTTTCGCGGACCCCGCGGTTAGTTTTTTCGGGTGATGGTAtt                                                 >  2:393404/1‑149 (MQ=255)
                                                                                                      gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATAtt                                                 >  1:81/1‑147 (MQ=255)
                                                                                                      gATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCATTGGGCGTGAGCTGAGAGCGGGGATATTTTGCCTTTCCTCACTGGTGTTTAAGTTTTTGCGTACCCCGCGTTTTGTTATTTCGGATGATGTTAttttac                                             >  2:450412/1‑147 (MQ=255)
                                                                                                         gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                           <  2:407098/150‑1 (MQ=255)
                                                                                                                        tcGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTTAGAAGTGGGGGCCGAGGGGGATTTTTTTTTCTGACCTAACCGGGGGTTAAGGTTCCCCGGACCCCGCCTTTTGGTATTTCGGGTGATGATATTGAAGACCGGCTGATTGGtt                              >  2:759992/1‑148 (MQ=255)
                                                                                                                         cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                            <  2:93599/149‑1 (MQ=255)
                                                                                                                          gTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGGGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGGGGGTGATGATATTGAAGCTCGGCTGATTGGTTTAtt                          >  2:960448/1‑150 (MQ=255)
                                                                                                                             tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                       <  2:54513/150‑1 (MQ=255)
                                                                                                                                ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCa                   <  1:840599/151‑1 (MQ=255)
                                                                                                                                 gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                  >  1:786656/1‑151 (MQ=255)
                                                                                                                                      cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg             <  1:515358/151‑1 (MQ=255)
                                                                                                                                           tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg         <  2:1041252/150‑1 (MQ=255)
                                                                                                                                            ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGa        <  2:503028/150‑1 (MQ=255)
                                                                                                                                            ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc       >  2:749662/1‑151 (MQ=255)
                                                                                                                                                  ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt  <  2:264124/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  NC_000913/4604202‑4604497

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