Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,675,646 C→G 5.2% G116A (GGC→GCC)  yceJ ← iron‑regulated cytochrome b561

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,675,6460CG 5.2% 284.7 / 2.8 89G116A (GGC→GCC) yceJiron‑regulated cytochrome b561
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base G (3/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCTGATGAAACGC  >  NC_002947/5675497‑5675786
                                                                                                                                                     |                                                                                                                                             
gCTTGAGCGCGGCAAGGGCATGCAGCACGGCCAACATGATCAGGCCAACCTCACGCTGCAAGAGAACGGCACGTGCCACATCTGCCTGGTAAGGCAGGTCGCAGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGc                                                                                                                                               >  1:1252321/1‑150 (MQ=255)
  ttGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCg                                                                                                                                             >  2:2359551/1‑150 (MQ=255)
  ttGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCg                                                                                                                                             >  2:178720/1‑150 (MQ=255)
          ggCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGt                                                                                                                                     <  2:1366994/150‑1 (MQ=255)
               gggCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCTCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGaaa                                                                                                                                >  1:2318591/1‑150 (MQ=255)
                 gCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATc                                                                                                                              >  2:1708943/1‑150 (MQ=255)
                  cATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCa                                                                                                                             >  2:392094/1‑150 (MQ=255)
                   aTGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAg                                                                                                                            >  2:471854/1‑150 (MQ=255)
                          aaggaaaaacTCGCCCGCCCCCAAGACCGAGAGAAATCAACCGCACCCGCCCCAACTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACGTCGAACGGGCCGAACACCGGAATGCCGACGCCGGCGGCGGTGGAAATCAGGTAACCg                                                                                                                     <  2:1768837/145‑1 (MQ=255)
                                aaaaTCACCAGCCCCCAGGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAAcc                                                                                                               >  2:532702/1‑150 (MQ=255)
                                 aaaTCACCAGCCCCCAAGCCAGCCAAAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCa                                                                                                              >  1:1531512/1‑150 (MQ=255)
                                      aCCAGCCCCCAAGCCAGCCACAAATGACCCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAaccaccacc                                                                                                         <  2:763075/150‑1 (MQ=255)
                                       ccAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCg                                                                                                        >  2:2081777/1‑150 (MQ=255)
                                        cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGGCGTcggcat                                                                                                                                      >  1:1634628/1‑117 (MQ=255)
                                        cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGGCGTcggcat                                                                                                                                      <  2:1634628/119‑3 (MQ=255)
                                        cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGt                                                                                                                          >  1:1608925/1‑131 (MQ=255)
                                        cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGt                                                                                                                          <  2:1608925/131‑1 (MQ=255)
                                        cAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGc                                                                                                       >  2:1236454/1‑150 (MQ=255)
                                            ccccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAAc                                                                                                   <  2:2230816/150‑1 (MQ=255)
                                             cccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACa                                                                                                  >  1:1422235/1‑150 (MQ=255)
                                              ccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATc                                                                                                                              >  2:1049710/1‑121 (MQ=255)
                                              ccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATc                                                                                                                              <  1:1049710/121‑1 (MQ=255)
                                              ccAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAg                                                                                                 >  2:1558160/1‑150 (MQ=255)
                                                   ccagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGTACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCaa                                                                                           <  1:1877113/150‑1 (MQ=255)
                                                   ccagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCa                                                                                                              >  2:1482985/1‑132 (MQ=255)
                                                   ccagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCa                                                                                                              <  1:1482985/132‑1 (MQ=255)
                                                    cagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCTACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAg                                                                                          <  2:1561457/150‑1 (MQ=255)
                                                    cagccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAg                                                                                          <  2:907661/150‑1 (MQ=255)
                                                     agccaCAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAgg                                                                                         >  2:896688/1‑150 (MQ=255)
                                                        cacaAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTAc                                                                                      >  2:1593601/1‑150 (MQ=255)
                                                        cacaAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTAc                                                                                      >  1:85765/1‑150 (MQ=255)
                                                           aaaTGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAgg                                                                                   <  1:2081777/150‑1 (MQ=255)
                                                              tGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAggggc                                                                                >  1:411436/1‑148 (MQ=255)
                                                              tGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGcc                                                                                <  1:1877522/150‑1 (MQ=255)
                                                               gAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCa                                                                               <  1:471854/150‑1 (MQ=255)
                                                                aaCCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTCCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAg                                                                              <  1:790941/150‑1 (MQ=255)
                                                                  ccACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACCGCCACCATGACAGGGCCAGGt                                                                            >  2:1044532/1‑150 (MQ=255)
                                                                    acacCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAg                                                                              >  1:671154/1‑146 (MQ=255)
                                                                    acacCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAg                                                                              <  2:671154/146‑1 (MQ=255)
                                                                      acCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGcc                                                                        >  2:481686/1‑150 (MQ=255)
                                                                        ccGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAgg                                                                                                                           >  2:1756908/1‑98 (MQ=255)
                                                                        ccGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAgg                                                                                                                           <  1:1756908/98‑1 (MQ=255)
                                                                        ccGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGcc                                                                                >  1:830271/1‑140 (MQ=255)
                                                                        ccGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGcc                                                                                <  2:830271/140‑1 (MQ=255)
                                                                        ccGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCa                                                                      <  1:2359551/150‑1 (MQ=255)
                                                                           ccACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAAACGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCt                                                                   >  2:1818722/1‑150 (MQ=255)
                                                                               aTCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgc                                                               >  1:1828447/1‑150 (MQ=255)
                                                                               aTCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgc                                                               >  1:174623/1‑150 (MQ=255)
                                                                                  tGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGgccgc                                                            >  2:2124062/1‑150 (MQ=255)
                                                                                     cTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTg                                                                 >  2:607842/1‑142 (MQ=255)
                                                                                     cTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCAGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTg                                                                 <  1:607842/142‑1 (MQ=255)
                                                                                     cTGGTCAGGCAGGTCGCTGATCAGTCCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAg                                                         >  2:2357829/1‑150 (MQ=255)
                                                                                          caggcaggTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAAACGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCgtg                                                    <  1:1818722/150‑1 (MQ=255)
                                                                                          caggcaggTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGGCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCgtg                                                    <  2:178685/150‑1 (MQ=255)
                                                                                           aggcaggTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCgtgt                                                   <  1:1558160/150‑1 (MQ=255)
                                                                                             gcaggTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGa                                                 >  2:543457/1‑150 (MQ=255)
                                                                                              caggTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGAc                                                <  1:178720/150‑1 (MQ=255)
                                                                                                    gCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCAcc                                          >  1:793364/1‑150 (MQ=255)
                                                                                                            gTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAAcc                                                                                                                      <  2:2207181/67‑1 (MQ=255)
                                                                                                            gTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAAcc                                                                                                                      >  1:2207181/1‑67 (MQ=255)
                                                                                                            gTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttgg                                  >  2:2156707/1‑150 (MQ=255)
                                                                                                                      ttCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttggttggt                             <  2:831510/145‑1 (MQ=255)
                                                                                                                      ttCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttggttggt                             >  1:831510/1‑145 (MQ=255)
                                                                                                                      ttCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCgg                        <  2:1159639/150‑1 (MQ=255)
                                                                                                                         gAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGTTTGGTTGGTGCCggtgg                     >  1:1394274/1‑150 (MQ=255)
                                                                                                                          aaCAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCaa                                                                                                                 >  1:1323832/1‑58 (MQ=255)
                                                                                                                          aaCAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCaa                                                                                                                 <  2:1323832/58‑1 (MQ=255)
                                                                                                                           aCAGGCCGAACACCGGAATGCCGCCGCCGTCGGCGGGGGAAACCAGGTAACCGGCAACCATCACCACAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGTCGCCAGGCGCGTGACCGCACCATGGTTGGTTGGTGCCGgtggtg                   >  2:1736512/1‑150 (MQ=255)
                                                                                                                             aGGCCGAACACCGGCAGGTCGACGCCGTCGGCGGTTGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGCTGGCCCAGCTTGGCGGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGc                 <  2:793364/150‑1 (MQ=255)
                                                                                                                             aGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTATGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGc                 <  2:174623/150‑1 (MQ=255)
                                                                                                                              ggCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCg                >  1:86624/1‑150 (MQ=255)
                                                                                                                               gCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCgg               <  1:1236454/150‑1 (MQ=255)
                                                                                                                                ccGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCggg              <  1:1708943/150‑1 (MQ=255)
                                                                                                                                 cGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGc             <  2:1422235/150‑1 (MQ=255)
                                                                                                                                  gAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACATGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCt            >  1:900508/1‑150 (MQ=255)
                                                                                                                                  gAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttggttg                               >  1:1301207/1‑131 (MQ=255)
                                                                                                                                  gAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttggttg                               <  2:1301207/131‑1 (MQ=255)
                                                                                                                                    acacCGGAATGCCGACGGCGTCGGCGGTGGAAATCGGGTACCCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtga          <  2:1092873/150‑1 (MQ=255)
                                                                                                                                    acacCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtga          <  2:85765/150‑1 (MQ=255)
                                                                                                                                     cacCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtgat         <  1:2075090/150‑1 (MQ=255)
                                                                                                                                     cacCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtgat         <  1:543457/150‑1 (MQ=255)
                                                                                                                                      acCGGAATGCCGACGCCGTCGGTGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtgatg        >  1:2329876/1‑150 (MQ=255)
                                                                                                                                      acCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtgatg        >  2:1860676/1‑150 (MQ=255)
                                                                                                                                       ccGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCtgatga       <  2:900508/150‑1 (MQ=255)
                                                                                                                                         ggAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAGCAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGGCTGGGTCCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCTGATGaaa     <  1:1044532/150‑1 (MQ=255)
                                                                                                                                          naaTGCCGACGCCGTCGGCGGTGGAAATNAGGTAACCNGCAACCATCNCCGCAAACAGCCCAAGGTACAGGGGCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCTGATGaaag    >  1:996815/2‑149 (MQ=255)
                                                                                                                                            ccgCCGACGGCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCACAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCTGATGAAACGc  >  1:1290434/3‑150 (MQ=255)
                                                                                                                                                cGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttgg                                  <  1:879810/114‑1 (MQ=255)
                                                                                                                                                cGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCAtggttgg                                  >  2:879810/1‑114 (MQ=255)
                                                                                                                                                  aCGGCGTCGGCAGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGccc                                                                       >  2:491960/1‑75 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GCTTGAGCGCGGCAAGGGCATGCAGCACGGCAAAAATCACCAGCCCCCAAGCCAGCCACAAATGAACCACACCCGCCACATCTGCCTGGTCAGGCAGGTCGCTGATCAGTGCCGGCACTTCGAACAGGCCGAACACCGGAATGCCGACGCCGTCGGCGGTGGAAATCAGGTAACCGGCAACCATCACCGCAAACAGCCCAAGGTACAGGGCCAGGTGGCCCAGCTTGGCCGCCAGGCGTGTGACCGCACCATGGTTGGTTGGTGCCGGTGGTGGCGGGCTGATGAAACGC  >  NC_002947/5675497‑5675786

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: