Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,286 A→T 5.4% L42H (CTC→CAC)  TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,2860AT 5.4% 418.4 / 8.3 138L42H (CTC→CAC) TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (71/59);  new base T (4/4);  total (75/63)
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.99e-01

Only 100 of 140 total aligned reads displayed.

GGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297138‑5297433
                                                                                                                                                    |                                                                                                                                                   
ggCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAg                                                                                                                                                    <  1:2167380/150‑1 (MQ=255)
ggCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGCCCTCTAGTTCCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAg                                                                                                                                                    <  1:1354702/150‑1 (MQ=255)
  ccTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                                                                                  >  1:313738/1‑150 (MQ=255)
   cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                 >  2:2152396/1‑150 (MQ=255)
   cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                 >  1:2236734/1‑150 (MQ=255)
   cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                 >  1:951780/1‑150 (MQ=255)
    tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGTGctcg                                                                                                                                                >  1:1626688/1‑147 (MQ=255)
        ccACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                                                                                                                                            >  2:1697115/1‑150 (MQ=255)
          aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGTGCTCGCCctgc                                                                                                                                          <  2:1569426/150‑1 (MQ=255)
          aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGTGCTCGCCctgc                                                                                                                                          <  2:1626688/150‑1 (MQ=255)
           cttcttCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagca                                                                                                                                         <  1:2251458/150‑1 (MQ=255)
            ttcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCTGGGCGTGc                                                                                                                                                   <  2:706269/139‑1 (MQ=255)
             tcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                                                                       >  2:1229564/1‑150 (MQ=255)
               ttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                                                                     <  2:961432/134‑1 (MQ=255)
               ttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGt                                                                                                                                     >  1:1717364/1‑150 (MQ=255)
               ttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCGGTCATCGAACACCACAGCGGCTTCTACGGTGTAGTTGATCGAGCGGGCGATAAACTGCAGGTTATGGTTGGCATAGACAAGCAGGGCGAGCCCGCCCAGAAgggg                                                                                                                                     >  2:1136282/1‑149 (MQ=255)
                  ggTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCGGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                  <  1:1708762/150‑1 (MQ=255)
                   gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 <  1:1229564/150‑1 (MQ=255)
                   gTACTGGCAATCAACGCAAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 <  2:1159975/150‑1 (MQ=255)
                     aCTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                                                                               >  2:357795/1‑150 (MQ=255)
                     aCTGGAAATCAACGCCAGCGATTCGTGGGCGGTGGCTCTATGATCGAGCACCACAGCCGCTTTTACGGTGTAGTTGATCGAGCTGGCTATAAACTGCAGGTTATGGTTGGCATAGACACGCAGGGGGAGGACGCCCAGCAGGGTGAGGga                                                                                                                               <  2:791549/150‑1 (MQ=255)
                       tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgaga                                                                                                                             >  1:903892/1‑150 (MQ=255)
                       tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCCCTGCACCCGATAAGGCTTGCCGCTGTCAAGGTGTAGTTTGTCGAGCGGGCGATCAACTGCAGGATATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgaga                                                                                                                             >  1:1252990/1‑150 (MQ=255)
                        ggCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGAAACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGCGAt                                                                                                                            >  2:1036964/1‑150 (MQ=255)
                         gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGACATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATa                                                                                                                           >  1:989452/1‑150 (MQ=255)
                           aaTCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAcc                                                                                                                         >  1:2141131/1‑150 (MQ=255)
                            aTCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                                                                                                                        >  1:2039577/1‑150 (MQ=255)
                                     aGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                 >  1:201887/1‑116 (MQ=255)
                                     aGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                 <  2:201887/116‑1 (MQ=255)
                                             gTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaag                                                                                                       >  2:979441/1‑150 (MQ=255)
                                                  cggcggcgCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCACCAGGCCAACGGCAAGCAAgg                                                                                                  <  2:665306/150‑1 (MQ=255)
                                                         gcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                                                                            <  1:1604451/133‑1 (MQ=255)
                                                         gcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                                                                            >  2:1604451/1‑133 (MQ=255)
                                                         gcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTCGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCt                                                                                           >  1:1532977/1‑150 (MQ=255)
                                                            gTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  1:357795/150‑1 (MQ=255)
                                                            gTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  1:1649605/150‑1 (MQ=255)
                                                            gTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  2:1252990/150‑1 (MQ=255)
                                                               aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCTAGGCCAACGGCAAGCAAGGCGACGCTGAGGTg                                                                                     >  1:27275/1‑150 (MQ=255)
                                                               aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        >  2:1319198/1‑147 (MQ=255)
                                                               aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  1:1319198/147‑1 (MQ=255)
                                                               aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTg                                                                                     <  2:1917870/150‑1 (MQ=255)
                                                                tCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAc                                                                                                                          >  1:510560/1‑112 (MQ=255)
                                                                tCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                                                                  <  2:2109015/120‑1 (MQ=255)
                                                                tCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                                                                  >  1:2109015/1‑120 (MQ=255)
                                                                  gAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAg                                                                                                                                                    <  2:1810822/84‑1 (MQ=255)
                                                                  gAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAg                                                                                                                                                    >  1:1810822/1‑84 (MQ=255)
                                                                    acacCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAcc                                                                                                                         <  1:222434/109‑1 (MQ=255)
                                                                    acacCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAcc                                                                                                                         >  2:222434/1‑109 (MQ=255)
                                                                     caccacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACg                                                                               >  2:432885/1‑150 (MQ=255)
                                                                            gCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGc                                                                        >  1:1195705/1‑150 (MQ=255)
                                                                             ccGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCa                                                                       >  2:2089643/1‑150 (MQ=255)
                                                                                   ttACGGTGTAGTTGATCGAGCGGGCGATCAACTGCATGTTATGGTTGGCATACACACGCAGTCTGAGCACGCCCAGAAGGGTGAGTGAGATCCCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAgc                                                                 <  2:963426/149‑1 (MQ=255)
                                                                                       ggTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                  <  2:314006/81‑1 (MQ=255)
                                                                                        gtgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg                                                            >  1:1840438/1‑150 (MQ=255)
                                                                                           tAGTTGATCCAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTc                                                         <  1:1613362/150‑1 (MQ=255)
                                                                                            aGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACgg                                                                              >  2:644831/1‑128 (MQ=255)
                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagcag                                                                                                                                        >  2:1470691/1‑69 (MQ=255)
                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAgg                                                                                       >  1:632030/1‑118 (MQ=255)
                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        >  1:790060/1‑117 (MQ=255)
                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       >  1:203731/1‑150 (MQ=255)
                                                                                               tGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc                                                     <  2:903892/150‑1 (MQ=255)
                                                                                               tGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc                                                     <  2:1840438/150‑1 (MQ=255)
                                                                                                 aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgcgc                                                   >  1:825757/1‑150 (MQ=255)
                                                                                                    gAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       >  2:2316255/1‑143 (MQ=255)
                                                                                                    gAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       <  1:2316255/143‑1 (MQ=255)
                                                                                                     aGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACgg                                               >  2:694031/1‑150 (MQ=255)
                                                                                                     aGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACgg                                               <  1:694031/150‑1 (MQ=255)
                                                                                                           cGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGNGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        >  1:986981/1‑103 (MQ=255)
                                                                                                           cGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  2:986981/103‑1 (MQ=255)
                                                                                                            gATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCACCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCtt                                        >  2:1822804/1‑150 (MQ=255)
                                                                                                                  cTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGt                                  >  2:648820/1‑150 (MQ=255)
                                                                                                                   tGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTc                                 >  1:947942/1‑150 (MQ=255)
                                                                                                                      aGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCTAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGc                              <  2:27275/150‑1 (MQ=255)
                                                                                                                      aGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGc                              <  2:2039577/150‑1 (MQ=255)
                                                                                                                         ttATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       <  1:536250/122‑1 (MQ=255)
                                                                                                                         ttATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       >  2:536250/1‑122 (MQ=255)
                                                                                                                           aTGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCAt                         <  1:979441/150‑1 (MQ=255)
                                                                                                                            tGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCAt                         <  2:1357938/149‑1 (MQ=255)
                                                                                                                             ggTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATgg                       >  2:1538867/1‑150 (MQ=255)
                                                                                                                              gTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgc                      >  2:2053795/1‑150 (MQ=255)
                                                                                                                              gTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgc                      >  1:1967875/1‑150 (MQ=255)
                                                                                                                               ttGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAg                                                             >  1:1414009/1‑110 (MQ=255)
                                                                                                                               ttGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgct                     >  1:1414942/1‑150 (MQ=255)
                                                                                                                               ttGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAACGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgct                     >  2:923048/1‑150 (MQ=255)
                                                                                                                                   cgatctcacGCAGGGAGAGTACCCCCAGCAGGGTGAGGGAGATACTCACCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCATCACCGAGCTCAGGGTCGGGCGTACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTc                 <  1:587503/144‑1 (MQ=255)
                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTcc                <  2:1414942/150‑1 (MQ=255)
                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTcc                >  1:205792/1‑150 (MQ=255)
                                                                                                                                     tAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                                                                            >  2:1363157/1‑57 (MQ=255)
                                                                                                                                     tAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                                                                            <  1:1363157/57‑1 (MQ=255)
                                                                                                                                          ctgctgggCGTGCTCGCCCTGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGGCGGGCGCACGGTACGCt                                         >  1:261004/6‑119 (MQ=255)
                                                                                                                                          acGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCgg          >  1:358185/1‑150 (MQ=255)
                                                                                                                                            gctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCCACGGCAAGCAAGGCGACGCTGAGATGGCCAAGCCCCAGCACCGCGCGCAGGGTCGTGCGCACGGTACGCTTGCGCGTCGGCTTCATGGCTTCTCCGCCCGgcg        >  2:573539/4‑150 (MQ=255)
                                                                                                                                            gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAATGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgcg        >  1:1907509/1‑150 (MQ=255)
                                                                                                                                            gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgcg        >  1:998455/1‑150 (MQ=255)
                                                                                                                                              tgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGcgcc      >  2:833891/2‑150 (MQ=255)
                                                                                                                                               ggtcGAGCACCCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGAACGCAATTCGACGCGGTGGTGGCCACGGCCTAGCAGCGTGCGGAGGGTCAGGCGCACTGTACGCTTGCGGGGCGGCTTCATGGCTGCTCTGCCAGGccccc     >  2:408796/4‑149 (MQ=255)
                                                                                                                                               gggCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGcgccg     >  2:107206/1‑150 (MQ=255)
                                                                                                                                               gggCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGcgccg     <  2:947942/150‑1 (MQ=255)
                                                                                                                                                ggCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGcgccgc    <  2:203731/150‑1 (MQ=255)
                                                                                                                                                  cGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGa  <  2:989452/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297138‑5297433

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: