Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,959,351 G→T 90.3% P77P (CCG→CCT trxA → thioredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,959,3510GT90.3% 788.2 / 79.2 288P77P (CCG→CCTtrxAthioredoxin 1
Reads supporting (aligned to +/- strand):  ref base G (16/12);  new base T (135/125);  total (151/137)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 288 total aligned reads displayed.

TGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACGCCC  >  CP009273/3959205‑3959495
                                                                                                                                                  |                                                                                                                                                
     agagTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgc                                                                                                                                           <  2:4254316/149‑1 (MQ=255)
      gagTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgct                                                                                                                                          >  1:3906063/1‑149 (MQ=255)
      gagTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgct                                                                                                                                          <  2:1648698/149‑1 (MQ=255)
        gtggtgCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgctgc                                                                                                                                        <  2:3020553/149‑1 (MQ=255)
         tggtgCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgctgct                                                                                                                                       <  2:918828/149‑1 (MQ=255)
          ggtgCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTctgctgctg                                                                                                                                      >  1:4078562/1‑149 (MQ=255)
           gtgCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGt                                                                                                                                     <  2:3122975/149‑1 (MQ=255)
                  ccGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCaaaaa                                                                                                                              <  2:2529525/149‑1 (MQ=255)
                        aaaaTGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGa                                                                                                                        >  2:4234030/1‑149 (MQ=255)
                         aaaTGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGaa                                                                                                                       <  2:3701703/149‑1 (MQ=255)
                          aaTGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAg                                                                                                                      >  2:4658556/1‑149 (MQ=255)
                             gATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTgg                                                                                                                   <  2:3549129/149‑1 (MQ=255)
                              aTCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTggc                                                                                                                  <  1:104924/149‑1 (MQ=255)
                                cGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTggcgg                                                                                                                <  2:4085012/149‑1 (MQ=255)
                                  ccccGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCa                                                                                                              >  1:4272464/1‑149 (MQ=255)
                                   cccGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCaa                                                                                                             <  1:1516356/149‑1 (MQ=255)
                                      gATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCa                                                                                                          >  1:889756/1‑149 (MQ=255)
                                        ttCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCaaa                                                                                                        >  1:667358/1‑149 (MQ=255)
                                             gATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTggg                                                                                                   >  1:2486358/1‑149 (MQ=255)
                                             gATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTggg                                                                                                   <  1:4617615/149‑1 (MQ=255)
                                             gATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTggg                                                                                                   >  2:1913216/1‑149 (MQ=255)
                                                gAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGc                                                                                                >  2:1628679/1‑149 (MQ=255)
                                                 aaaTCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCa                                                                                               >  2:3150751/1‑149 (MQ=255)
                                                 aaaTCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCa                                                                                               >  1:3833619/1‑149 (MQ=255)
                                                  aaTCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCAc                                                                                              >  2:945801/1‑149 (MQ=255)
                                                    tCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTg                                                                                            <  1:2438826/149‑1 (MQ=255)
                                                    tCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTg                                                                                            <  1:4780466/149‑1 (MQ=255)
                                                     cGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGt                                                                                           <  2:3298275/149‑1 (MQ=255)
                                                      gCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTc                                                                                          >  1:4149583/1‑149 (MQ=255)
                                                       cTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCt                                                                                         >  2:3900992/1‑149 (MQ=255)
                                                        tGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTa                                                                                        >  2:2871285/1‑149 (MQ=255)
                                                        tGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTa                                                                                        >  2:1120306/1‑149 (MQ=255)
                                                         gACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTaa                                                                                       >  2:3289155/1‑149 (MQ=255)
                                                          aCGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTaaa                                                                                      >  2:2367355/1‑149 (MQ=255)
                                                           cGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAg                                                                                     >  2:4395648/1‑149 (MQ=255)
                                                            gAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAgg                                                                                    >  1:2809797/1‑149 (MQ=255)
                                                            gAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAgg                                                                                    >  1:2369724/1‑149 (MQ=255)
                                                             aaTATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGt                                                                                   >  2:1022038/1‑149 (MQ=255)
                                                             aaTATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGt                                                                                   <  1:4672255/149‑1 (MQ=255)
                                                              atatCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTc                                                                                  >  1:1069064/1‑149 (MQ=255)
                                                              atatCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTc                                                                                  >  2:179567/1‑149 (MQ=255)
                                                              atatCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTc                                                                                  <  1:2826582/149‑1 (MQ=255)
                                                                 tCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAATGTCAGt                                                                               <  2:647712/149‑1 (MQ=255)
                                                                    gggCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGa                                                                            >  2:2294043/1‑149 (MQ=255)
                                                                       cAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAag                                                                         <  1:2171676/149‑1 (MQ=255)
                                                                       cAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAag                                                                         >  2:3193566/1‑149 (MQ=255)
                                                                        aaaCTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAaga                                                                        <  1:4442045/149‑1 (MQ=255)
                                                                        aaaCTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAaga                                                                        <  1:732173/149‑1 (MQ=255)
                                                                         aaCTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAagag                                                                       <  1:1797998/149‑1 (MQ=255)
                                                                         aaCTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAagag                                                                       >  2:844731/1‑149 (MQ=255)
                                                                            tGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTc                                                                    >  1:325203/1‑149 (MQ=255)
                                                                             gACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGACGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTcc                                                                   >  2:54394/1‑149 (MQ=255)
                                                                             gACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTcc                                                                   <  2:801600/149‑1 (MQ=255)
                                                                              aCCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCt                                                                  <  1:3796729/149‑1 (MQ=255)
                                                                              aCCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCt                                                                  <  1:4713963/149‑1 (MQ=255)
                                                                              aCCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCt                                                                  <  1:1379682/149‑1 (MQ=255)
                                                                               ccGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAGTATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTc                                                                 >  1:195874/1‑149 (MQ=255)
                                                                               ccGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTc                                                                 >  1:2589126/1‑149 (MQ=255)
                                                                               ccGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTc                                                                 >  2:625573/1‑149 (MQ=255)
                                                                                 gTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACCCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGa                                                               >  1:787175/1‑149 (MQ=255)
                                                                                    gCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGc                                                            <  2:2404992/149‑1 (MQ=255)
                                                                                    gCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGc                                                            >  1:828309/1‑149 (MQ=255)
                                                                                     cAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCt                                                           >  2:1047946/1‑149 (MQ=255)
                                                                                       aaaCTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTTAAAGAGTTCCTCGACGCTaa                                                         >  2:2920999/1‑149 (MQ=255)
                                                                                       aaaCTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTaa                                                         >  1:2721922/1‑149 (MQ=255)
                                                                                         aCTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAAcc                                                       <  1:2529132/149‑1 (MQ=255)
                                                                                         aCTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAAcc                                                       >  2:70761/1‑149 (MQ=255)
                                                                                          cTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCt                                                      <  1:3370913/149‑1 (MQ=255)
                                                                                           tGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTg                                                     >  1:1619094/1‑149 (MQ=255)
                                                                                            gAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTgg                                                    >  1:130101/1‑149 (MQ=255)
                                                                                            gAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTgg                                                    <  2:3269123/149‑1 (MQ=255)
                                                                                             aaCATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGc                                                   <  1:3442168/149‑1 (MQ=255)
                                                                                              aCATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCg                                                  >  2:1767209/1‑149 (MQ=255)
                                                                                                  cGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAg                                              <  2:751282/149‑1 (MQ=255)
                                                                                                   gATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAgg                                             >  1:167188/1‑149 (MQ=255)
                                                                                                     tCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGa                                           <  2:1241198/149‑1 (MQ=255)
                                                                                                        aaaCCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAAtt                                        <  1:2521296/149‑1 (MQ=255)
                                                                                                        aaaCCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTACAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAAtt                                        >  2:2584012/1‑149 (MQ=255)
                                                                                                        aaaCCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAAtt                                        >  1:4160595/1‑149 (MQ=255)
                                                                                                         caCCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAAttt                                       <  1:1502533/148‑1 (MQ=255)
                                                                                                             cTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATg                                   >  1:4413906/1‑149 (MQ=255)
                                                                                                             cTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATg                                   >  2:1815396/1‑149 (MQ=255)
                                                                                                              tGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGt                                  >  2:3065180/1‑149 (MQ=255)
                                                                                                              tGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGt                                  <  2:4789482/149‑1 (MQ=255)
                                                                                                                    tGCGCCGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCGTGTTCGGGt                            <  2:1375905/149‑1 (MQ=255)
                                                                                                                        ccGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGccc                        <  1:1903711/149‑1 (MQ=255)
                                                                                                                         cGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGcccc                       <  1:1670169/149‑1 (MQ=255)
                                                                                                                         cGAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGcccc                       <  1:296518/149‑1 (MQ=255)
                                                                                                                          gAAATATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCg                      >  1:2128970/1‑149 (MQ=255)
                                                                                                                            aaTATGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTc                    >  2:2176869/1‑149 (MQ=255)
                                                                                                                            aaTATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTc                    <  1:2762146/149‑1 (MQ=255)
                                                                                                                             atatGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCg                   >  1:1431618/1‑149 (MQ=255)
                                                                                                                               atGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCt                 >  1:2643842/1‑149 (MQ=255)
                                                                                                                               atGGCATCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCt                 >  1:2643782/1‑149 (MQ=255)
                                                                                                                                    aTCCGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTaaaaa            >  1:4302160/1‑149 (MQ=255)
                                                                                                                                      ccGTGGTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACt          >  2:1397580/1‑149 (MQ=255)
                                                                                                                                           gTATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACg     >  1:529277/1‑149 (MQ=255)
                                                                                                                                            tATCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACGc    <  2:630169/149‑1 (MQ=255)
                                                                                                                                             aTCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACGcc   >  1:3683925/1‑149 (MQ=255)
                                                                                                                                              tCCCTACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACGccc  >  2:3640317/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAACTGACCGTTGCAAAACTGAACATCGATCAAAACCCTGGCACTGCGCCGAAATATGGCATCCGTGGTATCCCGACTCTGCTGCTGTTCAAAAACGGTGAAGTGGCGGCAACCAAAGTGGGTGCACTGTCTAAAGGTCAGTTGAAAGAGTTCCTCGACGCTAACCTGGCGTAAGGGAATTTCATGTTCGGGTGCCCCGTCGCTAAAAACTGGACGCCC  >  CP009273/3959205‑3959495

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: