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

GAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  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)
                                                    cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGGGAAGGGGGGGCCGAGGGGGGTTTTCTTGACTGACTGAACTGGG‑GAt                                                                                               >  2:819319/1‑151 (MQ=255)
                                                     gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGGGGAGTGAGCGCAGAGAGGGATTTTTTTGCCTGACTGAACTgggg                                                                                                  >  2:794718/1‑146 (MQ=255)
                                                     gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GAtt                                                                                              <  1:155229/151‑1 (MQ=255)
                                                        tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGGGAGGGGTTATTTTTCCTGCCTTAAATGGC‑GAtt                                                                                              >  2:343665/1‑148 (MQ=255)
                                                            ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTGCATGGGAAGTGGGCGCCGAGGGGGGTTTTTTTTGCTGTCTGAACTGGC‑GGTTAAGgtt                                                                                        >  2:458803/1‑150 (MQ=255)
                                                             cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCGGAGAGGGGTTTTATTGCCTGACTGAACTGGG‑GGTTAAGgtt                                                                                        >  1:628963/1‑149 (MQ=255)
                                                             cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTCTTTTTTTTCCCCGAGGAGAGAGGGCGGTGGGGGATTATCTTTCCTGTGTGAACTTGG‑GATTAAGttt                                                                                        >  1:1075289/1‑149 (MQ=255)
                                                              ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGGGAGCGCGGAGGGGGATTATTTTTGCTGACTGAACTGGC‑GGTTAAGTTTcc                                                                                      >  2:156950/1‑150 (MQ=255)
                                                               tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTcc                                                                                      <  2:842234/149‑1 (MQ=255)
                                                                      tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTTCGCGTAccc                                                                              >  1:1056156/1‑150 (MQ=255)
                                                                      tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGG‑GGTTAAGTTTCCGCGTAccc                                                                              >  1:598765/1‑150 (MQ=255)
                                                                       ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCgg                                                                            <  1:764501/151‑1 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGGCATGATAAGGGAGCGCTGAGAGGGATTATCTTTGCTGACTGAAATGGC‑GCTTAAGATTCCGCGTGACCGGCGttt                                                                       >  1:408143/1‑151 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGC‑GATTAAGTTTTCGCGTACCCGGCGgtt                                                                       >  2:272447/1‑151 (MQ=255)
                                                                            ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATGTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGttt                                                                       <  1:946888/151‑1 (MQ=255)
                                                                               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                    <  1:688149/151‑1 (MQ=255)
                                                                                  ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGGGGGGGGTTTTTTTGTCTGAATGAACTGGC‑GGTTTAGTTTTCGCGTACCCCGCGGTTTGtttt                                                                 >  2:94227/1‑149 (MQ=255)
                                                                                   cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                                <  1:606207/151‑1 (MQ=255)
                                                                                    aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                               <  1:884032/151‑1 (MQ=255)
                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                             <  1:718902/150‑1 (MQ=255)
                                                                                        aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGGGGGGGATTTTTTTTCCTGGCTGGACTGGC‑GGTTAAGTTTTCGCGTACCCCGCGGTTTGTTATTGCgg                                                            >  2:151992/1‑150 (MQ=255)
                                                                                        aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGGCTGACTGAAATGGC‑GGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCgg                                                            >  2:878684/1‑150 (MQ=255)
                                                                                         ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                          <  2:809987/151‑1 (MQ=255)
                                                                                          gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGC‑GGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTTCgggt                                                          >  1:470181/1‑148 (MQ=255)
                                                                                            gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTTTTTTTCCTGACTGAAATGGC‑GATTAAGGTTTCGCGGACCCCGCGGTTAGTTAttttgg                                                            >  2:654987/1‑142 (MQ=255)
                                                                                             tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCCTGAGAAAGGAGAGCAGAGAGGGATTATCTTGCCTTACTGAAATGGG‑GATTTAGTTTTCCCGTAACCGGCGGTTTGATATTGCggg                                                           >  1:455413/1‑145 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTTCTGACTGAAATGGG‑GATTTAGTTTTCGCGTACCCGGGGGTTTGTTTTTGCgg                                                            >  2:840565/1‑143 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                    <  1:663776/151‑1 (MQ=255)
                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                    <  2:230606/151‑1 (MQ=255)
                                                                                                tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGGTTTAAGTTT‑CGCGGACCCGGCGGTTTGTTATTGCGggtgatgat                                                    >  1:326521/1‑150 (MQ=255)
                                                                                                   cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                 <  1:622579/150‑1 (MQ=255)
                                                                                                    aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGTGAGGGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGC‑GGTTAAGGTTTCGCGGACCCCGCGGTTAGTTTTTTCGGGTGATGGTAtt                                                 >  2:393404/1‑149 (MQ=255)
                                                                                                      gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGC‑GATTTAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATAtt                                                 >  1:81/1‑147 (MQ=255)
                                                                                                      gATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCATTGGGCGTGAGCTGAGAGCGGGGATATTTTGCCTTTCCTCACTGGT‑GTTTAAGTTTTTGCGTACCCCGCGTTTTGTTATTTCGGATGATGTTAttttac                                             >  2:450412/1‑147 (MQ=255)
                                                                                                         gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                           <  2:407098/150‑1 (MQ=255)
                                                                                                                        tcGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTTAGAAGTGGGGGCCGAGGGGGATTTTTTTTTCTGACCTAACCGGG‑GGTTAAGGTTCCCCGGACCCCGCCTTTTGGTATTTCGGGTGATGATATTGAAGACCGGCTGATTGGtt                              >  2:759992/1‑148 (MQ=255)
                                                                                                                         cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                            <  2:93599/149‑1 (MQ=255)
                                                                                                                          gTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGGGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGGGGGTGATGATATTGAAGCTCGGCTGATTGGTTTAtt                          >  2:960448/1‑150 (MQ=255)
                                                                                                                             tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                       <  2:54513/150‑1 (MQ=255)
                                                                                                                                ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCa                   <  1:840599/151‑1 (MQ=255)
                                                                                                                                 gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                  >  1:786656/1‑151 (MQ=255)
                                                                                                                                      cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg             <  1:515358/151‑1 (MQ=255)
                                                                                                                                           tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg         <  2:1041252/150‑1 (MQ=255)
                                                                                                                                            ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGa        <  2:503028/150‑1 (MQ=255)
                                                                                                                                            ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc       >  2:749662/1‑151 (MQ=255)
                                                                                                                                                  ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt  <  2:264124/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGT  >  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: ‑