Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,821,669 C→T 5.2% F78F (TTC→TTT metW → methionine biosynthesis protein MetW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,821,6690CT 5.2% 226.8 / 2.2 76F78F (TTC→TTTmetWmethionine biosynthesis protein MetW
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (2/2);  total (37/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 = 7.20e-01

CGGTAGCGGCGAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCTGCCGCTGGTAC  >  NC_002947/5821522‑5821816
                                                                                                                                                    |                                                                                                                                                   
cGGTAGCGGCGAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGc                                                                                                                                                   >  2:422079/1‑150 (MQ=255)
 ggTAGCGGCGAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGcc                                                                                                                                                  >  2:76930/1‑150 (MQ=255)
 ggTAGCGGCGAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGcc                                                                                                                                                  <  2:758610/150‑1 (MQ=255)
  gtcgcggcgAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGCCAACGTCATCGAACAAGACCTGGACAAAGCCCTGGGCAACTTCGCCa                                                                                                                                                 <  2:889816/147‑1 (MQ=255)
      cggcgAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCaa                                                                                                                                             <  1:915005/150‑1 (MQ=255)
           aaCTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCt                                                                                                                                        >  1:481815/1‑150 (MQ=255)
             ctgctgGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTc                                                                                                                                      <  1:1719554/150‑1 (MQ=255)
                ctgGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGAc                                                                                                                                   <  1:570986/150‑1 (MQ=255)
                 tgGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTTGCCAGCAACAGCTTCGACg                                                                                                                                  >  2:197410/1‑150 (MQ=255)
                   gCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCg                                                                                                                                     >  2:1801825/1‑145 (MQ=255)
                   gCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCg                                                                                                                                     <  1:1801825/145‑1 (MQ=255)
                   gCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACgtg                                                                                                                                <  1:1821763/150‑1 (MQ=255)
                   gCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACgtg                                                                                                                                <  2:493545/150‑1 (MQ=255)
                     cTCACTGCCTGACCGCAAACAGCTCCCAGGCTATGGCCTGGAAATCGACGCCTACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACgtggg                                                                                                                              >  2:1506883/1‑149 (MQ=255)
                      tCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACgtggtg                                                                                                                             >  1:1718005/1‑150 (MQ=255)
                       cGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGa                                                                                                                            <  1:1012400/150‑1 (MQ=255)
                             gTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGa                                                                                                                      <  1:369176/150‑1 (MQ=255)
                             gTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGa                                                                                                                      >  2:32528/1‑150 (MQ=255)
                             gTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGa                                                                                                                      <  1:32528/150‑1 (MQ=255)
                              tGACCGCAAACAGGTCACCGGCTATGGCCTGGAACCCGACGCCGAAAATCACGCCGCCTTCTTGTTCTAGGGTGTCAAAAACCTCGAACCTGCGGTTGGCATAAGCCGGGTCAACTGCGCCAGCAACAGCTTCGACGTGGTGATCATGAc                                                                                                                     <  1:36139/150‑1 (MQ=255)
                                 ccGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCa                                                                                                                  >  2:114261/1‑150 (MQ=255)
                                  cGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCaa                                                                                                                 <  1:1362466/150‑1 (MQ=255)
                                        cAGGTCACCGGCTATGGCCTGGAAAACGCAGCCTACAATCTCGGTCGGATCGAAGGCTTGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTGCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTg                                                                                                           <  2:49920/150‑1 (MQ=255)
                                                cGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGt                                                                                                   <  2:481815/150‑1 (MQ=255)
                                                  gCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTgg                                                                                                 <  2:904163/150‑1 (MQ=255)
                                                        gCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTTGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATAcc                                                                                           <  1:197410/150‑1 (MQ=255)
                                                        gCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATAcc                                                                                           <  2:1176680/150‑1 (MQ=255)
                                                        gCCTGGAAATCGACGCCGACAATACAGTCCCTTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATAcc                                                                                           <  1:1526/150‑1 (MQ=255)
                                                         ccTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTTGCCAGCAACAGCTTCGACg                                                                                                                                  >  1:257965/1‑110 (MQ=255)
                                                         ccTGGAAATCGACGCCGACAAGATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTTGCCAGCAACAGCTTCGACg                                                                                                                                  <  2:257965/110‑1 (MQ=255)
                                                          cTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTTGGCAACTTCGCCAGCAACAGCGTCGACGTGGTGATCATGAGCCAAGCCCTGCAGGCCGTGGAATAcccc                                                                                         >  2:860910/1‑150 (MQ=255)
                                                           tGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCa                                                                                                                         <  2:184496/117‑1 (MQ=255)
                                                           tGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCa                                                                                                                         >  1:184496/1‑117 (MQ=255)
                                                             gAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGAc                                                                                      <  1:1611523/150‑1 (MQ=255)
                                                             gAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGAc                                                                                      <  1:114261/150‑1 (MQ=255)
                                                                 tCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCa                                                                                  >  1:123792/1‑150 (MQ=255)
                                                                 tCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCa                                                                                  >  1:1310396/1‑150 (MQ=255)
                                                                   gACGCCGACAATATCGCCGCCTGCGTTGCGAAGGGTGTCAACGTCATCGATCAAGACCTTTTCAAAGGCCTGTACCACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATc                                                                                <  1:698435/150‑1 (MQ=255)
                                                                   gACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATc                                                                                >  2:830639/1‑150 (MQ=255)
                                                                    aCGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATcc                                                                               <  1:76930/150‑1 (MQ=255)
                                                                    aCGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATcc                                                                               >  1:382309/1‑150 (MQ=255)
                                                                      gCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTc                                                                             >  1:196767/1‑150 (MQ=255)
                                                                       ccGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGaa                                                                                               <  2:650847/131‑1 (MQ=255)
                                                                       ccGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGaa                                                                                               >  1:650847/1‑131 (MQ=255)
                                                                        cGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGaa                                                                                               <  1:1062771/130‑1 (MQ=255)
                                                                        cGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGaa                                                                                               >  2:1062771/1‑130 (MQ=255)
                                                                         gACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTcgac                                                                          >  1:1848525/1‑150 (MQ=255)
                                                                         gACAATATCGCCGCCTGCGTGGCCAAGGGCGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTcgac                                                                          >  1:1365071/1‑150 (MQ=255)
                                                                            ttcaccgccgcCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGaa                                                                       >  2:1100468/6‑150 (MQ=255)
                                                                            aaTATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGaa                                                                       >  1:1887445/1‑150 (MQ=255)
                                                                             atatCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGACAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGaaa                                                                      <  2:1310396/150‑1 (MQ=255)
                                                                                   gccgccTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAAtgctg                                                                 <  1:422079/150‑1 (MQ=255)
                                                                                    ccgccTGCGTGGCCAAGGGTGGCAACGTCATCGAACAAGACCTGGAGAAAGGCCTGGGCAACTTCGCCAGCAACAGCTGCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAAtgctgc                                                                <  2:173026/150‑1 (MQ=255)
                                                                                                ccAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTc                                                    >  2:1745133/1‑150 (MQ=255)
                                                                                                ccAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTc                                                    >  1:900120/1‑150 (MQ=255)
                                                                                                 cAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCa                                                   >  1:1844973/1‑150 (MQ=255)
                                                                                                  aaGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGt                                                     <  2:1130664/147‑1 (MQ=255)
                                                                                                  aaGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGt                                                     >  1:1130664/1‑147 (MQ=255)
                                                                                                  aaGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAg                                                  >  1:1265312/1‑150 (MQ=255)
                                                                                                   aGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGt                                                 <  1:830639/150‑1 (MQ=255)
                                                                                                         tCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAAtgctg                                                                 >  2:1561893/1‑128 (MQ=255)
                                                                                                         tCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAAtgctg                                                                 <  1:1561893/128‑1 (MQ=255)
                                                                                                              gTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCAcc                                      <  2:123792/150‑1 (MQ=255)
                                                                                                                cATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCtt                                    >  1:685561/1‑150 (MQ=255)
                                                                                                                  tCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTcc                                  >  1:1712438/1‑150 (MQ=255)
                                                                                                                     aaCAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCatcatc                                         >  1:1139682/1‑140 (MQ=255)
                                                                                                                     aaCAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCatcatc                                         <  2:1139682/140‑1 (MQ=255)
                                                                                                                          gACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTc                          >  2:507333/1‑150 (MQ=255)
                                                                                                                           aCCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCg                         <  2:475120/150‑1 (MQ=255)
                                                                                                                            ccTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTg                                                           >  1:1568369/1‑115 (MQ=255)
                                                                                                                            ccTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTg                                                           <  2:1568369/115‑1 (MQ=255)
                                                                                                                             cTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGc                       <  1:1745133/150‑1 (MQ=255)
                                                                                                                                       gCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCt             >  1:655510/1‑150 (MQ=255)
                                                                                                                                               aaCTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGaaa                                                                      >  1:328359/1‑85 (MQ=255)
                                                                                                                                               aaCTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGaaa                                                                      <  2:328359/85‑1 (MQ=255)
                                                                                                                                                  ttCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCCGCCGCTGGTAc  >  1:741131/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGTAGCGGCGAACTGCTGGCCTCGCTGCGTGACCGCAAACAGGTCACCGGCTATGGCCTGGAAATCGACGCCGACAATATCGCCGCCTGCGTGGCCAAGGGTGTCAACGTCATCGAACAAGACCTGGACAAAGGCCTGGGCAACTTCGCCAGCAACAGCTTCGACGTGGTGATCATGACCCAAGCCCTGCAGGCCGTGGAATACCCCGACCGCATCCTCGACGAAATGCTGCGAGTGGGCCGTCAGTGCATCATCACCTTCCCGAACTTCGGCCACTGGCGCTGCCGCTGGTAC  >  NC_002947/5821522‑5821816

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: