Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,990 T→G 8.5% intergenic (+30/‑134) katA → / → bfr‑I catalase/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,9900TG 8.5% 192.5 / 9.9 71intergenic (+30/‑134)katA/bfr‑Icatalase/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base T (28/36);  new base G (4/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

GTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGG  >  NC_002947/565845‑566131
                                                                                                                                                 |                                                                                                                                             
gTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTgg                                                                                                                                            >  7:75349/1‑149 (MQ=255)
 tACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCTGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTggg                                                                                                                                           <  5:4365/149‑1 (MQ=255)
 tACCATGCAGGGCGTCCGCGAACAAGTCATTCAGCGCCAGCTACCGTACTTCTTCAAGGCGGAGCCGTCTTACGGCGAAGGGAGTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAGTTggc                                                                                                                                           <  8:281305/149‑2 (MQ=255)
             cgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAAt                                                                                                                               <  7:145351/149‑1 (MQ=255)
                 ggCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttcttcttatca                                                                                                                            <  1:311424/147‑7 (MQ=255)
                  gCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAata                                                                                                                          >  6:251749/1‑149 (MQ=255)
                     aaGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCa                                                                                                                       <  1:92269/149‑1 (MQ=255)
                         tGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCacta                                                                                                                    <  8:75349/148‑1 (MQ=255)
                             aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGtt                                                                                                               >  2:117062/1‑149 (MQ=255)
                                    gTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGa                                                                                                        >  4:159854/1‑149 (MQ=255)
                                    gTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTAAATGAATATCACTAATGTTTAACCGa                                                                                                        >  3:46420/1‑149 (MQ=255)
                                           aCAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTg                                                                                                 >  1:181661/1‑149 (MQ=255)
                                            cAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGAATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGc                                                                                                >  8:64140/1‑149 (MQ=255)
                                             aGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTgag                                                                                               >  5:244384/1‑147 (MQ=255)
                                                  ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                          >  5:217010/1‑149 (MQ=255)
                                                  ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                          >  8:274842/1‑149 (MQ=255)
                                                    cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                        >  4:287951/1‑149 (MQ=255)
                                                      tcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                      >  2:168092/1‑149 (MQ=255)
                                                       cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                        <  6:128215/146‑1 (MQ=255)
                                                       cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                      <  5:251749/148‑1 (MQ=255)
                                                         aGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACGATTTAACCGATTAATTGCGTCAGATTGcgcg                                                                                   >  7:112727/1‑149 (MQ=255)
                                                           gctgATCAGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGAATCAATCTCGCCTAAGTCGAAGTTATAAGAAGAACCGCCCTCATTTGGGc                                                                                                                                          <  8:286707/89‑1 (MQ=255)
                                                           gCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGc                                                                                                                                          >  7:286707/1‑92 (MQ=255)
                                                            cGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttt                                                                                <  7:179841/149‑1 (MQ=255)
                                                                tCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGTTTTTTCAATTAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                            >  3:38608/1‑149 (MQ=255)
                                                                tCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                              >  6:146301/1‑147 (MQ=255)
                                                                    gCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGt                                                                          <  8:23004/147‑1 (MQ=255)
                                                                    gCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTAAGTGa                                                                        >  7:23004/1‑149 (MQ=255)
                                                                     cTTACTGCGATGGGTTTGCAAAGGCTTTGGGTTTCAATCTCGCCTAAGTCGTTTTGATTAGTTGAACCGCTCTGATTTGGGCGGTTTTTCAATGATTTTGACTTCTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGAc                                                                       <  8:112727/149‑1 (MQ=255)
                                                                        aCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCa                                                                     >  6:279064/1‑148 (MQ=255)
                                                                           gCGAAGGGATTGCAAATGCTTTGGGCATCAATCTCGCCTTAGTCTAAGTTATTAGAAGAACCGCCCTCATTTTGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                 <  7:130239/149‑1 (MQ=255)
                                                                           gCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGt                                                                                                                <  7:70057/102‑1 (MQ=255)
                                                                           gCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGt                                                                                                                >  8:70057/1‑102 (MQ=255)
                                                                                gggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                            <  1:117062/149‑1 (MQ=255)
                                                                                gggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                            <  2:323918/149‑1 (MQ=255)
                                                                                gggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                            <  2:4345/149‑1 (MQ=255)
                                                                                        aaaGGCTTTGGCCATCAATCTCGCCTAACTCGAAGAGATAAGAAGAACCGCCCTCATATGTGCCGTTTTTCAATGAATATCACTACTGTTTAACCGATTTAAAGCGACTGATTGCCCGTTTTTCAGTGACCATa                                                                   >  2:7605/1‑134 (MQ=255)
                                                                                           ggCTTTGGTCATCAATCTAGCCTAAGTCGAAGTGATAAGAAGAACCGCGCTCATTTGGGAGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcaca                                                 <  8:153291/149‑1 (MQ=255)
                                                                                                  ggCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtgt                                           <  1:168092/148‑1 (MQ=255)
                                                                                                    cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTgg                                         >  3:16285/1‑148 (MQ=255)
                                                                                                    cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCACATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGc                                        <  5:279064/149‑1 (MQ=255)
                                                                                                     aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtgtg                                          <  4:250466/146‑1 (MQ=255)
                                                                                                     aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtgtg                                          >  3:250466/1‑146 (MQ=255)
                                                                                                     aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                       <  7:274842/149‑1 (MQ=255)
                                                                                                       cAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGAGTTTTTAATTGACCATAGACAAAAGCTGGTTCACAATATGTGGCa                                       >  8:151052/1‑147 (MQ=255)
                                                                                                              gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGttt                              <  6:217010/149‑1 (MQ=255)
                                                                                                               ccTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                             <  4:115899/149‑1 (MQ=255)
                                                                                                                 tAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                   >  4:289148/1‑141 (MQ=255)
                                                                                                                 tAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                   <  3:289148/141‑1 (MQ=255)
                                                                                                                     tCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTAAATGAATATCACTACTGTTTAACCGATTTTTTGCGTATGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATg                       <  6:244384/149‑1 (MQ=255)
                                                                                                                       gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTAAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATg                       >  3:160558/1‑147 (MQ=255)
                                                                                                                       gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTATGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                     <  4:5463/149‑1 (MQ=255)
                                                                                                                         aGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgaga                   <  2:181661/149‑1 (MQ=255)
                                                                                                                            gATAATAATAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTTTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                 <  7:151052/148‑1 (MQ=255)
                                                                                                                            gATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                       <  1:212446/110‑1 (MQ=255)
                                                                                                                            gATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                       >  2:212446/1‑110 (MQ=255)
                                                                                                                             aTAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCactac                                                                                                                   >  7:39140/1‑49 (MQ=255)
                                                                                                                             aTAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCactac                                                                                                                   <  8:39140/49‑1 (MQ=255)
                                                                                                                              tAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAg              <  5:146301/149‑1 (MQ=255)
                                                                                                                               ttgaaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                     <  1:161966/139‑1 (MQ=255)
                                                                                                                               ttgaaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                     >  2:161966/3‑141 (MQ=255)
                                                                                                                                agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg            <  3:287951/149‑1 (MQ=255)
                                                                                                                                 gaaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc           >  7:131182/1‑149 (MQ=255)
                                                                                                                                   agaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGa         >  5:204503/1‑149 (MQ=255)
                                                                                                                                   aaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGa         >  8:177068/3‑149 (MQ=255)
                                                                                                                                     aaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATg       >  1:152342/1‑149 (MQ=255)
                                                                                                                                      aCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGc      <  3:159854/149‑1 (MQ=255)
                                                                                                                                         gCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAg   <  7:64140/149‑1 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCaa    >  2:273527/1‑147 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCaa    <  1:273527/147‑1 (MQ=255)
                                                                                                                                          cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAgg  <  6:315571/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGG  >  NC_002947/565845‑566131

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: ‑