Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,818 A→T 5.6% L111Q (CTG→CAG)  algK ← alginate biosynthesis protein AlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,8180AT 5.6% 266.0 / 5.9 90L111Q (CTG→CAG) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base T (3/2);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCG  >  NC_002947/1469669‑1469964
                                                                                                                                                     |                                                                                                                                                  
gcgcGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGccacca                                                                                                                                                    <  2:840316/150‑1 (MQ=255)
  gcgTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACcagc                                                                                                                                                  >  2:1254680/1‑150 (MQ=255)
       gACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggcg                                                                                                                                             <  1:928630/150‑1 (MQ=255)
          cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGcc                                                                                                                                          >  2:74651/1‑150 (MQ=255)
                   aCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgcgcgc                                                                                                                                 <  2:1026549/150‑1 (MQ=255)
                     tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                               >  2:757173/1‑150 (MQ=255)
                     tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                               <  1:1436759/150‑1 (MQ=255)
                     tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                               >  1:403618/1‑150 (MQ=255)
                     tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCGTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                               >  2:1821826/1‑150 (MQ=255)
                     tatGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                               <  1:74651/148‑1 (MQ=255)
                       tGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTgg                                                                                                                             >  2:33496/1‑150 (MQ=255)
                       tGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTAGGACACCAGCAGGCGGCCAAGGCGCGCTTgg                                                                                                                             <  2:1137619/150‑1 (MQ=255)
                         ggggTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGGATTGCTTTGGACCTGTTTCGTGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGAGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGgc                                                                                                                           <  2:979351/149‑1 (MQ=255)
                          gggTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGgcg                                                                                                                          >  2:10463/1‑150 (MQ=255)
                           ggTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGgcgc                                                                                                                         >  2:1673491/1‑150 (MQ=255)
                                 tGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCg                                                                                                                   >  2:1750415/1‑150 (MQ=255)
                                     gTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGt                                                                                                               <  2:898303/150‑1 (MQ=255)
                                       agagCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTgg                                                                                                             >  1:1743325/1‑150 (MQ=255)
                                       agagCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTgg                                                                                                             <  1:1818571/150‑1 (MQ=255)
                                           cATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGACTAGCCACCATCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTggcggc                                                                                                         >  2:528660/1‑150 (MQ=255)
                                                ccAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAg                                                                                                    <  2:1890142/150‑1 (MQ=255)
                                                 cAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGAc                                                                                                                 >  1:1399307/1‑136 (MQ=255)
                                                 cAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGAc                                                                                                                 <  2:1399307/136‑1 (MQ=255)
                                                  aGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcg                                                                                                  <  1:1762019/150‑1 (MQ=255)
                                                   gCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcgc                                                                                                 >  1:1687327/1‑150 (MQ=255)
                                                   gCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcgc                                                                                                 <  1:439298/150‑1 (MQ=255)
                                                      gAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGAt                                                                                              >  2:486303/1‑150 (MQ=255)
                                                       aaTCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCGTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATa                                                                                             <  2:796352/150‑1 (MQ=255)
                                                        aTCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAg                                                                                            <  2:352047/150‑1 (MQ=255)
                                                         tCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAgg                                                                                           >  2:866503/1‑150 (MQ=255)
                                                          cAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGt                                                                                          <  1:1801074/150‑1 (MQ=255)
                                                            gggTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCTTGCTTGGCAACAAACAGCAGGACACCACTCACGGGGTCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTgg                                                                                        >  1:425054/1‑150 (MQ=255)
                                                            gggTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTgg                                                                                        <  1:1777748/150‑1 (MQ=255)
                                                             ggTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGc                                                                                       >  1:847538/1‑150 (MQ=255)
                                                                aTTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                    >  1:569758/1‑150 (MQ=255)
                                                                   gCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTTGCGCCCGCAGCGCGATAGGTGGCCTCGGc                                                                                 <  1:528660/150‑1 (MQ=255)
                                                                   gCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGc                                                                                 <  2:403618/150‑1 (MQ=255)
                                                                     tttGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCtt                                                                               <  1:902062/150‑1 (MQ=255)
                                                                      ttGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCAACAAGCAGGCGGACAAGTTGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTc                                                                              <  1:1673491/150‑1 (MQ=255)
                                                                          cccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCGCCTGCTGGTGGCCAAgccg                                                                                                                                    >  2:56904/1‑89 (MQ=37)
                                                                          cccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCGCCTGCTGGTGGCCAAGccg                                                                                                                                    <  1:56904/92‑1 (MQ=38)
                                                                           ccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGa                                                                         >  2:1632745/1‑150 (MQ=255)
                                                                            cTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGAt                                                                        <  1:164728/150‑1 (MQ=255)
                                                                             tGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGAtt                                                                       >  2:856055/1‑150 (MQ=255)
                                                                             tGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGAtt                                                                       <  1:1750415/150‑1 (MQ=255)
                                                                              gTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGAAGCAATCGCGTCCTTGAttt                                                                      >  2:1295628/1‑150 (MQ=255)
                                                                                ttCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGGGTCGGCTTCCTTGATTTGc                                                                    <  2:274344/150‑1 (MQ=255)
                                                                                ttCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGc                                                                    <  2:41285/150‑1 (MQ=255)
                                                                                       ccTGCTTGAGCCGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGt                                                             <  1:856055/150‑1 (MQ=255)
                                                                                       ccTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGt                                                             >  2:1511460/1‑150 (MQ=255)
                                                                                        cTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTc                                                            >  1:167280/1‑150 (MQ=255)
                                                                                         tGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                    >  2:684787/1‑125 (MQ=255)
                                                                                         tGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                    <  1:684787/125‑1 (MQ=255)
                                                                                          gCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgc                                                          <  2:847538/150‑1 (MQ=255)
                                                                                           cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                         >  2:1713385/1‑150 (MQ=255)
                                                                                           cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                         >  2:1355023/1‑150 (MQ=255)
                                                                                               aGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                     <  1:211844/150‑1 (MQ=255)
                                                                                               aGCAGGGTTTCGGCTTCCGCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                     <  2:412084/150‑1 (MQ=255)
                                                                                                gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTCAAGCTCGTGGCGACGTGGAAGCCGTACAGCGAAAGGTGGCCGCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTgg                                                    >  1:888969/1‑150 (MQ=255)
                                                                                                 cAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGc                                                   >  2:1431538/1‑150 (MQ=255)
                                                                                                  aGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGcc                                                  >  1:67136/1‑150 (MQ=255)
                                                                                                        tCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTg                                            >  2:264996/1‑150 (MQ=255)
                                                                                                           gCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATg                                         >  1:403010/1‑150 (MQ=255)
                                                                                                                 tCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGcc                                   >  2:1301142/1‑150 (MQ=255)
                                                                                                                 tCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGcc                                   >  2:1167491/1‑150 (MQ=255)
                                                                                                                   aCGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAg                                 >  1:671980/1‑150 (MQ=255)
                                                                                                                     gTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                               <  1:1511460/150‑1 (MQ=255)
                                                                                                                     gTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                               <  1:1355023/150‑1 (MQ=255)
                                                                                                                          gCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACt                          >  1:625211/1‑150 (MQ=255)
                                                                                                                            cTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCAcc                                               <  1:1082059/127‑1 (MQ=255)
                                                                                                                            cTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCAcc                                               >  2:1082059/1‑127 (MQ=255)
                                                                                                                                  gCTGCCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTc                  <  2:707019/150‑1 (MQ=255)
                                                                                                                                    tGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcg                                                                                                  >  1:680834/1‑68 (MQ=255)
                                                                                                                                    tGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcg                                                                                                  <  2:680834/68‑1 (MQ=255)
                                                                                                                                    tGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGc                >  1:171790/1‑150 (MQ=255)
                                                                                                                                    tGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGc                >  1:564578/1‑150 (MQ=255)
                                                                                                                                      tacGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGGGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTg              >  1:892523/3‑150 (MQ=255)
                                                                                                                                       gcgccTTGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                         >  2:1637348/5‑106 (MQ=38)
                                                                                                                                       gcgccTTGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                         <  1:1637348/102‑1 (MQ=38)
                                                                                                                                         ggCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGAt                                          >  1:1411425/1‑119 (MQ=255)
                                                                                                                                         ggCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGAt                                          <  2:1411425/119‑1 (MQ=255)
                                                                                                                                          gCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCgg                                                                                  >  1:1043232/1‑78 (MQ=255)
                                                                                                                                          gCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCgg                                                                                  <  2:1043232/78‑1 (MQ=255)
                                                                                                                                          ccTTGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGc          >  2:1691191/2‑150 (MQ=255)
                                                                                                                                            ttGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAg                                                                                            <  2:1309995/66‑1 (MQ=255)
                                                                                                                                            ttGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAg                                                                                            >  1:1309995/1‑66 (MQ=255)
                                                                                                                                              ggCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                      >  1:160134/1‑118 (MQ=255)
                                                                                                                                              ggCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                      <  2:160134/118‑1 (MQ=255)
                                                                                                                                              ggCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAg      <  1:757173/150‑1 (MQ=255)
                                                                                                                                              ggCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAg      >  1:1840752/1‑150 (MQ=255)
                                                                                                                                               gCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAgg     <  1:886150/150‑1 (MQ=255)
                                                                                                                                               gCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCAGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCGGTAGCCCacg     >  1:72334/1‑148 (MQ=255)
                                                                                                                                                 caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTc   <  2:1064706/150‑1 (MQ=255)
                                                                                                                                                 caccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTc   >  1:733925/1‑150 (MQ=255)
                                                                                                                                                  accaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCg  >  2:1537075/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGCGTTGACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCG  >  NC_002947/1469669‑1469964

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