Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,125,407 C→A 87.9% H549N (CAT→AAT)  katG → catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,125,4070CA87.9% 457.8 / 64.8 191H549N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base A (75/93);  total (86/105)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 191 total aligned reads displayed.

GGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCG  >  CP009273/4125262‑4125551
                                                                                                                                                 |                                                                                                                                                
ggATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAAt                                                                                                                                                >  2:1637991/1‑148 (MQ=255)
  aTGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGt                                                                                                                                              >  1:908622/1‑148 (MQ=255)
  aTGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGt                                                                                                                                              <  2:2512974/148‑1 (MQ=255)
    gtgAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTAc                                                                                                                                            <  1:381205/148‑1 (MQ=255)
       aaCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGt                                                                                                                                         <  1:2910569/148‑1 (MQ=255)
        aCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGtt                                                                                                                                        <  2:2481679/148‑1 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCggg                                                                                                                               >  2:198124/1‑148 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCggg                                                                                                                               >  2:1584922/1‑148 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCggg                                                                                                                               >  2:1585002/1‑148 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCggg                                                                                                                               >  1:147147/1‑148 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCggg                                                                                                                               >  2:2723972/1‑148 (MQ=255)
                       gTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGt                                                                                                                         <  2:569924/148‑1 (MQ=255)
                        tGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGtt                                                                                                                        <  2:3007008/148‑1 (MQ=255)
                           tctGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGAt                                                                                                                     <  1:2621942/148‑1 (MQ=255)
                               ccTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATgcgc                                                                                                                 >  2:638229/1‑148 (MQ=255)
                                 tGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGt                                                                                                               <  2:1858687/148‑1 (MQ=255)
                                 tGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGt                                                                                                               >  2:2244107/1‑148 (MQ=255)
                                  gTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTc                                                                                                              >  1:3084078/1‑148 (MQ=255)
                                     cTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAgg                                                                                                           <  1:1026803/148‑1 (MQ=255)
                                      tGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGa                                                                                                          <  2:782897/148‑1 (MQ=255)
                                       ggAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGAt                                                                                                         >  1:768624/1‑148 (MQ=255)
                                        gagaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATc                                                                                                        <  2:2935430/148‑1 (MQ=255)
                                         agaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCa                                                                                                       <  1:413861/148‑1 (MQ=255)
                                          gaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAg                                                                                                      <  2:1793941/148‑1 (MQ=255)
                                          gaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAg                                                                                                      <  1:1337799/148‑1 (MQ=255)
                                           aaaaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGa                                                                                                     <  2:1524232/148‑1 (MQ=255)
                                            aaaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGAc                                                                                                    <  2:2635191/148‑1 (MQ=255)
                                             aaTCCAGAAAGAGTCTGGTAAAGCCTGGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACt                                                                                                   >  2:63002/1‑148 (MQ=255)
                                             aaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACt                                                                                                   >  1:762222/1‑148 (MQ=255)
                                                ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                <  2:2228871/148‑1 (MQ=255)
                                                 cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACa                                                                                               <  2:768624/148‑1 (MQ=255)
                                                 cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACa                                                                                               <  2:2647360/148‑1 (MQ=255)
                                                   gAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAtt                                                                                             >  1:359432/1‑148 (MQ=255)
                                                    aaaGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTg                                                                                            <  1:1308301/148‑1 (MQ=255)
                                                      agagTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgag                                                                                          >  2:2866584/1‑148 (MQ=255)
                                                       gagTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgaga                                                                                         >  2:1262748/1‑148 (MQ=255)
                                                        agTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGAt                                                                                        <  1:1211814/148‑1 (MQ=255)
                                                          tCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGt                                                                                      <  2:1173970/148‑1 (MQ=255)
                                                               tAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAg                                                                                 >  1:2987430/1‑148 (MQ=255)
                                                                aaaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCTGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgc                                                                                <  2:2816769/148‑1 (MQ=255)
                                                                aaaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgc                                                                                >  1:1849650/1‑148 (MQ=255)
                                                                 aaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgct                                                                               >  1:1667557/1‑148 (MQ=255)
                                                                  aggcTCGCTGGCGTATATGATTGTGGTGGGTGGTGTTGTTGGTTTTGAGAAAGCCGCAAGCGCCGCTGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCATGATCAGACTGACATTGAGATGTTTGAgctg                                                                              <  1:1145917/145‑1 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  2:3099879/1‑148 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  2:111285/1‑148 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  2:111044/1‑148 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  1:2329537/1‑148 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCCCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                              >  1:2814725/1‑148 (MQ=255)
                                                                     cTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGctgctg                                                                           >  1:1059092/1‑148 (MQ=255)
                                                                      tCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTgg                                                                          <  1:831175/148‑1 (MQ=255)
                                                                           ggCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCa                                                                     >  2:1333352/1‑148 (MQ=255)
                                                                               gATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGATTGAGCTGCTGGAGCCAATTg                                                                 <  1:2022018/148‑1 (MQ=255)
                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                <  1:1807726/148‑1 (MQ=255)
                                                                                 tatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCt                                                               <  1:111044/148‑1 (MQ=255)
                                                                                 tatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCt                                                               <  1:111285/148‑1 (MQ=255)
                                                                                 tatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCt                                                               >  1:1147640/1‑148 (MQ=255)
                                                                                  atCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                                              >  1:1298968/1‑148 (MQ=255)
                                                                                  atCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                                              <  1:1194230/148‑1 (MQ=255)
                                                                                  atCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                                              <  1:560245/148‑1 (MQ=255)
                                                                                  atCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTg                                                              <  1:1171047/148‑1 (MQ=255)
                                                                                      tAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACgg                                                          <  1:3305381/148‑1 (MQ=255)
                                                                                       aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                         <  1:1585002/148‑1 (MQ=255)
                                                                                       aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                         <  1:2723972/148‑1 (MQ=255)
                                                                                       aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                         <  1:1584922/148‑1 (MQ=255)
                                                                                       aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                         <  2:654102/148‑1 (MQ=255)
                                                                                        gTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGtt                                                        >  2:237930/1‑148 (MQ=255)
                                                                                        gTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGtt                                                        >  2:2783748/1‑148 (MQ=255)
                                                                                        gTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGtt                                                        <  1:22398/148‑1 (MQ=255)
                                                                                          gctggctgGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTc                                                      <  1:660916/148‑1 (MQ=255)
                                                                                          gctggctgGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTc                                                      <  1:660814/148‑1 (MQ=255)
                                                                                           ctggctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTcc                                                     >  2:3202145/1‑148 (MQ=255)
                                                                                              gctggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTa                                                  >  2:1997772/1‑148 (MQ=255)
                                                                                                 ggTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACt                                               <  1:3170616/148‑1 (MQ=255)
                                                                                                   tgtgGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTAt                                             <  1:2270620/148‑1 (MQ=255)
                                                                                                   tgtgGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTAt                                             >  2:1238496/1‑148 (MQ=255)
                                                                                                       gttggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATcgcg                                         >  1:316732/1‑148 (MQ=255)
                                                                                                         tggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCt                                       >  1:2110052/1‑148 (MQ=255)
                                                                                                         tggtGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCt                                       >  1:440997/1‑148 (MQ=255)
                                                                                                              ttGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCt                                  <  1:2727388/148‑1 (MQ=255)
                                                                                                                 agaAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGa                               <  1:3057901/148‑1 (MQ=255)
                                                                                                                    aaGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGt                            <  2:214202/148‑1 (MQ=255)
                                                                                                                      gCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATTTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                          >  1:2457325/1‑148 (MQ=255)
                                                                                                                      gCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                          <  1:2142651/148‑1 (MQ=255)
                                                                                                                      gCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                          >  1:2824505/1‑148 (MQ=255)
                                                                                                                         gCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTcca                       <  2:78835/148‑1 (MQ=255)
                                                                                                                          cAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTccac                      >  1:1706328/1‑148 (MQ=255)
                                                                                                                            aGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTccacca                    <  2:142425/148‑1 (MQ=255)
                                                                                                                             gcgcCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCcaccac                   >  1:859348/1‑148 (MQ=255)
                                                                                                                             gcgcCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCcaccac                   >  2:1002091/1‑148 (MQ=255)
                                                                                                                              cgccgcAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCaccacc                  >  2:2885069/1‑148 (MQ=255)
                                                                                                                              cgccgcAGGTTTGAGCATTAATGTACCGTTTGCCCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCaccacc                  >  2:941825/1‑148 (MQ=255)
                                                                                                                                    aGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCa            >  1:3241763/1‑148 (MQ=255)
                                                                                                                                     ggTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCAc           >  2:478502/1‑148 (MQ=255)
                                                                                                                                     ggTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCAc           >  1:2095609/1‑148 (MQ=255)
                                                                                                                                      gTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCAct          >  1:1079810/1‑148 (MQ=255)
                                                                                                                                        ttGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActgc        >  2:3054040/1‑148 (MQ=255)
                                                                                                                                           aGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGa     >  1:2542618/1‑148 (MQ=255)
                                                                                                                                           aGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGa     >  1:1716392/1‑148 (MQ=255)
                                                                                                                                              aTTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCg  <  2:146736/148‑1 (MQ=255)
                                                                                                                                              aTTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCg  >  2:324895/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCG  >  CP009273/4125262‑4125551

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: