Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,691,3340TC49.2% 19.3 / inf 334D425G (GAC→GGC) bcsBregulator of cellulose synthase, cyclic di‑GMP binding
Reads supporting (aligned to +/- strand):  ref base T (77/92);  new base C (78/86);  total (156/178)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 335 total aligned reads displayed.

TCAGCGCCGGAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTC  >  NC_000913/3691185‑3691483
                                                                                                                                                     |                                                                                                                                                     
 cAGCGCCGGAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTc                                                                                                                                                      <  3:448783/150‑1 (MQ=255)
 cAGCGCCGGAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCt                                                                                                                                                     >  4:810209/1‑151 (MQ=255)
  aGCGCCGGAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCtt                                                                                                                                                    >  4:1691963/1‑151 (MQ=255)
         gAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACcgg                                                                                                                                             >  3:254473/1‑151 (MQ=255)
          aaTAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACcggc                                                                                                                                            <  2:443843/151‑1 (MQ=255)
           aTAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACcggcg                                                                                                                                           <  3:432873/151‑1 (MQ=255)
           aTAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACcggcg                                                                                                                                           >  1:823745/1‑151 (MQ=255)
             agagACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCggcggc                                                                                                                                         >  2:10847/1‑151 (MQ=255)
                gaCATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATg                                                                                                                                      <  3:315639/151‑1 (MQ=255)
                gaCATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATg                                                                                                                                      <  2:396705/151‑1 (MQ=255)
                 aCATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATgg                                                                                                                                     >  2:88483/1‑151 (MQ=255)
                  cATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATgg                                                                                                                                     >  1:498995/1‑150 (MQ=255)
                      tGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGgtgt                                                                                                                                  >  1:552777/1‑149 (MQ=255)
                       gTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGgtgt                                                                                                                                  <  2:539142/148‑1 (MQ=255)
                       gTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTa                                                                                                                                 <  4:1338983/149‑1 (MQ=255)
                        ttttGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGATATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCg                                                                                                                              >  3:398023/1‑151 (MQ=255)
                           tGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGt                                                                                                                            <  3:271094/150‑1 (MQ=255)
                            gCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGta                                                                                                                           <  1:1591195/150‑1 (MQ=255)
                            gCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGta                                                                                                                           <  2:1391/150‑1 (MQ=255)
                            gCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGta                                                                                                                           <  1:913460/150‑1 (MQ=255)
                             ccATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGtaat                                                                                                                         >  4:1105666/1‑151 (MQ=255)
                              cATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGtaatt                                                                                                                        <  4:475034/151‑1 (MQ=255)
                                 ccAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGtaattaat                                                                                                                     <  1:360036/151‑1 (MQ=255)
                                 ccAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGtaattaat                                                                                                                     <  4:50869/151‑1 (MQ=255)
                                   aGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAatat                                                                                                                   <  1:1920073/151‑1 (MQ=255)
                                     cAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATcc                                                                                                                 <  2:17208/151‑1 (MQ=255)
                                      aaaCCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCa                                                                                                                <  4:1362361/151‑1 (MQ=255)
                                       aaCCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCat                                                                                                               <  4:1517243/151‑1 (MQ=255)
                                       aaCCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCat                                                                                                               <  4:54362/151‑1 (MQ=255)
                                        aCCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCat                                                                                                               >  4:419318/1‑150 (MQ=255)
                                         ccTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCatat                                                                                                             <  4:921296/151‑1 (MQ=255)
                                         ccTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCatat                                                                                                             >  4:1826560/1‑151 (MQ=255)
                                             gTAATACCGGAATCCGCAGCCGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCa                                                                                                           <  1:685267/149‑1 (MQ=255)
                                             gTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAAt                                                                                                         <  1:1139161/151‑1 (MQ=255)
                                                aTACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGcc                                                                                                      <  1:5216/151‑1 (MQ=255)
                                                    cGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTa                                                                                                  <  2:422245/151‑1 (MQ=255)
                                                     ggAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTAc                                                                                                 <  4:1302274/151‑1 (MQ=255)
                                                       aaTCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTg                                                                                               <  1:1712955/151‑1 (MQ=255)
                                                       aaTCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTg                                                                                               >  1:327115/1‑151 (MQ=255)
                                                        aTCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTAc                                                                                                 >  4:574814/1‑148 (MQ=255)
                                                         tCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTgcg                                                                                             >  3:474099/1‑151 (MQ=255)
                                                            gcagcagcaGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCAt                                                                                          <  4:1188486/151‑1 (MQ=255)
                                                              agcagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATc                                                                                         >  1:1869008/1‑150 (MQ=255)
                                                              agcagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCa                                                                                        <  4:1905018/151‑1 (MQ=255)
                                                               gcagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAg                                                                                       <  4:1111344/151‑1 (MQ=255)
                                                                cagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAgg                                                                                      >  1:1744154/1‑151 (MQ=255)
                                                                 agcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGt                                                                                     >  4:1423011/1‑151 (MQ=255)
                                                                  gcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTa                                                                                    <  1:621611/151‑1 (MQ=255)
                                                                   cagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGt                                                                                     >  1:246002/1‑149 (MQ=255)
                                                                   cagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTag                                                                                   <  2:278216/151‑1 (MQ=255)
                                                                      gCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAgaga                                                                                <  1:118321/151‑1 (MQ=255)
                                                                        ggTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTag                                                                                   >  1:473834/1‑146 (MQ=255)
                                                                        ggTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAt                                                                               <  1:1835275/150‑1 (MQ=255)
                                                                           tCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcgg                                                                           >  4:364232/1‑151 (MQ=255)
                                                                            cGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggc                                                                          <  4:623556/151‑1 (MQ=255)
                                                                            cGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggc                                                                          <  2:230695/151‑1 (MQ=255)
                                                                               ctcctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCggcggc                                                                       >  1:1451628/1‑151 (MQ=255)
                                                                                 cctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCggcggc                                                                       >  4:223517/1‑149 (MQ=255)
                                                                                 cctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCa                                                                      <  3:136200/150‑1 (MQ=255)
                                                                                    gTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGtt                                                                  <  4:852283/151‑1 (MQ=255)
                                                                                        gctgctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGc                                                              <  2:512435/151‑1 (MQ=255)
                                                                                          tgctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCg                                                             >  1:1420735/1‑150 (MQ=255)
                                                                                           gctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGaa                                                           <  1:227804/151‑1 (MQ=255)
                                                                                           gctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCg                                                             <  3:371354/149‑1 (MQ=255)
                                                                                             tCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAAc                                                         >  1:1119484/1‑151 (MQ=255)
                                                                                              cAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGaaa                                                          <  4:1295756/149‑1 (MQ=255)
                                                                                                      aGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgc                                                <  1:889740/151‑1 (MQ=255)
                                                                                                        gATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctg                                              <  4:1255183/151‑1 (MQ=255)
                                                                                                         aTTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctgc                                             >  4:656786/1‑151 (MQ=255)
                                                                                                           tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctgct                                            <  4:256725/150‑1 (MQ=255)
                                                                                                                gAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCa                                      <  1:1222540/151‑1 (MQ=255)
                                                                                                                  aCTGGTTATTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCa                                      >  4:1562116/1‑149 (MQ=255)
                                                                                                                         aTTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCt                              <  1:404059/150‑1 (MQ=255)
                                                                                                                         aTTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                             >  4:1458160/1‑151 (MQ=255)
                                                                                                                         aTTCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCt                              <  4:736835/150‑1 (MQ=255)
                                                                                                                          ttCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTgg                            >  1:60317/1‑151 (MQ=255)
                                                                                                                          ttCAGGCTGATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTCAATCGCTGCTGGCTCAAGACCGCTgg                            >  3:58956/1‑151 (MQ=255)
                                                                                                                            cAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGa                           <  4:934853/150‑1 (MQ=255)
                                                                                                                                 tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                      >  4:1792406/1‑150 (MQ=255)
                                                                                                                                  gATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAAc                    <  4:888055/151‑1 (MQ=255)
                                                                                                                                  gATATCCATCCGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                      <  1:935545/149‑1 (MQ=255)
                                                                                                                                   atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACt                   <  2:327292/151‑1 (MQ=255)
                                                                                                                                   atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACt                   >  1:78256/1‑151 (MQ=255)
                                                                                                                                   atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACt                   <  1:1543492/151‑1 (MQ=255)
                                                                                                                                     atccatccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGt                 <  1:1092521/151‑1 (MQ=255)
                                                                                                                                     atccatccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGt                 <  4:105362/151‑1 (MQ=255)
                                                                                                                                     atccatccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGt                 <  3:88483/151‑1 (MQ=255)
                                                                                                                                         atccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGttctt             <  4:23756/151‑1 (MQ=255)
                                                                                                                                           ccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGttcttc            <  4:1081038/150‑1 (MQ=255)
                                                                                                                                           ccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGttcttc            <  1:1386747/150‑1 (MQ=255)
                                                                                                                                           ccGCGAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCa           <  2:254480/151‑1 (MQ=255)
                                                                                                                                            cgcgAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCAt          >  3:197513/1‑151 (MQ=255)
                                                                                                                                            cgcgAACTGCCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCAt          >  1:1043888/1‑151 (MQ=255)
                                                                                                                                              cgAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAg        <  4:1152371/151‑1 (MQ=255)
                                                                                                                                                aaCTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGt      <  1:723432/151‑1 (MQ=255)
                                                                                                                                                aaCTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGt      <  1:1657758/151‑1 (MQ=255)
                                                                                                                                                  cTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGttt    >  4:1875742/1‑151 (MQ=255)
                                                                                                                                                   tGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGtttt   >  1:1722191/1‑151 (MQ=255)
                                                                                                                                                    gTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGtttt   <  4:1928302/150‑1 (MQ=255)
                                                                                                                                                    gTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGtttt   <  3:173530/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGCGCCGGAATAGAGACATCTGTTTTGCCATCCAGCAAACCTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTC  >  NC_000913/3691185‑3691483

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