Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,964,757 A→T 5.6% Q309L (CAG→CTG)  yegE → phosphodiesterase domain‑containing diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,964,7570AT 5.6% 220.3 / 4.4 72Q309L (CAG→CTG) yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base T (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

GGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGTGT  >  NC_002947/1964623‑1964906
                                                                                                                                      |                                                                                                                                                      
ggAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATTCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTc                                                                                                                                        <  2:1692735/150‑1 (MQ=255)
    gaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCa                                                                                                                                    <  2:648079/150‑1 (MQ=255)
    gaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCAGGCCCTGGGCCTggccg                                                                                                                                            >  2:1464283/1‑137 (MQ=255)
    gaaAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCAGGCCCTGGGCCTggccg                                                                                                                                            <  1:1464283/141‑5 (MQ=255)
                 cgcAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCtgg                                                                                                                       >  1:1040289/1‑150 (MQ=255)
                   cAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGTTTGATTCTGGAGGAGCTGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCtggtg                                                                                                                     <  2:50142/150‑1 (MQ=255)
                     gATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGa                                                                                                                   <  1:1482623/150‑1 (MQ=255)
                      aTTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGc                                                                                                                                               <  2:206980/121‑1 (MQ=255)
                      aTTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGc                                                                                                                                               >  1:206980/1‑121 (MQ=255)
                        tACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCACAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGa                                                                                                                <  1:945188/150‑1 (MQ=255)
                         aCCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAg                                                                                                               <  2:1040289/150‑1 (MQ=255)
                          cccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGc                                                                                                                           <  1:164761/137‑1 (MQ=255)
                          cccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGc                                                                                                                           >  2:164761/1‑137 (MQ=255)
                          cccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAt                                                                                                                  <  2:6331/146‑1 (MQ=255)
                          cccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGAt                                                                                                                  >  1:6331/1‑146 (MQ=255)
                          cccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGc                                                                                                              <  2:1117603/150‑1 (MQ=255)
                           ccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGc                                                                                                                                                    <  2:921550/112‑1 (MQ=255)
                           ccTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGc                                                                                                                                                    >  1:921550/1‑112 (MQ=255)
                             tGCAGGCCATCGGGGCTGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCCACCAGGCCCAGGGCCTGCCGCTATCGGCACTGTTCACCCTGGTGGCTGAGCACg                                                                                                           >  2:457606/1‑150 (MQ=255)
                             tGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACg                                                                                                           >  1:567370/1‑150 (MQ=255)
                                    cATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCgagga                                                                                                    <  2:567370/150‑1 (MQ=255)
                                      tCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGaggagg                                                                                                  >  2:875212/1‑150 (MQ=255)
                                        ggggATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGGGCAGTTGATCAACTGGCATTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGAt                                                                                                <  2:1426132/150‑1 (MQ=255)
                                           gATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGt                                                                                             <  1:399719/150‑1 (MQ=255)
                                           gATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGt                                                                                             <  1:1889560/150‑1 (MQ=255)
                                                ggTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAg                                                                                        >  2:729361/1‑150 (MQ=255)
                                                    aTTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCCTGCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTg                                                                                    <  1:1090031/150‑1 (MQ=255)
                                                       aCAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                                 <  1:323311/150‑1 (MQ=255)
                                                        cAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCg                                                                                >  1:555029/1‑150 (MQ=255)
                                                         aaCCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGa                                                                               >  1:64050/1‑150 (MQ=255)
                                                          aCCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAg                                                                              <  2:1335669/150‑1 (MQ=255)
                                                          aCCGATGTCGAGCGCAATATCGGCTACATGAACCCGGCCGCCGATCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTATCGGCAGTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAg                                                                              <  1:457606/150‑1 (MQ=255)
                                                              ctgacGAGCGCAATATCGGCTCCATCAACCCGCCCGCCGAGCAATTCACCCACTGCCAATCCGGCCAGGCCCATGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAgg                                                                          <  2:693816/146‑1 (MQ=255)
                                                                gTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTg                                                                        <  1:875212/150‑1 (MQ=255)
                                                                        cAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCgg                                                                >  1:1773953/1‑150 (MQ=255)
                                                                            atCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGc                                                            <  2:64050/150‑1 (MQ=255)
                                                                             tCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCa                                                                    >  2:495099/1‑141 (MQ=255)
                                                                             tCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCa                                                                    <  1:495099/141‑1 (MQ=255)
                                                                               ggCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCGAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGAAGGTCGCAGATCTGTCGATCAGGTGCTCAGCGGCAGCCTc                                                         >  1:195/1‑150 (MQ=255)
                                                                                  tACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAg                                                      >  2:811152/1‑150 (MQ=255)
                                                                                     aTGAACCCGGCCGCCGAGCAGTTGACCCCCTGGCACTACGGCGGTGCCCAGGGCCTGCCGCTGTCGGCAATGTTTAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGggc                                                   <  1:434140/150‑1 (MQ=255)
                                                                                     aTGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGggc                                                   >  2:1152138/1‑150 (MQ=255)
                                                                                             ggCCGCCGAGCAGTTGACCCGCTTGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCGGCACTGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGa                                           <  1:612741/150‑1 (MQ=255)
                                                                                             ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTg                                                                                    >  2:134970/1‑109 (MQ=255)
                                                                                             ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTg                                                                                    <  1:134970/109‑1 (MQ=255)
                                                                                             ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGggc                                                   <  1:314277/142‑1 (MQ=255)
                                                                                             ggCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGggc                                                   >  2:314277/1‑142 (MQ=255)
                                                                                                    gAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                       >  1:349760/1‑131 (MQ=255)
                                                                                                    gAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCaa                                                       <  2:349760/131‑1 (MQ=255)
                                                                                                         gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                                 >  2:65180/1‑100 (MQ=255)
                                                                                                         gTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTc                                                                                 <  1:65180/100‑1 (MQ=255)
                                                                                                          ttGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTg                              >  1:1734887/1‑150 (MQ=255)
                                                                                                              cccACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATcc                          >  1:1379255/1‑150 (MQ=255)
                                                                                                                cACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCt                                               >  2:478342/1‑127 (MQ=255)
                                                                                                                cACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCt                                               <  1:478342/127‑1 (MQ=255)
                                                                                                                 aCTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAgc                       >  1:302794/1‑150 (MQ=255)
                                                                                                                   tGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAgcgc                     >  2:1040396/1‑150 (MQ=255)
                                                                                                                    ggCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGTATGCCCGGCTGATCCAGCGcc                    >  2:564099/1‑150 (MQ=255)
                                                                                                                     gCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCt                   <  2:1734887/150‑1 (MQ=255)
                                                                                                                      cAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTg                  <  2:302794/150‑1 (MQ=255)
                                                                                                                      cAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTg                  <  1:811152/150‑1 (MQ=255)
                                                                                                                       aaTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGc                                                               >  1:748013/1‑104 (MQ=255)
                                                                                                                       aaTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGc                                                               <  2:748013/104‑1 (MQ=255)
                                                                                                                        aTCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGa                <  2:1773953/150‑1 (MQ=255)
                                                                                                                         tCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATg                                                                                               >  2:630340/1‑70 (MQ=255)
                                                                                                                         tCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATg                                                                                               <  1:630340/70‑1 (MQ=255)
                                                                                                                            ggCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGc            >  2:1625168/1‑150 (MQ=255)
                                                                                                                            cggcAGGCCCTGGGCCTGGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGa                                           >  1:1067878/4‑120 (MQ=255)
                                                                                                                            cggcAGGCCCTGGGCCTGGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGa                                           <  2:1067878/117‑1 (MQ=255)
                                                                                                                             gccagccaCAGCGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGACGAGGAGGATGGTCGCAGCTTTGTCGAGCAGGTGCTCAGCGGCAGCCTTAAGGGCTGCTCCGAGAATGCCCGGCTGATCCAGCGCCTGGACGGCa           <  2:195/150‑1 (MQ=255)
                                                                                                                                  gCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTAcc      >  1:1319887/1‑150 (MQ=255)
                                                                                                                                    ccAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCgt    >  1:583829/1‑150 (MQ=255)
                                                                                                                                      aGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGTTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGTTCCAGCGCCTGGACGGCAGTACCgtgt  <  2:135498/150‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                      
GGAAGAAAAAGAGCGTGCGCAGATTACCCTGCAGGCCATCGGGGATGCGGTGATTACAACCGATGTCGAGGGCAATATCGGCTACATGAACCCGGCCGCCGAGCAGTTGACCCACTGGCAATCCGGCCAGGCCCAGGGCCTGCCGCTGTCGGCACTGTTCAGCCTGGTGGATGAGCACGCCGAGGAGGATGGTCGCAGCTTGGTCGAGCAGGTGCTCAGCGGCAGCCTCAAGGGCGGCTCCGAGCATGCCCGGCTGATCCAGCGCCTGGACGGCAGTACCGTGT  >  NC_002947/1964623‑1964906

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