Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,125,407 C→A 100% 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,4070CA95.8% 570.3 / 15.3 190H549N (CAT→AAT) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (95/87);  total (100/90)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 190 total aligned reads displayed.

GGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCGAC  >  CP009273/4125262‑4125553
                                                                                                                                                 |                                                                                                                                                  
   tgtgAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTa                                                                                                                                               <  2:1447213/148‑1 (MQ=255)
    gtgAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTAc                                                                                                                                              <  1:2921986/148‑1 (MQ=255)
      gAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCg                                                                                                                                            >  2:542214/1‑148 (MQ=255)
      gAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCg                                                                                                                                            <  1:351947/148‑1 (MQ=255)
      gAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGGGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCg                                                                                                                                            >  2:1993275/1‑148 (MQ=255)
       aaCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGt                                                                                                                                           <  2:1138770/148‑1 (MQ=255)
        aCGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGtt                                                                                                                                          <  2:1901318/148‑1 (MQ=255)
               aGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATATTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAATCGCCGCAGGTTTGAGCATTAATGTAACGTTTGCGCCg                                                                                                                                   >  1:465041/1‑148 (MQ=255)
                 ccGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCggg                                                                                                                                 >  2:1482069/1‑148 (MQ=255)
                   gTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTc                                                                                                                               >  1:2669301/1‑148 (MQ=255)
                    ttCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTcg                                                                                                                              >  2:2299424/1‑148 (MQ=255)
                    ttCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTcg                                                                                                                              >  2:2299382/1‑148 (MQ=255)
                        tGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGtt                                                                                                                          <  1:283659/148‑1 (MQ=255)
                           tctGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGAt                                                                                                                       <  2:1041684/148‑1 (MQ=255)
                                 tGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGt                                                                                                                 <  2:3251176/148‑1 (MQ=255)
                                 tGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGt                                                                                                                 >  2:391014/1‑148 (MQ=255)
                                       ggAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGAt                                                                                                           <  1:3151323/148‑1 (MQ=255)
                                       ggAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGAt                                                                                                           <  1:2840658/148‑1 (MQ=255)
                                        gagaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATc                                                                                                          >  2:62854/1‑148 (MQ=255)
                                        gagaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATc                                                                                                          <  1:2679689/148‑1 (MQ=255)
                                         agaAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCa                                                                                                         >  1:1421973/1‑148 (MQ=255)
                                           aaaaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCATGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGa                                                                                                       >  2:1413941/1‑148 (MQ=255)
                                             aaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACt                                                                                                     >  2:2983894/1‑148 (MQ=255)
                                             aaTCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTCGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACt                                                                                                     <  1:1058607/148‑1 (MQ=255)
                                               tCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGa                                                                                                   <  1:3103753/148‑1 (MQ=255)
                                                ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                  >  1:305230/1‑148 (MQ=255)
                                                ccAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGATTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                  >  1:2116749/1‑148 (MQ=255)
                                                 cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACa                                                                                                 <  2:2910864/148‑1 (MQ=255)
                                                 cAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACa                                                                                                 >  1:1630908/1‑148 (MQ=255)
                                                   gAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACAtt                                                                                               >  1:2683055/1‑148 (MQ=255)
                                                     aaGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTga                                                                                             >  2:3423006/1‑148 (MQ=255)
                                                      agagTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgag                                                                                            >  1:1918747/1‑148 (MQ=255)
                                                      agagTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTgag                                                                                            <  1:2369427/148‑1 (MQ=255)
                                                        agTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGAt                                                                                          >  1:1158501/1‑148 (MQ=255)
                                                         gTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCTTCAGGATCAGACTGACATTGAGATg                                                                                         <  2:3376188/148‑1 (MQ=255)
                                                         gTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATg                                                                                         >  2:1648272/1‑148 (MQ=255)
                                                          tCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGt                                                                                        <  2:2947979/148‑1 (MQ=255)
                                                             ggTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTg                                                                                     <  2:1090055/148‑1 (MQ=255)
                                                             ggTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTg                                                                                     >  2:1278632/1‑148 (MQ=255)
                                                               tAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAg                                                                                   >  2:2526403/1‑148 (MQ=255)
                                                               tAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAg                                                                                   <  2:2596742/148‑1 (MQ=255)
                                                                 aaGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgct                                                                                 >  1:3004958/1‑148 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                                <  2:1451233/148‑1 (MQ=255)
                                                                  aGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctg                                                                                <  1:1659466/148‑1 (MQ=255)
                                                                   gCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctgc                                                                               >  1:2528996/1‑148 (MQ=255)
                                                                   gCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctgc                                                                               >  1:3135132/1‑148 (MQ=255)
                                                                    ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                  <  2:3177409/128‑1 (MQ=255)
                                                                    ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGAc                                                                                                  >  1:3177409/1‑128 (MQ=255)
                                                                    ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctgct                                                                              >  1:709117/1‑148 (MQ=255)
                                                                    ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctgct                                                                              >  1:648978/1‑148 (MQ=255)
                                                                    ccTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAgctgct                                                                              >  1:507202/1‑148 (MQ=255)
                                                                      tCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTgg                                                                            <  1:1605460/148‑1 (MQ=255)
                                                                      tCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTgg                                                                            <  2:2713303/148‑1 (MQ=255)
                                                                      tCGCTGGCGGAGATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCTCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTgg                                                                            >  2:825065/1‑148 (MQ=255)
                                                                        gCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAg                                                                          <  2:649120/148‑1 (MQ=255)
                                                                          tGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGcc                                                                        >  1:1035799/1‑148 (MQ=255)
                                                                          tGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGcc                                                                        <  1:2588546/148‑1 (MQ=255)
                                                                             cGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAAt                                                                     >  2:2316268/1‑148 (MQ=255)
                                                                               gATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTg                                                                   <  2:2651171/148‑1 (MQ=255)
                                                                               gATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTg                                                                   <  2:2651722/148‑1 (MQ=255)
                                                                               gATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTg                                                                   <  2:2172381/148‑1 (MQ=255)
                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                  >  1:592432/1‑148 (MQ=255)
                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                  >  1:2139429/1‑148 (MQ=255)
                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                  >  1:2848446/1‑148 (MQ=255)
                                                                                atatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGc                                                                  >  2:2120036/1‑148 (MQ=255)
                                                                                 tatCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCt                                                                 <  1:542214/148‑1 (MQ=255)
                                                                                   tCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCTGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGa                                                               <  2:2817068/148‑1 (MQ=255)
                                                                                   tCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGa                                                               <  1:1795431/148‑1 (MQ=255)
                                                                                    cATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGAc                                                              >  1:3243773/1‑148 (MQ=255)
                                                                                    cATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGAc                                                              >  1:2343099/1‑148 (MQ=255)
                                                                                      tAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACgg                                                            >  1:3195027/1‑148 (MQ=255)
                                                                                      tAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACgg                                                            >  2:3432036/1‑148 (MQ=255)
                                                                                       aGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGt                                                           <  1:3172606/148‑1 (MQ=255)
                                                                                        gTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGtt                                                          <  2:1076615/148‑1 (MQ=255)
                                                                                                  gtgtGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTa                                                <  1:3127491/148‑1 (MQ=255)
                                                                                                             gtTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTc                                     >  2:736458/1‑148 (MQ=255)
                                                                                                                gagaAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTgg                                  <  1:62640/148‑1 (MQ=255)
                                                                                                                  gaAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGAc                                >  1:3020001/1‑148 (MQ=255)
                                                                                                                  gaAAGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGAc                                >  2:730061/1‑148 (MQ=255)
                                                                                                                   aaaGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACg                               >  1:3275839/1‑148 (MQ=255)
                                                                                                                     aGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGtt                             >  1:726407/1‑148 (MQ=255)
                                                                                                                     aGCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGtt                             >  2:31461/1‑148 (MQ=255)
                                                                                                                      gCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                            <  1:3010538/148‑1 (MQ=255)
                                                                                                                      gCCGCAAGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGttt                            <  1:1233134/148‑1 (MQ=255)
                                                                                                                           aaGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTccacc                       <  2:22939/148‑1 (MQ=255)
                                                                                                                            aGCGCCGCAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTccacca                      >  1:1463605/1‑148 (MQ=255)
                                                                                                                                 cgcAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAg                 >  2:2372598/1‑148 (MQ=255)
                                                                                                                                   cAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTc               >  1:1595861/1‑148 (MQ=255)
                                                                                                                                   cAGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTc               >  1:1379549/1‑148 (MQ=255)
                                                                                                                                    aGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCa              <  1:1267031/148‑1 (MQ=255)
                                                                                                                                    aGGTTTGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCa              <  1:2719174/148‑1 (MQ=255)
                                                                                                                                       tttGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActg           <  2:362146/148‑1 (MQ=255)
                                                                                                                                       tttGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActg           >  1:679607/1‑148 (MQ=255)
                                                                                                                                        ttGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActgc          <  2:2688489/148‑1 (MQ=255)
                                                                                                                                        ttGAGCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCActgc          >  2:313114/1‑148 (MQ=255)
                                                                                                                                            gCATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGAt      <  1:2560873/148‑1 (MQ=255)
                                                                                                                                             cATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATc     >  2:2824714/1‑148 (MQ=255)
                                                                                                                                             cATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATc     >  2:1839145/1‑148 (MQ=255)
                                                                                                                                             cATTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATc     <  2:1398724/148‑1 (MQ=255)
                                                                                                                                              aTTAATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACCTTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCg    >  2:1256780/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTGGCTGGTGTGGTTGGTGTTGAGAAAGCCGCAAGCGCCGCAGGTTTGAGCATTCATGTACCGTTTGCGCCGGGTCGCGTTGATGCGCGTCAGGATCAGACTGACATTGAGATGTTTGAGCTGCTGGAGCCAATTGCTGACGGTTTCCGTAACTATCGCGCTCGTCTGGACGTTTCCACCACCGAGTCACTGCTGATCGAC  >  CP009273/4125262‑4125553

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: