Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,821 A→T 5.4% L110Q (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,8210AT 5.4% 292.5 / 4.8 93L110Q (CTG→CAG) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base A (47/41);  new base T (3/2);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGC  >  NC_002947/1469676‑1469970
                                                                                                                                                 |                                                                                                                                                     
gACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAggcg                                                                                                                                                   <  1:928630/150‑1 (MQ=255)
   cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGcc                                                                                                                                                >  2:74651/1‑150 (MQ=255)
            aCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgcgcgc                                                                                                                                       <  2:1026549/150‑1 (MQ=255)
              tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                                     >  1:403618/1‑150 (MQ=255)
              tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                                     >  2:757173/1‑150 (MQ=255)
              tctGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCtt                                                                                                                                     <  1:1436759/150‑1 (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                                                                                                                       <  2:1399307/136‑1 (MQ=255)
                                          cAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGAc                                                                                                                       >  1:1399307/1‑136 (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                                                                                          <  1:684787/125‑1 (MQ=255)
                                                                                  tGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                          >  2:684787/1‑125 (MQ=255)
                                                                                   gCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgc                                                                <  2:847538/150‑1 (MQ=255)
                                                                                    cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                               >  2:1355023/1‑150 (MQ=255)
                                                                                    cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgcg                                                               >  2:1713385/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:1167491/1‑150 (MQ=255)
                                                                                                          tCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGcc                                         >  2:1301142/1‑150 (MQ=255)
                                                                                                            aCGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAg                                       >  1:671980/1‑150 (MQ=255)
                                                                                                              gTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                                     <  1:1355023/150‑1 (MQ=255)
                                                                                                              gTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                                     <  1:1511460/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                                                                                        <  2:1043232/78‑1 (MQ=255)
                                                                                                                                   gCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCgg                                                                                        >  1:1043232/1‑78 (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)
                                                                                                                                               gcagGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCa    >  1:1130815/1‑150 (MQ=255)
                                                                                                                                                 agGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAgc  >  1:1231759/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACCTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGC  >  NC_002947/1469676‑1469970

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: