Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,686,552 T→G 5.2% R874R (CGA→CGC)  bcsC ← cellulose synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,686,5520TG 5.2% 238.2 / 3.8 77R874R (CGA→CGC) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base G (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTG  >  NC_000913/3686406‑3686699
                                                                                                                                                  |                                                                                                                                                   
aGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACAGATGCCGCCGACCACATCCACACCGTTGAAGCCCAACGGCGTGGTACCGATATCCCAGCTCCTGGCGTCATTTCg                                                                                                                                                   >  8:155510/1‑149 (MQ=255)
 gCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGc                                                                                                                                                  <  1:296224/149‑1 (MQ=255)
     aGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGc                                                                                                                                              <  2:91871/149‑1 (MQ=255)
       aaCTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGcc                                                                                                                                             >  7:198799/1‑148 (MQ=255)
          tGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGAcc                                                                                                                                         >  4:196009/1‑149 (MQ=255)
           ggAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGCCCACATCAACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCg                                                                                                                                        >  2:21351/1‑149 (MQ=255)
              gaTGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCa                                                                                                                                     <  7:95005/149‑1 (MQ=255)
                tGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCAc                                                                                                                                    <  8:198799/148‑1 (MQ=255)
                  ggCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGTCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCt                                                                                                                                 <  2:23395/149‑1 (MQ=255)
                  ggCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCt                                                                                                                                 >  1:118141/1‑149 (MQ=255)
                    ccGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCCCATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTgg                                                                                                                               >  8:110893/1‑149 (MQ=255)
                      gACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGAGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCa                                                                                                                             <  2:150979/149‑1 (MQ=255)
                           tGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGa                                                                                                                        <  2:118141/149‑1 (MQ=255)
                              gCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATc                                                                                                                     <  8:201102/149‑1 (MQ=255)
                               cGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGCCCGCCACGCTGGCACCGGAATcc                                                                                                                    >  3:8502/1‑149 (MQ=255)
                                 ttAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGa                                                                                                                  >  4:288356/1‑149 (MQ=255)
                                   aaCGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGAc                                                                                                                 <  3:288356/148‑1 (MQ=255)
                                      ggTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGCCCGCCACGCTGGCACCGGAATCCGACTGGc                                                                                                             >  3:30320/1‑149 (MQ=255)
                                        tgtAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGAc                                                                                                                 <  5:316974/143‑1 (MQ=255)
                                        tgtAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGAc                                                                                                                 >  6:316974/1‑143 (MQ=255)
                                          tAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCg                                                                                                         >  6:92109/1‑149 (MQ=255)
                                           aaCCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAg                                                                                                                                               <  4:227957/110‑1 (MQ=255)
                                            aCCCAGCGGCCCGATATCATCTCTGTAACTGATGCTGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCgg                                                                                                        <  6:105043/148‑1 (MQ=255)
                                            aCCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGc                                                                                                                                                  >  1:308106/1‑106 (MQ=255)
                                            aCCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGc                                                                                                                                                  <  2:308106/106‑1 (MQ=255)
                                                     cccGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                               >  8:57510/1‑148 (MQ=255)
                                                     cccGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTAc                                                                                              <  7:240895/149‑1 (MQ=255)
                                                      ccGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACa                                                                                             <  8:287740/149‑1 (MQ=255)
                                                      ccGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCACATCGGCGTGTTACAGATATCCCAGCTCCAGACGTCATTTCGCCAGCCTACCGCAACGCTGGCCCCGGAATCCGACTGGATGCGGTTGCCGCTACa                                                                                             >  2:64317/1‑149 (MQ=255)
                                                         atatCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                          <  4:30320/149‑1 (MQ=255)
                                                           atcatcGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCt                                                                                        >  1:210288/1‑149 (MQ=255)
                                                            tcatcGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTCCCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCTCCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCt                                                                                        <  5:12930/148‑1 (MQ=255)
                                                              atcGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGt                                                                                      >  6:150982/1‑148 (MQ=255)
                                                              atcGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGt                                                                                      >  7:297507/1‑148 (MQ=255)
                                                               tcGCTTTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTaa                                                                                    >  5:309202/1‑149 (MQ=255)
                                                                 gCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAAt                                                                                   <  5:92109/148‑1 (MQ=255)
                                                                     tAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTa                                                                                >  4:166178/1‑147 (MQ=255)
                                                                     tAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTa                                                                                >  6:10926/1‑147 (MQ=255)
                                                                      aaCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCAc                                                                                                                                    <  6:133152/94‑1 (MQ=255)
                                                                      aaCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCAc                                                                                                                                    >  5:133152/1‑94 (MQ=255)
                                                                       aCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAgg                                                                            >  3:158506/1‑149 (MQ=255)
                                                                       aCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTAGCCAGCAGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAg                                                                             <  1:64317/148‑1 (MQ=255)
                                                                        cTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGt                                                                           >  7:66671/1‑149 (MQ=255)
                                                                         tGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTCGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTg                                                                          >  3:247621/1‑149 (MQ=255)
                                                                           aTGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGCCGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGcc                                                                        >  2:377827/1‑149 (MQ=255)
                                                                               cgccgACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTAc                                                                               <  7:181535/138‑1 (MQ=255)
                                                                               cgccgACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAAACCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTAc                                                                               >  8:181535/1‑138 (MQ=255)
                                                                                gccgACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGt                                                                   >  8:48449/1‑149 (MQ=255)
                                                                                       acaTCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATc                                                            <  3:90748/149‑1 (MQ=255)
                                                                                        caTCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATcc                                                           >  4:300776/1‑149 (MQ=255)
                                                                                             accacTTTTAAGCCCATCGGCGTGGTACCGTTATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCGGACTGGCTGCGGTTTCCGCTACAGTCCTGTAATTTACAGGTGCCCCAGTTGTCATCCCAttt                                                      <  7:155510/149‑1 (MQ=255)
                                                                                             accacGTTGAAGCCCATCGGCGTGGTACCGATATCCCTGCTCAAGACGTCATTTCGCCAGCCGACCGCCCCGCTAGCACCGGATTCCGACTGGCTGCGGTTTCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTAGTCATCCCAttt                                                      >  1:317143/1‑149 (MQ=255)
                                                                                             accacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCAttt                                                      <  7:110893/149‑1 (MQ=255)
                                                                                              ccacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACa                                                                                             >  3:85581/1‑109 (MQ=255)
                                                                                              ccacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACa                                                                                             <  4:85581/109‑1 (MQ=255)
                                                                                               cacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGc                                                    >  2:147718/1‑149 (MQ=255)
                                                                                               cacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTGCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGc                                                    >  2:313753/1‑149 (MQ=255)
                                                                                               cacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTGCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGc                                                    >  6:184742/1‑149 (MQ=255)
                                                                                                     gAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAAt                                                                                   <  8:81977/112‑1 (MQ=255)
                                                                                                          ccATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCAt                                           >  1:200409/1‑147 (MQ=255)
                                                                                                          ccATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTa                                         <  5:150982/149‑1 (MQ=255)
                                                                                                           cATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGTCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAg                                        <  6:309202/149‑1 (MQ=255)
                                                                                                             tCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTg                                      <  3:196009/149‑1 (MQ=255)
                                                                                                             tCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTGCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTg                                      <  1:313753/149‑1 (MQ=255)
                                                                                                             tCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTGCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTg                                      <  5:184742/149‑1 (MQ=255)
                                                                                                                   tGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTGGTCATCCCATTTGCCATCGGCATTAGTGGAGaaa                                <  5:81253/149‑1 (MQ=255)
                                                                                                                    ggTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                               <  3:67503/85‑1 (MQ=255)
                                                                                                                    ggTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                               >  4:67503/1‑85 (MQ=255)
                                                                                                                    ggTACCGATATCCCAGCTCCAGACGTCATTTAGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAAc                               >  5:224904/1‑149 (MQ=255)
                                                                                                                     gTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAAc                               <  1:21351/148‑1 (MQ=255)
                                                                                                                     gTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACt                              >  5:267225/1‑149 (MQ=255)
                                                                                                                            tatCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACg                       >  4:206528/1‑149 (MQ=255)
                                                                                                                                  aGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCata                 <  7:57510/149‑1 (MQ=255)
                                                                                                                                       cAGACGTCATTTCGCCATCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGAc            >  3:116625/1‑149 (MQ=255)
                                                                                                                                         gACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGa          <  5:10926/149‑1 (MQ=255)
                                                                                                                                               aTTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCata                 >  5:81122/1‑136 (MQ=255)
                                                                                                                                               aTTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCata                 <  6:81122/136‑1 (MQ=255)
                                                                                                                                                tttCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACt   >  2:296160/1‑149 (MQ=255)
                                                                                                                                                  tCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTg  <  8:66671/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGTCATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCACTG  >  NC_000913/3686406‑3686699

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑