Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,987 C→A 7.3% intergenic (+27/‑137) katA → / → bfr‑I catalase/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,9870CA 7.3% 234.3 / 7.9 82intergenic (+27/‑137)katA/bfr‑Icatalase/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base A (3/3);  total (40/42)
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.70e-01

CCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCA  >  NC_002947/565842‑566134
                                                                                                                                                 |                                                                                                                                                    
ccGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAt                                                                                                                                                    <  1:315608/148‑1 (MQ=255)
 cGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAttt                                                                                                                                                  <  2:138686/149‑1 (MQ=255)
      ccATGGCGGGCGTCAGCGAGGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAAACTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCACATTTGGGCg                                                                                                                                             >  2:92567/1‑149 (MQ=255)
           gcgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttc                                                                                                                                        <  5:130689/148‑2 (MQ=255)
           gcgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttc                                                                                                                                        <  7:249244/148‑2 (MQ=255)
            cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGttttt                                                                                                                                      <  3:126415/149‑1 (MQ=255)
            cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGttttt                                                                                                                                      <  4:112850/149‑1 (MQ=255)
             gggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAAAACCGCCCAAATGAGGGCGGttcttc                                                                                                                                     <  2:69718/149‑1 (MQ=255)
             gggcgTCAGCGAAGATGGCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGCGATTGCAAAGTCTTTGGGCATCATTCTCGCCTAAGTCGAAGTCATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTc                                                                                                                                     <  7:307124/149‑1 (MQ=255)
                cgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAAt                                                                                                                                  <  7:153074/149‑1 (MQ=255)
                cgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAAt                                                                                                                                  <  5:2900/149‑1 (MQ=255)
                    aGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGaa                                                                                                                               <  6:156622/148‑1 (MQ=255)
                     gCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAata                                                                                                                             >  5:33138/1‑149 (MQ=255)
                      cGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAatat                                                                                                                            <  4:217359/149‑1 (MQ=255)
                      cGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAata                                                                                                                             >  3:210043/1‑148 (MQ=255)
                          gATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAATGGCTTTGGGCATCAATCTCGCCTAAGTCTAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCac                                                                                                                         <  1:92567/148‑1 (MQ=255)
                                aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                  >  5:179433/1‑149 (MQ=255)
                                aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                  >  2:60838/1‑149 (MQ=255)
                                aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                  <  8:133548/149‑1 (MQ=255)
                                 ttCAGCGTCAGCTACAGTACTTCTTCAAGGCGGAACCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                  >  5:211377/1‑148 (MQ=255)
                                    aGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAc                                                                                                              <  8:205722/149‑1 (MQ=255)
                                           gCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAtttt                                                                                                       >  2:231944/1‑149 (MQ=255)
                                           gCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATAACTACTGTGTAACCGAtttt                                                                                                       >  8:335438/1‑149 (MQ=255)
                                               cAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGc                                                                                                   >  3:244354/1‑149 (MQ=255)
                                               cAGTACTTCTTCAAGGCGGAACCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAtttttt                                                                                                     <  6:211377/147‑1 (MQ=255)
                                                  tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCAGTACTGTTTAACGGATTTTTTGAGTc                                                                                                >  6:301977/1‑149 (MQ=255)
                                                  tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                >  2:359701/1‑149 (MQ=255)
                                                    cttcttCAAGGCTGATCCGTCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTTATAAGAAGAACCGCACTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                               <  7:335438/148‑1 (MQ=255)
                                                    cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCTATTTTTTGAGTCt                                                                                               >  8:83231/1‑148 (MQ=255)
                                                    cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                               >  6:193360/1‑148 (MQ=255)
                                                     ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                             >  1:264351/1‑149 (MQ=255)
                                                       cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                           >  7:267523/1‑149 (MQ=255)
                                                       cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                           >  2:23041/1‑149 (MQ=255)
                                                          cAAGGCGGATGCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATGAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                         <  5:206723/148‑1 (MQ=255)
                                                          cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAtt                                                                                                                                                   <  3:223835/91‑1 (MQ=255)
                                                          cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAtt                                                                                                                                                   >  4:223835/1‑91 (MQ=255)
                                                           aaGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgcgc                                                                                       >  3:155246/1‑149 (MQ=255)
                                                            aGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTCTTTAACCGATTTTTTGCGTCTGATTGcgcg                                                                                      >  4:122754/1‑149 (MQ=255)
                                                               cGGATCCGGCTTACTGCGAAGGGATTGCAAAGGCTTTGGGGATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCACTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttt                                                                                   <  6:112406/149‑1 (MQ=255)
                                                                 gATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAATAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACAGATTTTTTGCGTCTGATTGCGCGtttt                                                                                  >  7:195552/1‑148 (MQ=255)
                                                                     cGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAg                                                                              <  1:231944/148‑1 (MQ=255)
                                                                      ggCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGt                                                                             >  1:130954/1‑148 (MQ=255)
                                                                        cTTACGGAGAATGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGAc                                                                          <  5:301977/149‑1 (MQ=255)
                                                                         ttACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGAc                                                                          >  6:93317/1‑148 (MQ=255)
                                                                           aCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCa                                                                        >  4:184440/1‑148 (MQ=255)
                                                                           aCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCa                                                                        >  5:96516/1‑148 (MQ=255)
                                                                            cGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATa                                                                      <  4:210043/149‑1 (MQ=255)
                                                                              gCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                    <  5:193360/149‑1 (MQ=255)
                                                                                   gggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                               <  1:359701/149‑1 (MQ=255)
                                                                                      aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCt                                                            >  3:223773/1‑149 (MQ=255)
                                                                                      aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCt                                                            <  5:93317/149‑1 (MQ=255)
                                                                                          cAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaa                                                                 <  8:163719/140‑1 (MQ=255)
                                                                                          cAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaa                                                                 >  7:163719/1‑140 (MQ=255)
                                                                                              ggCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcaca                                                    <  4:155246/149‑1 (MQ=255)
                                                                                                     ggCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtgtg                                             <  6:33138/149‑1 (MQ=255)
                                                                                                         tCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                          >  1:346704/1‑148 (MQ=255)
                                                                                                         tCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCAt                                         <  2:264351/149‑1 (MQ=255)
                                                                                                            aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTTGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                      <  6:27202/149‑1 (MQ=255)
                                                                                                            aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                      >  5:367752/1‑149 (MQ=255)
                                                                                                            aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                      >  8:48834/1‑149 (MQ=255)
                                                                                                               tcGCCTAAGTCGAAGTGATAAGAAGAACCGCCGTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGt                                   >  2:24623/1‑149 (MQ=255)
                                                                                                                 gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACg                                    <  2:130954/146‑1 (MQ=255)
                                                                                                                   cTAAGTCGAAGTGATAAGAAGAACCGCCGTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTc                               <  6:367752/149‑1 (MQ=255)
                                                                                                                      aGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTcac                             <  6:198519/148‑1 (MQ=255)
                                                                                                                          gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAATAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATg                          >  5:336588/1‑147 (MQ=255)
                                                                                                                           aaGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATg                          >  6:172143/1‑146 (MQ=255)
                                                                                                                               gATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                    <  3:184440/148‑1 (MQ=255)
                                                                                                                                aTAAGAAGAACCGCCCTCATTTGGGCGGTATTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                    >  1:215575/1‑148 (MQ=255)
                                                                                                                                 tAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAg                 <  4:244354/149‑1 (MQ=255)
                                                                                                                                  ttgaaaaaCCGCCCAAATGAGGGCGGTTCTTCAATGAATATCACTTCTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                >  5:58760/3‑149 (MQ=255)
                                                                                                                                   agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg               <  2:346704/149‑1 (MQ=255)
                                                                                                                                   agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg               <  6:179433/149‑1 (MQ=255)
                                                                                                                                   agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                <  4:223773/148‑1 (MQ=255)
                                                                                                                                    gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc              >  4:176287/1‑149 (MQ=255)
                                                                                                                                    gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc              >  5:182839/1‑149 (MQ=255)
                                                                                                                                        aaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGAt           >  7:132922/1‑148 (MQ=255)
                                                                                                                                         aCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATg          >  6:110153/1‑148 (MQ=255)
                                                                                                                                           cGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCaa       <  7:83231/149‑1 (MQ=255)
                                                                                                                                            gCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAg      <  3:176287/149‑1 (MQ=255)
                                                                                                                                            gCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAg      <  1:23041/149‑1 (MQ=255)
                                                                                                                                            gCCCTCATTTGGGCGGTTTTTCAATGAAAATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAg      >  1:210781/1‑149 (MQ=255)
                                                                                                                                                tCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATCGAGAATCAGGGCGATGCAAGGcaa  >  2:274342/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCA  >  NC_002947/565842‑566134

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: