Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,462,291 G→C 8.0% A214G (GCC→GGC)  algI ← alginate O‑acetylase AlgI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,462,2910GC 8.0% 234.7 / 9.5 75A214G (GCC→GGC) algIalginate O‑acetylase AlgI
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (2/4);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TCCATGAAGCGGAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAGCGCAGCACGGGGCCTGCGA  >  NC_002947/1462144‑1462438
                                                                                                                                                   |                                                                                                                                                   
tCCATGAAGCGGAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCg                                                                                                                                                   <  1:1012858/150‑1 (MQ=255)
  cATGAAGCGGAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGAc                                                                                                                                                 >  2:1849680/1‑150 (MQ=255)
    tGAAGCGGAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCa                                                                                                                                               >  2:461269/1‑150 (MQ=255)
        gCGGAAGCCCATCATCAGGCCCAGGCCGATGGCAATGTCGCTGTAGCCGTTGAAGTCGATGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGCCGACCACTGc                                                                                                                                           <  1:1365640/150‑1 (MQ=255)
           gAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAg                                                                                                                                        <  1:317287/150‑1 (MQ=255)
           gAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGCCGAacactggttc                                                                                                                                        <  2:1902504/150‑7 (MQ=255)
              gCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGGGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGGGGTCGGCGACGACGGCCAgtgt                                                                                                                                     >  2:1091135/1‑150 (MQ=255)
                 catcatCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGACAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTc                                                                                                                                  >  1:1676493/1‑150 (MQ=255)
                  atcatcAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCg                                                                                                                                 >  2:1521923/1‑150 (MQ=255)
                  atcatcAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCg                                                                                                                                 <  1:461269/150‑1 (MQ=255)
                       caggcccaggccGATGGCCATGTCGCTGTAGCCGCAAACGTGGCGTTAGAGCGGTTCGGTGTAGGCCAGGGCGCAGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGAt                                                                                                                            <  1:1739451/150‑1 (MQ=255)
                           cccaggccGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAac                                                                                                                        <  1:1849680/150‑1 (MQ=255)
                                    aTGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTtg                                                                                                               <  2:1215290/150‑1 (MQ=255)
                                     tGGCCATGTCGCTGTAGCCGTTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTtga                                                                                                              <  1:1091135/150‑1 (MQ=255)
                                     tGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTtga                                                                                                              >  1:2129917/1‑150 (MQ=255)
                                       gCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTtgatg                                                                                                            >  1:1444316/1‑150 (MQ=255)
                                             tCGCTGTAGCCGAGGAAGACGGGGTACAGCGGTGCGGCGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGccc                                                                                                      >  2:215839/1‑150 (MQ=255)
                                                     gCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGaa                                                                                                          >  2:399576/1‑138 (MQ=255)
                                                     gCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGaa                                                                                                          <  1:399576/138‑1 (MQ=255)
                                                         cTGAAGTCGAAGTACAGCTGTGCGGTGTAGGACAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCgg                                                                                          <  2:1676493/150‑1 (MQ=255)
                                                          tGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCCGGGTGTCGGCGATGAACGCTTTCTTGATGAAGCCCTGCATGAAGCggg                                                                                         >  1:872282/1‑150 (MQ=255)
                                                          tGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCggg                                                                                         <  1:1521923/150‑1 (MQ=255)
                                                          tGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCggg                                                                                         >  1:247616/1‑150 (MQ=255)
                                                              gTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTc                                                                                                                                  >  1:672105/1‑105 (MQ=255)
                                                              gTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTc                                                                                                                                  <  2:672105/105‑1 (MQ=255)
                                                                     tACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCg                                                                              <  2:936718/150‑1 (MQ=255)
                                                                       cAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGga                                                                            <  1:1595839/150‑1 (MQ=255)
                                                                           tgtgCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAAc                                                                        <  1:81490/150‑1 (MQ=255)
                                                                             tgCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACtt                                                                      >  2:1311394/1‑150 (MQ=255)
                                                                                ggTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACttt                                                                                                                   >  1:239293/1‑102 (MQ=255)
                                                                                ggTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACttt                                                                                                                   <  2:239293/102‑1 (MQ=255)
                                                                                ggTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTc                                                                   >  1:2164597/1‑150 (MQ=255)
                                                                                  tgtAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCa                                                                 <  2:1444316/150‑1 (MQ=255)
                                                                                    tAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAgg                                                               >  1:2171961/1‑150 (MQ=255)
                                                                                      ggCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGgt                                                             <  2:247616/150‑1 (MQ=255)
                                                                                       gCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAaca                                                                                                                       >  1:1385533/1‑91 (MQ=255)
                                                                                       gCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAaca                                                                                                                       <  2:1385533/91‑1 (MQ=255)
                                                                                       gCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGgtg                                                            >  2:767005/1‑150 (MQ=255)
                                                                                        ccAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGgtgt                                                           >  1:2187637/1‑150 (MQ=255)
                                                                                               gcCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCg                                                    <  2:1885018/150‑1 (MQ=255)
                                                                                                ccGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACTGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCgg                                                   <  1:1235360/150‑1 (MQ=255)
                                                                                                       aGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAgg                                                                                                                                       <  2:754139/59‑1 (MQ=255)
                                                                                                       aGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAgg                                                                                                                                       >  1:754139/1‑59 (MQ=255)
                                                                                                        ggCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGaa                                           <  2:2164597/150‑1 (MQ=255)
                                                                                                              gCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGga                                                                            >  1:1703209/1‑111 (MQ=255)
                                                                                                              gCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGga                                                                            <  2:1703209/111‑1 (MQ=255)
                                                                                                              gCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTc                                     <  1:1311394/150‑1 (MQ=255)
                                                                                                                ccGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGgt                                                             <  1:1924651/124‑1 (MQ=255)
                                                                                                                ccGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGgt                                                             >  2:1924651/1‑124 (MQ=255)
                                                                                                                ccGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGa                                   <  2:1419818/150‑1 (MQ=255)
                                                                                                                      tggGGTTTTGCGGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCCGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGt                                      >  1:1672100/1‑141 (MQ=255)
                                                                                                                      tggGGTTTTGCGGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCCGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGt                                      <  2:1672100/141‑1 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTc                                                                               >  2:1733422/1‑100 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTc                                                                               >  2:1735073/1‑100 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTc                                                                               <  1:1733422/100‑1 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTc                                                                               <  1:1735073/100‑1 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAAt                             >  2:658577/1‑150 (MQ=255)
                                                                                                                      tggGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAAt                             >  1:1837552/1‑150 (MQ=255)
                                                                                                                       ggggTTTTGCAGGACGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGACGACCAAATc                            >  1:974258/1‑150 (MQ=255)
                                                                                                                        gggTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATcc                           >  2:191849/1‑150 (MQ=255)
                                                                                                                          gTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGAc                                  <  2:1681416/141‑1 (MQ=255)
                                                                                                                          gTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGAc                                  >  1:1681416/1‑141 (MQ=255)
                                                                                                                           ttttGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTg                        <  2:2129917/150‑1 (MQ=255)
                                                                                                                               gCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAgc                    >  2:1704217/1‑150 (MQ=255)
                                                                                                                                 aGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTc                                                                   >  1:844203/1‑101 (MQ=255)
                                                                                                                                 aGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTc                                                                   <  2:844203/101‑1 (MQ=255)
                                                                                                                                 aGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAgcgc                  <  1:215839/150‑1 (MQ=255)
                                                                                                                                     cGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGAGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAGCgcagca              <  2:1707621/150‑1 (MQ=255)
                                                                                                                                         gccgTGGTCGCCGACCACGGCCAGTGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAggg                                                              <  1:540576/95‑1 (MQ=255)
                                                                                                                                         gccgTGGTCGCCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAggg                                                              >  2:540576/4‑98 (MQ=255)
                                                                                                                                          ccgTGGTCGCCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGaa                                                                         <  2:332810/84‑1 (MQ=255)
                                                                                                                                          ccgTGGTCGCCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGaa                                                                         >  1:332810/3‑86 (MQ=255)
                                                                                                                                                tCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAGCGCAGCACGGGGCCTGCg   >  2:2250184/1‑150 (MQ=255)
                                                                                                                                                tCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGGGGGTGCGGTTGATGAACTGGTCGACCAAAATCATTAAGCTCTACACGGGGccacgt   >  2:521297/1‑146 (MQ=255)
                                                                                                                                                 cGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAGCGCAGCACGGGGACTGCGa  <  1:191849/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCCATGAAGCGGAAGCCCATCATCAGGCCCAGGCCGATGGCCATGTCGCTGTAGCCGCTGAAGTCGAAGTACAGCTGTGCGGTGTAGGCCAGGGCGCCGAGCCAGGCATCGCCCGTGGTGGGGTTTTGCAGGGCGAAGCAGTGGTCGGCGACCACGGCCAGGGTGTCGGCGATGAACACTTTCTTGATGAAGCCCTGCATGAAGCGGGTGCAGCCTTCGGAGAACTTGTCCAGGGTGTGGGTGCGGTTGTTGAACTGGTCGACCAAATCCTTGAAGCGCAGCACGGGGCCTGCGA  >  NC_002947/1462144‑1462438

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