Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,991 T→A 7.7% intergenic (+31/‑133) katA → / → bfr‑I catalase/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,9910TA 7.7% 267.9 / 9.0 90intergenic (+31/‑133)katA/bfr‑Icatalase/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base T (39/44);  new base A (4/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACC  >  NC_002947/565845‑566136
                                                                                                                                                  |                                                                                                                                                 
gTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTgg                                                                                                                                                 <  5:13425/149‑1 (MQ=255)
   ccATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCg                                                                                                                                              >  1:8433/1‑149 (MQ=255)
    cATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGTGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCgg                                                                                                                                             <  6:13187/149‑1 (MQ=255)
           ggcgTCAGCGAAGATGTCATTCAGCGTCAGCTTCAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCa                                                                                                                                      <  3:231248/149‑1 (MQ=255)
                   cGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAata                                                                                                                               >  5:106704/1‑148 (MQ=255)
                       gATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCac                                                                                                                           >  7:46749/1‑148 (MQ=255)
                             aTTCAGCGTCAGCTACATTAGTTCTTCAAGGAGTTGCCGGCTTACGGAGAAGGGATTGGAAAGGCTTTGGGCATCAATCTCGCCTAATTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                    <  7:249064/149‑1 (MQ=255)
                             aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                                    <  2:274111/149‑1 (MQ=255)
                               tCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAat                                                                                                                                >  6:324826/1‑135 (MQ=255)
                               tCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAat                                                                                                                                <  5:324826/135‑1 (MQ=255)
                                 aGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAc                                                                                                                <  2:8433/149‑1 (MQ=255)
                                      cAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGttttt                                                                                                                                        <  4:316575/120‑1 (MQ=255)
                                      cAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGttttt                                                                                                                                        >  3:316575/1‑120 (MQ=255)
                                       aGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAttt                                                                                                          >  3:229271/1‑149 (MQ=255)
                                          tACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAtttttt                                                                                                       >  1:310562/1‑149 (MQ=255)
                                            cAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGc                                                                                                     >  4:42483/1‑149 (MQ=255)
                                            cAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGc                                                                                                     >  4:189943/1‑149 (MQ=255)
                                               tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                  <  2:121238/149‑1 (MQ=255)
                                               tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGACATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCggt                                                                                                  >  7:32496/1‑147 (MQ=255)
                                                  ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                               >  1:180599/1‑149 (MQ=255)
                                                    cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAcc                                                                                                               >  7:116342/1‑131 (MQ=255)
                                                    cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAcc                                                                                                               <  8:116342/131‑1 (MQ=255)
                                                    cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAc                                                                                                                <  3:127237/130‑1 (MQ=255)
                                                    cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAAc                                                                                                                >  4:127237/1‑130 (MQ=255)
                                                      tcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                                 <  7:145635/143‑1 (MQ=255)
                                                      tcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                                 >  8:145635/1‑143 (MQ=255)
                                                      tcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                           >  8:173735/1‑149 (MQ=255)
                                                        aaGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgcg                                                                                          <  3:245162/148‑1 (MQ=255)
                                                         aGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGcgcg                                                                                        <  2:180599/149‑1 (MQ=255)
                                                              gATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                   <  3:42483/149‑1 (MQ=255)
                                                                 ccGGCGTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAACTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAg                                                                                <  8:267588/149‑1 (MQ=255)
                                                                   ggCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTc                                                                                  >  4:31237/1‑145 (MQ=255)
                                                                   ggCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTc                                                                                  <  3:31237/145‑1 (MQ=255)
                                                                    gCTTACGGCGAAGGGATTGCATAGGCTTTGGGCATAAATCTAGACTAATTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCTATGAATATCACTACTGTTTAACCGAtttt                                                                                                         <  5:357121/121‑1 (MQ=255)
                                                                    gCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATAAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAATTGGGCGATTTTTCAATGAATATCAAAACTGTTTAAACGAtttt                                                                                                         >  6:357121/1‑121 (MQ=255)
                                                                          ggCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAg                                                                       <  4:68335/149‑1 (MQ=255)
                                                                           gataaGTGATTGCAAAGTCTTTGTGCATCAATCTAGCCTAAGTCGTAGTGATAAGTAGAACCGGCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                      <  5:127363/146‑1 (MQ=255)
                                                                           gCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                      <  7:173735/149‑1 (MQ=255)
                                                                             gAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCACTCATTTGGGCGGTTTTTCAATGAATTTCAAGACTGTTTAACCGATTTTTTGCGTCTGATTGAGCGTTTTTCAGTGACCTTAGAc                                                                     >  6:351261/1‑148 (MQ=255)
                                                                                gggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaa                                                                   >  3:111455/1‑147 (MQ=255)
                                                                                 ggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGAcc                                                                           >  7:143320/1‑138 (MQ=255)
                                                                                 ggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGAcc                                                                           <  8:143320/138‑1 (MQ=255)
                                                                                  gATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGc                                                               >  7:227735/1‑149 (MQ=255)
                                                                                  gATTGCAAAGGCTTTGGGCATCAATCTAGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAg                                                                <  7:123555/148‑1 (MQ=255)
                                                                                      gCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                            <  2:310562/148‑1 (MQ=255)
                                                                                      gCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                            <  4:229271/148‑1 (MQ=255)
                                                                                      gCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                            <  6:106704/148‑1 (MQ=255)
                                                                                          aGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcac                                                       <  5:363830/149‑1 (MQ=255)
                                                                                                tGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAt                                                  <  8:46749/148‑1 (MQ=255)
                                                                                                  ggCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtatg                                               <  3:189943/149‑1 (MQ=255)
                                                                                                   gCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTgg                                              >  3:4180/1‑149 (MQ=255)
                                                                                                   gCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTgg                                              <  8:227735/149‑1 (MQ=255)
                                                                                                     aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                   <  7:60065/110‑1 (MQ=255)
                                                                                                     aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                   >  8:60065/1‑110 (MQ=255)
                                                                                                             cGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcaca                                                      <  7:113884/131‑1 (MQ=255)
                                                                                                               ccTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                  >  5:58489/1‑149 (MQ=255)
                                                                                                                cTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                        >  4:192695/1‑142 (MQ=255)
                                                                                                                cTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                        <  3:192695/142‑1 (MQ=255)
                                                                                                                cTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCCATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                  >  8:105067/1‑148 (MQ=255)
                                                                                                                 tAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCAt                                                                         >  6:233225/1‑108 (MQ=255)
                                                                                                                 tAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCAt                                                                         <  5:233225/108‑1 (MQ=255)
                                                                                                                      cGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATg                            >  2:193338/1‑148 (MQ=255)
                                                                                                                       gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                 >  3:344580/1‑94 (MQ=255)
                                                                                                                       gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                 <  4:344580/94‑1 (MQ=255)
                                                                                                                       gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                          >  8:185086/1‑149 (MQ=255)
                                                                                                                         aGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgaga                        <  4:4180/149‑1 (MQ=255)
                                                                                                                          gTGATAAGAATAACCGCCCTCATTTGTGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGaa                       <  7:105067/149‑1 (MQ=255)
                                                                                                                          attcattgaaaaaCCGCCCAAATGAGGGCGGTTCTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGaa                       >  7:316682/8‑149 (MQ=38)
                                                                                                                            gATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                      >  7:174549/1‑148 (MQ=255)
                                                                                                                              tAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAg                   <  6:58489/149‑1 (MQ=255)
                                                                                                                               aagaagAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                  >  8:115538/1‑149 (MQ=255)
                                                                                                                                agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                  <  7:115538/148‑1 (MQ=255)
                                                                                                                                 gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc                >  6:313062/1‑149 (MQ=255)
                                                                                                                                  aagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                       >  2:100677/1‑125 (MQ=255)
                                                                                                                                  aagaaCCGCCCTAATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                       <  1:100677/125‑1 (MQ=255)
                                                                                                                                  aaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtt                                    <  7:2270/125‑1 (MQ=255)
                                                                                                                                  aaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtt                                    >  8:2270/4‑128 (MQ=255)
                                                                                                                                  aaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc                >  2:220416/4‑148 (MQ=255)
                                                                                                                                     aaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                       >  8:273826/1‑122 (MQ=255)
                                                                                                                                     aaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                       <  7:273826/122‑1 (MQ=255)
                                                                                                                                        cGCCCTCATTTGGGCGTTTTTTCAATGAATATCACTACTGTTTAACAGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGTTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCaa         <  5:351261/149‑1 (MQ=255)
                                                                                                                                        cGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCCCACTATGTTGCATAGACGTTTTCACATGGAGAATCAGGGCGATGc           >  7:194448/1‑147 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTc                                 <  2:318141/123‑1 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTc                                 >  1:318141/1‑123 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAgg       <  5:313062/149‑1 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAgg       <  1:193338/149‑1 (MQ=255)
                                                                                                                                           ccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                  <  2:301103/57‑1 (MQ=255)
                                                                                                                                           ccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                  >  1:301103/1‑57 (MQ=255)
                                                                                                                                             tCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCa    >  6:347819/1‑149 (MQ=255)
                                                                                                                                               aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCAcc  >  6:256768/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACC  >  NC_002947/565845‑566136

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: