Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,691,406 T→G D401A (GAT→GCT)  bcsB ← regulator of cellulose synthase, cyclic di‑GMP binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,691,4060TG98.7% 276.1 / ‑4.2 74D401A (GAT→GCT) bcsBregulator of cellulose synthase, cyclic di‑GMP binding
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (43/30);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

CCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTACGCGG  >  NC_000913/3691263‑3691548
                                                                                                                                               |                                                                                                                                              
cctcctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCcggcgg                                                                                                                                         <  2:923012/151‑1 (MQ=255)
cctcctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCcggcgg                                                                                                                                         >  1:393979/1‑151 (MQ=255)
cctcctGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCcggcgg                                                                                                                                         >  2:427066/1‑151 (MQ=255)
 ctcctGTTTGCTGCTCAGGGTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCTTCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCggcggc                                                                                                                                        >  1:10537/1‑151 (MQ=255)
      gTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAgg                                                                                                                                     <  1:262755/149‑1 (MQ=255)
         tgctgcTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTa                                                                                                                                 <  1:564951/150‑1 (MQ=255)
          gctgctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAg                                                                                                                                <  2:670646/150‑1 (MQ=255)
                 aGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGt                                                                                                                        <  2:188837/151‑1 (MQ=255)
                  ggTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGtt                                                                                                                       <  1:835307/151‑1 (MQ=255)
                     tGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAAt                                                                                                                    >  1:12458/1‑151 (MQ=255)
                     tGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAAt                                                                                                                    >  2:277261/1‑151 (MQ=255)
                       aaGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCg                                                                                                                  >  2:1193909/1‑151 (MQ=255)
                        aGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATACCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCgc                                                                                                                 >  1:763137/1‑151 (MQ=255)
                               cAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCTGCGGCATGGGGTTGTGGTAATTAATATCCATATCAATGCCGGTACTGCCCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGCAACGCTAATCGCTGCTgg                                                                                                           >  1:1176037/1‑150 (MQ=255)
                               cAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGc                                                                                                          >  1:1213421/1‑151 (MQ=255)
                                aGGAACTGGTTATTCAGGCCGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGctgctg                                                                                                            >  1:1199871/1‑148 (MQ=255)
                                    aCTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAg                                                                                                     <  2:1208545/151‑1 (MQ=255)
                                      tGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAc                                                                                                   >  2:1263270/1‑151 (MQ=255)
                                      tGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAc                                                                                                   >  2:1029522/1‑151 (MQ=255)
                                       ggTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAcc                                                                                                  <  1:620406/151‑1 (MQ=255)
                                        gTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCg                                                                                                 <  2:1086663/151‑1 (MQ=255)
                                           aTTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                                                                                              >  2:1153012/1‑151 (MQ=255)
                                            ttCAGCCTGATATCCATCCGCGAACTGTCTTTCAACGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGAGCATCAGGCAGAGAGCCGGCGGCAGGTTTAGAGAAACGTTAATCGCTGCTGGCTCAAGACCGCt                                                                                               >  1:1209424/1‑149 (MQ=255)
                                             tCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTgg                                                                                             <  2:850991/150‑1 (MQ=255)
                                              cAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGAt                                                                                           <  1:992117/151‑1 (MQ=255)
                                                  cTGATATCCATCCGCGAACTGTCTTTCCCCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTgg                                                                                             >  1:198878/1‑145 (MQ=255)
                                                   tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                                                                                       >  1:635957/1‑150 (MQ=255)
                                                   tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                                                                                       >  2:1070015/1‑150 (MQ=255)
                                                   tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                                                                                       >  2:650335/1‑150 (MQ=255)
                                                    gATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTa                                                                                       <  2:229370/149‑1 (MQ=255)
                                                            tccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGttcttc                                                                             >  1:1069121/1‑151 (MQ=255)
                                                               gcgAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAAATGTTCTTCAt                                                                           >  1:631116/1‑150 (MQ=255)
                                                                 gAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAgg                                                                        <  2:128237/151‑1 (MQ=255)
                                                                   aCTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTACCTGTTCGTCATAGGt                                                                       <  1:664873/150‑1 (MQ=255)
                                                                      gTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGtttt                                                                    >  2:520695/1‑150 (MQ=255)
                                                                      gTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGtttt                                                                    <  1:1032936/150‑1 (MQ=255)
                                                                        cTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTc                                                                   >  1:217067/1‑149 (MQ=255)
                                                                          ttCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAg                                                                 <  2:1162588/149‑1 (MQ=255)
                                                                           tCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGtt                                                               <  2:732889/150‑1 (MQ=255)
                                                                                   ggcATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAAtg                                                      <  1:110009/151‑1 (MQ=255)
                                                                                   ggcATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAAtg                                                      >  2:468255/1‑151 (MQ=255)
                                                                                    gcATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAAtgt                                                     >  2:542655/1‑151 (MQ=255)
                                                                                    gcATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAAtgt                                                     <  2:172068/151‑1 (MQ=255)
                                                                                          tgtAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCgg                                               >  1:878984/1‑151 (MQ=255)
                                                                                             aGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACg                                            >  2:534328/1‑151 (MQ=255)
                                                                                             aGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACg                                            >  1:723902/1‑151 (MQ=255)
                                                                                               cGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGAt                                          >  1:59861/1‑151 (MQ=255)
                                                                                               cGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGAt                                          <  2:735393/151‑1 (MQ=255)
                                                                                                ggTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATc                                         <  2:1170927/151‑1 (MQ=255)
                                                                                                 gcaattaatATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCg                                        >  2:606520/3‑151 (MQ=255)
                                                                                                 gTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCg                                        >  1:516224/1‑151 (MQ=255)
                                                                                                  taattaatATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCgg                                       >  2:760965/1‑151 (MQ=255)
                                                                                                    attaatATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGgt                                      >  1:1125152/1‑150 (MQ=255)
                                                                                                      taatATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGgtac                                    >  1:531151/1‑150 (MQ=255)
                                                                                                            ccATATCAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTAccc                              >  1:18329/1‑150 (MQ=255)
                                                                                                            ccATATCAATACCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGATCGCCAAATGCGACCGGACGATCGGTACGTAcc                               >  2:472907/1‑149 (MQ=255)
                                                                                                                atCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTACTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTc                         >  2:482135/1‑151 (MQ=255)
                                                                                                                  cAATGCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCgg                       <  1:795492/151‑1 (MQ=255)
                                                                                                                      gCCGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGgcgc                    >  2:838852/1‑150 (MQ=255)
                                                                                                                        cGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATc                 >  2:29911/1‑151 (MQ=255)
                                                                                                                        cGGTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATc                 >  2:1013976/1‑151 (MQ=255)
                                                                                                                         ggTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCg                <  2:1035269/151‑1 (MQ=255)
                                                                                                                          gTACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCTAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGt               <  1:1021120/151‑1 (MQ=255)
                                                                                                                           tACTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGt               >  2:508736/1‑150 (MQ=255)
                                                                                                                             cTGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCCGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTAc             <  1:96878/150‑1 (MQ=255)
                                                                                                                              tGCGCATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACgg           >  2:1234183/1‑151 (MQ=255)
                                                                                                                               gcgcATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACg            <  1:538978/149‑1 (MQ=255)
                                                                                                                                cgcATCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCt         <  2:855676/151‑1 (MQ=255)
                                                                                                                                   aTCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTAc      >  2:35952/1‑151 (MQ=255)
                                                                                                                                   aTCAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTAc      <  1:427056/151‑1 (MQ=255)
                                                                                                                                     cAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTAcgc    <  2:393989/151‑1 (MQ=255)
                                                                                                                                     cAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTAcgc    >  1:85505/1‑151 (MQ=255)
                                                                                                                                     cAGGTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTAcgc    <  2:1130091/151‑1 (MQ=255)
                                                                                                                                       ggTAGAGAGCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTACGCgg  >  1:382741/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGCCAAATGTGACCGGACGATCGGTACGTACCCAGTTCGGCGCATCGTACGGCTTACGCGG  >  NC_000913/3691263‑3691548

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