Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,996,057 C→G 9.8% intergenic (‑260/+25) PP_mr26 ← / ← rfbC C4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,996,0570CG 9.8% 239.1 / 16.6 82intergenic (‑260/+25)PP_mr26/rfbCC4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base G (5/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

TCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAATGCAGCGCTC  >  NC_002947/1995916‑1996204
                                                                                                                                             |                                                                                                                                                   
tCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGcca                                                                                                                                             <  2:55655/150‑2 (MQ=255)
  tttGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACCAGGCTGCGCTACCGCAGCCCtt                                                                                                                                           >  1:245805/1‑150 (MQ=255)
  tttGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCtt                                                                                                                                           >  1:1148319/1‑150 (MQ=255)
  tttGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGACACAAGGCTGCGCTACCGCAGCCCtt                                                                                                                                           <  1:1098654/150‑1 (MQ=255)
      ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTcc                                                                                                                                       <  1:1088294/150‑1 (MQ=255)
              ccTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAAc                                                                                                                               >  1:716932/1‑150 (MQ=255)
                tACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCg                                                                                                                             >  1:1911693/1‑150 (MQ=255)
                  cAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGTAGCCCTTCTCCGCTAAAACCGtt                                                                                                                           >  2:336706/1‑150 (MQ=255)
                  cAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGtt                                                                                                                           >  1:264415/1‑150 (MQ=255)
                     ggaaggGGCAAATTCTGATTATTGGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGt                                                                                                                        <  1:574565/150‑1 (MQ=255)
                         ggggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGaa                                                                                                                    <  1:651995/150‑1 (MQ=255)
                          gggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCTCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAAc                                                                                                                   >  1:100052/1‑150 (MQ=255)
                          gggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAAc                                                                                                                   >  2:723002/1‑150 (MQ=255)
                          gggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAAc                                                                                                                   >  2:758890/1‑150 (MQ=255)
                            gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAg                                                                                                                 >  1:183339/1‑150 (MQ=255)
                               aaTTCTGATTATTTGCCAGGCCATGCAGCGGTTGCCTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTc                                                                                                              <  2:1301354/150‑1 (MQ=255)
                               aaTTCTGATTATTTGCCAGGCCATGCAGCGGTTGCCTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTc                                                                                                              <  2:1301649/150‑1 (MQ=255)
                                aTTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTTTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTGGTCGAACAGTTCg                                                                                                             >  2:492217/1‑150 (MQ=255)
                                aTTCTGATTATTCGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGGTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCg                                                                                                             <  1:686715/150‑1 (MQ=255)
                                 ttCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCgg                                                                                                            >  1:545940/1‑150 (MQ=255)
                                  tCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGc                                                                                                           <  2:1924218/150‑1 (MQ=255)
                                    tGATTATTTGCCAGCCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTGAGTCTACCAGTCCGGAGt                                                                                                         >  2:339857/1‑150 (MQ=255)
                                     gATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTc                                                                                                        >  2:1503566/1‑150 (MQ=255)
                                      attattTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCa                                                                                                       <  2:1487068/150‑1 (MQ=255)
                                          tttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATTGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCa                                                                                                   >  2:516201/1‑150 (MQ=255)
                                           ttGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCaa                                                                                                  >  1:641578/1‑150 (MQ=255)
                                            tGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAAt                                                                                                 <  1:1909783/150‑1 (MQ=255)
                                             gCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATg                                                                                                >  1:622643/1‑150 (MQ=255)
                                                        cAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTCTCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGc                                                                                     <  1:1259777/150‑1 (MQ=255)
                                                         aGCGGTTGCTTCTTAGTGGCGTGTTAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCt                                                                                    >  1:376219/1‑150 (MQ=255)
                                                          gCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTg                                                                                   <  2:1497753/150‑1 (MQ=255)
                                                          gCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTg                                                                                   <  1:748455/150‑1 (MQ=255)
                                                                 cttcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTg                                                                            <  2:545940/150‑1 (MQ=255)
                                                                 cttcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                                                                           <  1:294398/150‑1 (MQ=255)
                                                                    cttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCa                                                                         <  2:1745952/150‑1 (MQ=255)
                                                                     ttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTACACCACATTGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAg                                                                        >  2:259467/1‑150 (MQ=255)
                                                                             gtgtGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGTAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCtt                                                                <  1:336706/150‑1 (MQ=255)
                                                                             gtgtGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCtt                                                                >  2:641574/1‑150 (MQ=255)
                                                                              tgtgAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCACGCCGGTCCTTGCCAGACAAGCTTg                                                               <  2:1911693/150‑1 (MQ=255)
                                                                                 gAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTa                                                            <  1:641574/150‑1 (MQ=255)
                                                                                 aaaaaaCGCCACTATCGGCAAGCCCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCGGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTa                                                            <  2:1688228/149‑1 (MQ=255)
                                                                                  aaaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc                                                           <  2:1631536/150‑1 (MQ=255)
                                                                                     aaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGa                                                        <  1:827789/150‑1 (MQ=255)
                                                                                     aaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGa                                                        >  2:743267/1‑150 (MQ=255)
                                                                                      aCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGATACCCCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAg                                                       >  2:1606755/1‑150 (MQ=255)
                                                                                        gCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTc                                                     <  1:723002/150‑1 (MQ=255)
                                                                                          cACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAc                                                                                            <  2:240036/109‑1 (MQ=255)
                                                                                          cACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAc                                                                                            >  1:240036/1‑109 (MQ=255)
                                                                                           aCTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGAt                                                  <  2:716932/150‑1 (MQ=255)
                                                                                           aCTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGAt                                                  <  1:758890/150‑1 (MQ=255)
                                                                                             tATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTGGGTACGGAGTCGATgg                                                >  2:972296/1‑150 (MQ=255)
                                                                                               ttggCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGc                                              <  1:516201/148‑1 (MQ=255)
                                                                                               tCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGc                                              <  2:135323/150‑1 (MQ=255)
                                                                                                  gCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGc                                                                  <  2:632032/127‑1 (MQ=255)
                                                                                                  gCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGc                                                                  >  1:632032/1‑127 (MQ=255)
                                                                                                  gCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCAt                                           >  1:737649/1‑150 (MQ=255)
                                                                                                    aaGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTc                                         <  2:62263/150‑1 (MQ=255)
                                                                                                      gACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCGACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATTCGCCCCCAAGCTGGTCCTTGACAGACACGCTTGTTACGGAGTCGATGGGCCATTCGa                                       <  1:42687/150‑1 (MQ=255)
                                                                                                         ccGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATAc                                    <  1:693697/150‑1 (MQ=255)
                                                                                                          cGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                               >  2:1524549/1‑122 (MQ=255)
                                                                                                          cGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTg                                                               <  1:1524549/122‑1 (MQ=255)
                                                                                                               tGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCAtt                                          <  1:81391/138‑1 (MQ=255)
                                                                                                               tGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCAtt                                          >  2:81391/1‑138 (MQ=255)
                                                                                                                  cATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCg                           >  2:893545/1‑150 (MQ=255)
                                                                                                                   aTGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCgg                          >  2:1395089/1‑150 (MQ=255)
                                                                                                                   aTGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCgg                          >  2:923383/1‑150 (MQ=255)
                                                                                                                      tgtTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATc                       <  1:923383/150‑1 (MQ=255)
                                                                                                                       gtTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCa                      <  2:183339/150‑1 (MQ=255)
                                                                                                                           agagACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTcc                  >  1:1576871/1‑150 (MQ=255)
                                                                                                                            ggagaagggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCa                 >  1:1690898/8‑150 (MQ=255)
                                                                                                                            gagaCAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCa                 <  2:100052/150‑1 (MQ=255)
                                                                                                                             gagaagggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAg                <  1:931534/144‑1 (MQ=255)
                                                                                                                             gagaagggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAg                >  2:931534/7‑150 (MQ=255)
                                                                                                                             agaCAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAg                <  2:641578/150‑1 (MQ=255)
                                                                                                                              gaCAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAgg               >  1:1808056/1‑150 (MQ=255)
                                                                                                                               aCAAGGCCGCGTTACTGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAACCAAAGCGACCCCCAGCTGGACCTTGCCAGACCAGCTTCCTACGGAGTCGATGGCCTATTCGATACCCACTGCCGGATCATTCCAGGc              <  1:1077283/150‑1 (MQ=255)
                                                                                                                                aagggCTGCGGTAGCGCAGCCTTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCa             >  2:1397232/4‑150 (MQ=255)
                                                                                                                                 agggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc                                                           >  1:1504939/3‑103 (MQ=255)
                                                                                                                                 agggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc                                                           <  2:1504939/101‑1 (MQ=255)
                                                                                                                                  gggCTGCGGTAGCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAAt           >  2:1471957/2‑150 (MQ=255)
                                                                                                                                    gCTGCGCTACCGCAGCCCATCTACGCTAAAAACGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCACAAGCTGGTCCTTGCCAGACAAGCTTTTTACGGCGACGATGGGCGAGTCGATAACCACTGTCTGTTCATTCCAGGGACTGc         >  2:562767/1‑150 (MQ=255)
                                                                                                                                           tACCGCAGCCCTTCTTCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTACGGTATGATTCCAGGCAATGCAGCGCTc  <  2:376219/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACGGAGTCGATGGGCCATTCGATACCCACTGCCGGATCATTCCAGGCAATGCAGCGCTC  >  NC_002947/1995916‑1996204

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: