Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,809,935 G→T 6.8% A604A (GCG→GCT katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,809,9350GT 6.8% 697.4 / 17.9 189A604A (GCG→GCTkatEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base G (94/82);  new base T (5/8);  total (99/90)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 189 total aligned reads displayed.

CAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTAC  >  CP009273/1809789‑1810081
                                                                                                                                                  |                                                                                                                                                  
cAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGa                                                                                                                                                   <  1:451716/148‑1 (MQ=255)
 aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAt                                                                                                                                                  >  2:489053/1‑148 (MQ=255)
  ggCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTCGCGAtt                                                                                                                                                 >  2:3276894/1‑148 (MQ=255)
      gTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTAc                                                                                                                                             <  2:2070114/148‑1 (MQ=255)
         gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTa                                                                                                                                          <  1:1379420/148‑1 (MQ=255)
         gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGATTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTa                                                                                                                                          <  1:2920517/148‑1 (MQ=255)
          cGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTaa                                                                                                                                         <  2:423428/148‑1 (MQ=255)
           gAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                                                                                                                        <  2:183220/148‑1 (MQ=255)
            aaaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatg                                                                                                                                       <  2:2424550/148‑1 (MQ=255)
               aaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgatg                                                                                                                                    >  2:1847425/1‑148 (MQ=255)
                    cGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                                                                                                               >  2:1437201/1‑148 (MQ=255)
                      gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgag                                                                                                                             <  2:3434000/148‑1 (MQ=255)
                      gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgag                                                                                                                             <  1:3470447/148‑1 (MQ=255)
                       tATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgaga                                                                                                                            >  1:1627820/1‑148 (MQ=255)
                            aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                       >  1:229727/1‑148 (MQ=255)
                            aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                       <  1:299355/148‑1 (MQ=255)
                            aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                       >  2:2511035/1‑148 (MQ=255)
                            aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                       <  2:3097475/148‑1 (MQ=255)
                            aaCTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                       >  1:110719/1‑148 (MQ=255)
                             actgactgACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAGAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                                                      >  2:1017146/1‑148 (MQ=255)
                             actgactgACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                                                      >  1:1030126/1‑148 (MQ=255)
                             actgactgACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                                                      >  2:2617598/1‑148 (MQ=255)
                             actgactgACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                                                      >  1:3423516/1‑148 (MQ=255)
                              ctgactgaCGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAg                                                                                                                     >  1:1926868/1‑148 (MQ=255)
                              ctgactgaCGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAg                                                                                                                     >  1:1402536/1‑148 (MQ=255)
                              ctgactgaCGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAg                                                                                                                     >  2:485362/1‑148 (MQ=255)
                                 actgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                                                                                                                  >  2:1486806/1‑148 (MQ=255)
                                  ctgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCt                                                                                                                 <  2:1001943/148‑1 (MQ=255)
                                  ctgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCt                                                                                                                 <  2:977717/148‑1 (MQ=255)
                                   tgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCtt                                                                                                                <  2:14926/148‑1 (MQ=255)
                                     acgacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                                                              <  2:3348052/148‑1 (MQ=255)
                                        acCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGc                                                                                                           >  2:3022327/1‑148 (MQ=255)
                                        acCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGc                                                                                                           >  1:1259022/1‑148 (MQ=255)
                                             cTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtc                                                                                                      >  1:2927662/1‑148 (MQ=255)
                                               gAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtctc                                                                                                    >  2:3215776/1‑148 (MQ=255)
                                                aaTATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCa                                                                                                   >  1:941219/1‑148 (MQ=255)
                                                 atatCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCaa                                                                                                  <  2:1835073/148‑1 (MQ=255)
                                                  tatCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAg                                                                                                 >  2:2964905/1‑148 (MQ=255)
                                                   atCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAgg                                                                                                <  1:424508/148‑1 (MQ=255)
                                                    tCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgc                                                                                               <  1:881769/148‑1 (MQ=255)
                                                    tCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTGAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgc                                                                                               >  1:1622814/1‑148 (MQ=255)
                                                     cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                                                              >  2:1227047/1‑148 (MQ=255)
                                                     cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                                                              <  1:640112/148‑1 (MQ=255)
                                                     cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                                                              <  1:3478515/148‑1 (MQ=255)
                                                     cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCGTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                                                              >  2:916395/1‑148 (MQ=255)
                                                        cccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg                                                                                           >  1:3135337/1‑148 (MQ=255)
                                                         ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGa                                                                                          <  2:2895902/148‑1 (MQ=255)
                                                         ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGa                                                                                          <  2:1121316/148‑1 (MQ=255)
                                                             cTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGc                                                                                      <  1:3352159/148‑1 (MQ=255)
                                                              tCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGcc                                                                                     <  1:1878586/148‑1 (MQ=255)
                                                                  gACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAg                                                                                 >  1:835434/1‑148 (MQ=255)
                                                                    cGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGc                                                                               <  2:886159/148‑1 (MQ=255)
                                                                        aaCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                                                           <  2:2642761/148‑1 (MQ=255)
                                                                        aaCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                                                           <  2:809386/148‑1 (MQ=255)
                                                                          cGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCAt                                                                         >  1:1666051/1‑148 (MQ=255)
                                                                           ggTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATg                                                                        <  1:2405150/148‑1 (MQ=255)
                                                                           ggTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATg                                                                        <  2:1796149/148‑1 (MQ=255)
                                                                                 aaaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAAc                                                                  <  2:228173/148‑1 (MQ=255)
                                                                                  aaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                                                                 >  1:400573/1‑148 (MQ=255)
                                                                                  aaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                                                                 <  2:2214397/148‑1 (MQ=255)
                                                                                     aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGct                                                              <  1:2471892/148‑1 (MQ=255)
                                                                                      ggATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctc                                                             >  1:1470401/1‑148 (MQ=255)
                                                                                      ggATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctc                                                             >  1:3235696/1‑148 (MQ=255)
                                                                                      ggATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctc                                                             <  2:1655312/148‑1 (MQ=255)
                                                                                       gATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                                                            <  1:2508587/148‑1 (MQ=255)
                                                                                        atccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTa                                                           >  2:1247761/1‑148 (MQ=255)
                                                                                        atccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTa                                                           >  1:1951016/1‑148 (MQ=255)
                                                                                             tccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                                                      <  1:3435103/148‑1 (MQ=255)
                                                                                             tccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                                                      >  2:3374776/1‑148 (MQ=255)
                                                                                                ttAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGaa                                                   <  2:1126577/148‑1 (MQ=255)
                                                                                                 tAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAAt                                                  <  2:2842956/148‑1 (MQ=255)
                                                                                                     tttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGt                                              >  2:3403404/1‑148 (MQ=255)
                                                                                                      ttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTg                                             >  2:223019/1‑148 (MQ=255)
                                                                                                        gTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGaa                                           >  1:2872088/1‑148 (MQ=255)
                                                                                                           cGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTg                                        >  2:2997014/1‑148 (MQ=255)
                                                                                                            gCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGa                                       <  2:704673/148‑1 (MQ=255)
                                                                                                             ccATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGAc                                      <  2:1174066/148‑1 (MQ=255)
                                                                                                              cATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACt                                     >  1:842611/1‑148 (MQ=255)
                                                                                                              cATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACt                                     >  1:324147/1‑148 (MQ=255)
                                                                                                               aTTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTg                                    <  2:564815/148‑1 (MQ=255)
                                                                                                                taccTGACGGTGATGTGAAAGGTCGCGTCGTAGCGATTTTACTTAAAGATGAAGTGAGATCGGCAGACCTACTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAAAGGGTGAAGTGACTGc                                   >  2:364344/3‑148 (MQ=255)
                                                                                                                 tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                  >  1:174129/1‑148 (MQ=255)
                                                                                                                 tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCg                                  >  1:283256/1‑148 (MQ=255)
                                                                                                                     gACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATg                              >  1:2147967/1‑148 (MQ=255)
                                                                                                                      aCGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGa                             >  1:3125886/1‑148 (MQ=255)
                                                                                                                      aCGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGa                             >  1:3005463/1‑148 (MQ=255)
                                                                                                                      aCGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGa                             <  1:499023/148‑1 (MQ=255)
                                                                                                                          tGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggt                         >  2:2836458/1‑148 (MQ=255)
                                                                                                                              gtgAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggtacgg                     <  1:2259784/148‑1 (MQ=255)
                                                                                                                                 aaaGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGgtgt                  >  1:1332394/1‑148 (MQ=255)
                                                                                                                                    ggTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGc               <  1:2034050/148‑1 (MQ=255)
                                                                                                                                    ggTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGc               <  2:659244/148‑1 (MQ=255)
                                                                                                                                    ggTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGc               <  2:2640638/148‑1 (MQ=255)
                                                                                                                                     gTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGcc              <  2:275648/148‑1 (MQ=255)
                                                                                                                                     gTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGcc              <  2:1718280/148‑1 (MQ=255)
                                                                                                                                      tCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCt             >  1:2405382/1‑148 (MQ=255)
                                                                                                                                         cgTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCtata          >  1:1188134/1‑148 (MQ=255)
                                                                                                                                          gTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAg         <  1:719088/148‑1 (MQ=255)
                                                                                                                                           tGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGc        >  1:133515/1‑148 (MQ=255)
                                                                                                                                               aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCt    <  2:1402536/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTAC  >  CP009273/1809789‑1810081

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: