Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,993 T→G 5.2% 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,9930TG 5.2% 1131.7 / 36.8 344intergenic (+30/‑134)katA/bfr‑Icatalase/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base T (159/167);  new base G (11/7);  total (170/174)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 351 total aligned reads displayed.

GCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACCCGG  >  NC_002947/565844‑566142
                                                                                                                                                     |                                                                                                                                                     
gCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAtt                                                                                                                                                       <  1:4540855/150‑1 (MQ=255)
gCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCAtt                                                                                                                                                       >  2:3896121/1‑150 (MQ=255)
        cATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGt                                                                                                                                               <  2:4835251/150‑1 (MQ=255)
           ggcgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCTAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGtttt                                                                                                                                            <  2:1109499/150‑1 (MQ=255)
           ggcgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGtttt                                                                                                                                            >  2:767564/1‑150 (MQ=255)
             cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTc                                                                                                                                          <  1:2175203/150‑1 (MQ=255)
             cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttcttc                                                                                                                                          <  1:1310728/150‑1 (MQ=255)
              gggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCa                                                                                                                                         <  1:6348021/150‑1 (MQ=255)
               ggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGATGGGATTGCAAAGGCTTTGGGCATCTATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttcttctt                                                                                                                                        <  1:295127/150‑3 (MQ=255)
                gcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAAt                                                                                                                                       <  1:4043280/150‑1 (MQ=255)
                gcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAAt                                                                                                                                       <  1:2250514/150‑1 (MQ=255)
                  gTCAGCGAAGATGTCATTCAGCGTCAGATACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGa                                                                                                                                     <  1:6446819/150‑1 (MQ=255)
                    cAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAat                                                                                                                                   <  1:2326974/150‑1 (MQ=255)
                     aGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAata                                                                                                                                  >  1:5794081/1‑150 (MQ=255)
                      gCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAatat                                                                                                                                 >  1:4407087/1‑150 (MQ=255)
                        gAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCa                                                                                                                               <  2:5789177/150‑1 (MQ=255)
                              gTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTg                                                                                                                         >  2:2065580/1‑150 (MQ=255)
                              gTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTg                                                                                                                         >  1:5235332/1‑150 (MQ=255)
                                 aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGttt                                                                                                                      >  2:2643773/1‑150 (MQ=255)
                                 aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGttt                                                                                                                      >  1:5092540/1‑150 (MQ=255)
                                             cTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAtttttt                                                                                                          <  1:1154001/150‑1 (MQ=255)
                                              tACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTg                                                                                                         >  1:4120563/1‑150 (MQ=255)
                                               aCAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGc                                                                                                        <  2:4033700/150‑1 (MQ=255)
                                                 attACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGt                                                                                                      <  1:551625/148‑1 (MQ=255)
                                                 aGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGt                                                                                                      >  2:2177593/1‑150 (MQ=255)
                                                   tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAGGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                                    >  2:2925903/1‑150 (MQ=255)
                                                   tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                                    >  1:5341106/1‑150 (MQ=255)
                                                   tACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCt                                                                                                    >  1:6848131/1‑150 (MQ=255)
                                                    aCTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTg                                                                                                   >  2:1391435/1‑150 (MQ=255)
                                                      ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAt                                                                                                 >  1:2428125/1‑150 (MQ=255)
                                                      ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAt                                                                                                 >  2:4931113/1‑150 (MQ=255)
                                                       tcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                                <  1:6168457/150‑1 (MQ=255)
                                                         ttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                              >  1:322255/1‑150 (MQ=255)
                                                           cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgcgc                                                                                            >  2:1932820/1‑150 (MQ=255)
                                                           cAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgcgc                                                                                            <  1:3658472/150‑1 (MQ=255)
                                                            aaGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGcgcg                                                                                           >  2:2172429/1‑150 (MQ=255)
                                                               gCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttt                                                                                        >  2:3320533/1‑150 (MQ=255)
                                                               gCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttt                                                                                        >  2:2866550/1‑150 (MQ=255)
                                                               gCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttt                                                                                        <  1:5927842/150‑1 (MQ=255)
                                                                cGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGtttt                                                                                       <  1:377651/150‑1 (MQ=255)
                                                                 ggATCCGGCTTTCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                      <  1:1027055/150‑1 (MQ=255)
                                                                 ggATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                      <  1:6034094/150‑1 (MQ=255)
                                                                 ggATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                      <  1:2023144/150‑1 (MQ=255)
                                                                 ggATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                      <  2:3160940/150‑1 (MQ=255)
                                                                   aTCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                    <  2:5989125/150‑1 (MQ=255)
                                                                          ttACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCa                                                                             <  2:43614/150‑1 (MQ=255)
                                                                              ggCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                         <  2:6076884/150‑1 (MQ=255)
                                                                                 gAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaa                                                                      <  1:1525071/150‑1 (MQ=255)
                                                                                   aGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAGGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                                    <  1:2925903/150‑1 (MQ=255)
                                                                                   aGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaaa                                                                    <  1:3896495/150‑1 (MQ=255)
                                                                                      gATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                         <  2:6847480/142‑1 (MQ=255)
                                                                                       aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTg                                                                <  2:5972686/150‑1 (MQ=255)
                                                                                       aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTg                                                                <  2:6398322/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                               <  1:792325/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                               <  2:6551387/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                               <  1:179711/150‑1 (MQ=255)
                                                                                           cAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTc                                                            >  1:274530/1‑150 (MQ=255)
                                                                                           cAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTc                                                            >  1:3586825/1‑150 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                           <  2:5329752/150‑1 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                           >  1:1952089/1‑150 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                           >  1:4672854/1‑150 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                           >  1:2493095/1‑150 (MQ=255)
                                                                                              aGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcaca                                                         >  2:3797879/1‑150 (MQ=255)
                                                                                              aGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTCcaca                                                         <  2:1321755/150‑1 (MQ=255)
                                                                                                     gggCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAtgtg                                                  <  2:1272647/150‑1 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                               >  2:4480917/1‑150 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                               >  2:2447070/1‑150 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                               >  2:3749052/1‑150 (MQ=255)
                                                                                                         aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAGCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCAt                                              >  2:31189/1‑150 (MQ=255)
                                                                                                                  gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                     <  1:978630/150‑1 (MQ=255)
                                                                                                                  gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                     >  2:3693872/1‑150 (MQ=255)
                                                                                                                   ccTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGATTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTc                                    <  2:5088721/150‑1 (MQ=255)
                                                                                                                     tAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGACTATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTcac                                  >  2:1497881/1‑150 (MQ=255)
                                                                                                                       aGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACAt                                <  2:6715780/150‑1 (MQ=255)
                                                                                                                         tCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATgg                              <  1:1445718/150‑1 (MQ=255)
                                                                                                                          cGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                             <  2:6237032/150‑1 (MQ=255)
                                                                                                                          cGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                             <  1:3837367/150‑1 (MQ=255)
                                                                                                                          cGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                             <  1:2848927/150‑1 (MQ=255)
                                                                                                                            aaGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgaga                           <  2:5535232/150‑1 (MQ=255)
                                                                                                                            aaGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgaga                           <  1:671581/150‑1 (MQ=255)
                                                                                                                                 aTAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAg                      >  1:828976/1‑150 (MQ=255)
                                                                                                                                  tAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                     <  2:1824626/150‑1 (MQ=255)
                                                                                                                                  tAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAgg                     <  1:6103102/150‑1 (MQ=255)
                                                                                                                                   aagaagAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACGCTATGTGGCATAGACGTTTTCACATGGAGAATCAggg                    >  1:1106356/1‑150 (MQ=255)
                                                                                                                                   aagaagAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg                    <  2:4120563/150‑1 (MQ=255)
                                                                                                                                   aagaagAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg                    >  1:4851895/1‑150 (MQ=255)
                                                                                                                                   aagaagAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg                    >  1:6747398/1‑150 (MQ=255)
                                                                                                                                    agaagaACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGc                   >  2:3569870/1‑150 (MQ=255)
                                                                                                                                     gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCg                  >  2:5575746/1‑150 (MQ=255)
                                                                                                                                     gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCg                  >  1:3367551/1‑150 (MQ=255)
                                                                                                                                       aaaaCCGCCCAAATGAGGGCGGTTCTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAAGTTTTCACATGGAGAATCAGGTCGAt                >  2:3342077/3‑150 (MQ=255)
                                                                                                                                        aaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATg               >  1:2529322/2‑150 (MQ=255)
                                                                                                                                         aaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGc              >  2:45784/1‑150 (MQ=255)
                                                                                                                                         aaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGc              >  1:305558/1‑150 (MQ=255)
                                                                                                                                           ccGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCaa            <  1:2607961/150‑1 (MQ=255)
                                                                                                                                                cTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCa       >  2:1748711/1‑150 (MQ=255)
                                                                                                                                                   aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgaga                           >  1:5010857/1‑127 (MQ=255)
                                                                                                                                                   aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCAccc    <  2:5713991/150‑1 (MQ=255)
                                                                                                                                                   aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCAccc    >  1:5687736/1‑150 (MQ=255)
                                                                                                                                                   aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCAccc    >  1:1708037/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACCCGG  >  NC_002947/565844‑566142

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