Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,298 T→A 94.5% L652Q (CTG→CAG)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,2980TA94.5% 1016.7 / 35.6 290L652Q (CTG→CAG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (3/13);  new base A (151/123);  total (154/136)
Fisher's exact test for biased strand distribution p-value = 7.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 292 total aligned reads displayed.

GGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGG  >  CP009273/456151‑456447
                                                                                                                                                   |                                                                                                                                                     
  aCGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTCTACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTGCGTGCGCGTGCGGAAAAACAgggg                                                                                                                                                   >  1:842493/1‑150 (MQ=255)
   cGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGa                                                                                                                                                  <  2:4171689/150‑1 (MQ=255)
    ggAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGAt                                                                                                                                                 <  2:1390372/150‑1 (MQ=255)
          gggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAccc                                                                                                                                           >  1:266250/1‑150 (MQ=255)
          gggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAAccc                                                                                                                                           <  1:4165428/150‑1 (MQ=255)
            gCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTg                                                                                                                                         <  1:4342803/150‑1 (MQ=255)
               gTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAtt                                                                                                                                      <  1:3569770/150‑1 (MQ=255)
                tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                     <  2:3202184/150‑1 (MQ=255)
                tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttt                                                                                                                                     <  2:1017597/150‑1 (MQ=255)
                 gACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAtttt                                                                                                                                    >  2:3961330/1‑150 (MQ=255)
                  aCTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGAttttt                                                                                                                                   >  2:2087229/1‑150 (MQ=255)
                   cTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTa                                                                                                                                  >  2:753854/1‑150 (MQ=255)
                       cTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaa                                                                                                                              >  2:773954/1‑150 (MQ=255)
                        tGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaaa                                                                                                                             >  2:372714/1‑150 (MQ=255)
                         gACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaaaa                                                                                                                            >  2:1737917/1‑150 (MQ=255)
                          aCCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGaaaaa                                                                                                                           >  1:4047374/1‑150 (MQ=255)
                           ccATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAAc                                                                                                                          <  2:683494/150‑1 (MQ=255)
                           ccATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAAc                                                                                                                          <  2:415468/150‑1 (MQ=255)
                            cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACg                                                                                                                         >  2:369321/1‑150 (MQ=255)
                            cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACg                                                                                                                         <  1:1882527/150‑1 (MQ=255)
                                gAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGAc                                                                                                                     >  1:4144762/1‑150 (MQ=255)
                                 aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACa                                                                                                                    >  2:3407343/1‑150 (MQ=255)
                                        aTGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgt                                                                                                             >  2:657635/1‑150 (MQ=255)
                                          gtgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcc                                                                                                           <  2:574949/150‑1 (MQ=255)
                                           tgtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcca                                                                                                          >  2:3014800/1‑150 (MQ=255)
                                            gtTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtccac                                                                                                         <  2:1053243/150‑1 (MQ=255)
                                                  ggTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCg                                                                                                   <  1:69732/150‑1 (MQ=255)
                                                   gTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCgg                                                                                                  >  1:3760357/1‑150 (MQ=255)
                                                   gTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCgg                                                                                                  >  1:3930805/1‑150 (MQ=255)
                                                    tAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGa                                                                                                 <  2:1663502/150‑1 (MQ=255)
                                                       aGGCAAACTGACCAATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACAAACACGTCCACGTACCGGAAgg                                                                                              >  2:1250334/1‑150 (MQ=255)
                                                        ggCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGt                                                                                             <  1:2456970/150‑1 (MQ=255)
                                                        ggCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGt                                                                                             >  2:4059473/1‑150 (MQ=255)
                                                         gCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTg                                                                                            <  2:33600/150‑1 (MQ=255)
                                                         gCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTg                                                                                            <  2:4105099/150‑1 (MQ=255)
                                                          cAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGc                                                                                           <  2:3876102/150‑1 (MQ=255)
                                                           aaaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCg                                                                                          >  2:2442680/1‑150 (MQ=255)
                                                            aaCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGa                                                                                         <  2:1459735/150‑1 (MQ=255)
                                                             aCTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGAc                                                                                        <  1:3694503/150‑1 (MQ=255)
                                                               tGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGc                                                                                      >  2:2270241/1‑150 (MQ=255)
                                                                gACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGcc                                                                                     >  2:207882/1‑150 (MQ=255)
                                                                 aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCg                                                                                    <  2:915503/150‑1 (MQ=255)
                                                                  ccTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGa                                                                                   <  2:330893/150‑1 (MQ=255)
                                                                   cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaa                                                                                  >  2:2175089/1‑150 (MQ=255)
                                                                   cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaa                                                                                  >  2:1719454/1‑150 (MQ=255)
                                                                   cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaa                                                                                  <  1:1823914/150‑1 (MQ=255)
                                                                      taCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGa                                                                               >  2:1160297/1‑150 (MQ=255)
                                                                      taCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGa                                                                               >  1:1950103/1‑150 (MQ=255)
                                                                       aCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGAt                                                                              <  2:3590892/150‑1 (MQ=255)
                                                                        ccGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATg                                                                             <  2:3486986/150‑1 (MQ=255)
                                                                        ccGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATg                                                                             <  1:3510386/150‑1 (MQ=255)
                                                                         cGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATgg                                                                            <  1:3976004/150‑1 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                           >  1:1231379/1‑150 (MQ=255)
                                                                          ggTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                           >  2:416042/1‑150 (MQ=255)
                                                                          ggTACGCTCGGCGAAGTGATGCAGTAGTCCATTCAGGCGGCGATAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGTGATCAACCCTGATTTATACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGt                                                                           >  2:416090/1‑150 (MQ=255)
                                                                           gTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                          >  2:1521770/1‑150 (MQ=255)
                                                                              cGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGa                                                                       >  1:3902281/1‑150 (MQ=255)
                                                                              cGCTCGGCGAAGTGATGCAGGAGTCCAGTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGa                                                                       <  2:2992235/150‑1 (MQ=255)
                                                                               gCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAg                                                                      <  1:3447976/150‑1 (MQ=255)
                                                                                   ggCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGcc                                                                  <  1:106584/150‑1 (MQ=255)
                                                                                      gAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGACGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGt                                                               >  2:2342142/1‑150 (MQ=255)
                                                                                      gAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGACGGt                                                               >  1:3488731/1‑150 (MQ=255)
                                                                                          tGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTg                                                           <  2:2101597/150‑1 (MQ=255)
                                                                                          tGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTg                                                           >  1:4197311/1‑150 (MQ=255)
                                                                                            aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCt                                                         <  1:298982/150‑1 (MQ=255)
                                                                                            aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCt                                                         <  1:3667689/150‑1 (MQ=255)
                                                                                             tGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTa                                                        <  2:4314637/150‑1 (MQ=255)
                                                                                             tGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTa                                                        <  2:1131070/150‑1 (MQ=255)
                                                                                               cAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtg                                                      >  1:2159399/1‑150 (MQ=255)
                                                                                                aGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtgt                                                     >  1:3215945/1‑150 (MQ=255)
                                                                                                aGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtgt                                                     >  1:1027586/1‑150 (MQ=255)
                                                                                                aGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtgt                                                     <  2:1182655/150‑1 (MQ=255)
                                                                                                   aGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCa                                                  >  2:2912682/1‑150 (MQ=255)
                                                                                                     tCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAcc                                                >  1:1444195/1‑150 (MQ=255)
                                                                                                      ccATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACcg                                               >  1:3779696/1‑150 (MQ=255)
                                                                                                       cATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACcgc                                              >  2:3128312/1‑150 (MQ=255)
                                                                                                        aTTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACcgcg                                             >  1:4036588/1‑150 (MQ=255)
                                                                                                          tCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCt                                           >  2:1938132/1‑150 (MQ=255)
                                                                                                                  cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                                   <  1:3306257/150‑1 (MQ=255)
                                                                                                                  cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                                   <  1:2072030/150‑1 (MQ=255)
                                                                                                                  cgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                                   <  2:220584/150‑1 (MQ=255)
                                                                                                                    ttAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGc                                 >  1:3543674/1‑150 (MQ=255)
                                                                                                                     tAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                                <  2:2514382/150‑1 (MQ=255)
                                                                                                                     tAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                                <  1:1130498/150‑1 (MQ=255)
                                                                                                                     tAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                                >  2:3145555/1‑150 (MQ=255)
                                                                                                                     tAACGGTGGTTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                                >  2:1893346/1‑150 (MQ=255)
                                                                                                                         ggtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                            <  1:2387755/150‑1 (MQ=255)
                                                                                                                          gtggtTCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAcc                           <  1:2232973/150‑1 (MQ=255)
                                                                                                                               tCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTaa                      <  1:3262286/150‑1 (MQ=255)
                                                                                                                               tCGTGCGCGTGCGGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTaa                      >  1:1578639/1‑150 (MQ=255)
                                                                                                                                cgtgcgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAc                     >  1:666016/1‑150 (MQ=255)
                                                                                                                                cgtgcgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAc                     >  1:1665731/1‑150 (MQ=255)
                                                                                                                                   gcgcgtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCg                  <  1:3836373/150‑1 (MQ=255)
                                                                                                                                       gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc              <  2:1469807/150‑1 (MQ=255)
                                                                                                                                       gtgcgGAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc              <  2:266250/150‑1 (MQ=255)
                                                                                                                                           ggAAAAACAGGGGATCAACCCTGATTTTTACGAAAATCGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGc          >  2:2303/1‑150 (MQ=255)
                                                                                                                                            gAAAAACAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGcc         >  2:329481/1‑150 (MQ=255)
                                                                                                                                             aaaaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCg        >  1:2258190/1‑150 (MQ=255)
                                                                                                                                             aaaaaCAGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCg        >  2:3036374/1‑150 (MQ=255)
                                                                                                                                                   tGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTgg  <  2:504577/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGG  >  CP009273/456151‑456447

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: ‑