Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,588 G→T 100% D749Y (GAC→TAC)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,5880GT97.2% 849.0 / 12.8 253D749Y (GAC→TAC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (110/136);  total (116/137)
Fisher's exact test for biased strand distribution p-value = 5.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 255 total aligned reads displayed.

GATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGC  >  CP009273/456441‑456737
                                                                                                                                                   |                                                                                                                                                     
gATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTAc                                                                                                                                                     >  1:1392370/1‑150 (MQ=255)
gATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATAGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTAc                                                                                                                                                     <  2:1689850/150‑1 (MQ=255)
    tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACg                                                                                                                                                 <  2:1871007/150‑1 (MQ=255)
    tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACg                                                                                                                                                 <  2:467664/150‑1 (MQ=255)
       cAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTaa                                                                                                                                              <  2:2675741/150‑1 (MQ=255)
        aaTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAt                                                                                                                                             >  2:2332899/1‑150 (MQ=255)
        aaTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAAt                                                                                                                                             >  2:758920/1‑150 (MQ=255)
          tGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCTCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTg                                                                                                                                           <  1:4470052/150‑1 (MQ=255)
           gACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGc                                                                                                                                          <  1:2566377/150‑1 (MQ=255)
           gACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGc                                                                                                                                          <  1:2417699/150‑1 (MQ=255)
           gACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGc                                                                                                                                          <  1:2167049/150‑1 (MQ=255)
              cGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGACGAGATTCCTTACAACGTAATTGCCGa                                                                                                                                       >  1:4311957/1‑150 (MQ=255)
               ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      <  1:4321394/150‑1 (MQ=255)
               ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      <  1:841171/150‑1 (MQ=255)
               ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      >  1:1968774/1‑150 (MQ=255)
                 tGAGATCACTCTGCGTGGTCGGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCt                                                                                                                                    >  1:2801300/1‑150 (MQ=255)
                    gaTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGa                                                                                                                                 >  1:755921/1‑150 (MQ=255)
                      tCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACa                                                                                                                               <  2:223879/150‑1 (MQ=255)
                      tCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACa                                                                                                                               >  1:4381800/1‑150 (MQ=255)
                         ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTc                                                                                                                            <  2:1597378/150‑1 (MQ=255)
                          tctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCa                                                                                                                           <  1:1900086/150‑1 (MQ=255)
                           ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                          >  2:1166898/1‑150 (MQ=255)
                              cGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCt                                                                                                                       <  1:3540143/150‑1 (MQ=255)
                              cGTGGTCAGGTACTGCCGATCGGTGGATTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCt                                                                                                                       <  2:1533308/150‑1 (MQ=255)
                                 ggTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtgtg                                                                                                                    >  2:2423332/1‑150 (MQ=255)
                                 ggTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCtgtg                                                                                                                    <  2:3267225/150‑1 (MQ=255)
                                    cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                                 <  1:2084624/150‑1 (MQ=255)
                                     aGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCAACGCGGCTGGATTAAAACAGTGCTAATTACGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTCTGAAgc                                                                                                                >  1:3091019/1‑150 (MQ=255)
                                       gTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcgc                                                                                                              >  2:1274889/1‑150 (MQ=255)
                                          cTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAtt                                                                                                           <  1:3728444/150‑1 (MQ=255)
                                           tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTg                                                                                                          <  2:3748154/150‑1 (MQ=255)
                                           tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTg                                                                                                          <  2:832081/150‑1 (MQ=255)
                                           tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAGTg                                                                                                          >  1:3508761/1‑150 (MQ=255)
                                             ccGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAg                                                                                                        <  1:1192012/150‑1 (MQ=255)
                                              cGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAgg                                                                                                       >  2:1040047/1‑150 (MQ=255)
                                               gATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGACGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGa                                                                                                      >  1:185192/1‑150 (MQ=255)
                                                aTCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGaa                                                                                                     <  1:215244/150‑1 (MQ=255)
                                                   ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGtt                                                                                                  <  1:3466794/150‑1 (MQ=255)
                                                   ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGtt                                                                                                  <  1:2754902/150‑1 (MQ=255)
                                                     tggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCt                                                                                                <  2:1353320/150‑1 (MQ=255)
                                                        tttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGAc                                                                                             <  2:2583734/150‑1 (MQ=255)
                                                        tttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGAc                                                                                             <  2:1673006/150‑1 (MQ=255)
                                                           gaaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGActct                                                                                          <  2:2890118/150‑1 (MQ=255)
                                                            aaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTg                                                                                         >  2:3642919/1‑150 (MQ=255)
                                                                 aaaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCt                                                                                    >  2:2059956/1‑150 (MQ=255)
                                                                  aaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTg                                                                                   >  1:3252640/1‑150 (MQ=255)
                                                                   aaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTgg                                                                                  <  1:2938745/150‑2 (MQ=255)
                                                                    aCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCa                                                                                 <  1:2549626/150‑1 (MQ=255)
                                                                     cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                                >  2:4322164/1‑150 (MQ=255)
                                                                     cTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaa                                                                                >  2:1785087/1‑150 (MQ=255)
                                                                       ccTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCaaaa                                                                              >  2:950369/1‑150 (MQ=255)
                                                                        cTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAAt                                                                             >  1:3607060/1‑150 (MQ=255)
                                                                           gCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGaa                                                                          <  1:1144427/150‑1 (MQ=255)
                                                                            cAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAc                                                                         >  2:1566633/1‑150 (MQ=255)
                                                                             aGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAcc                                                                        <  1:692298/150‑1 (MQ=255)
                                                                              gcgcATCGCGGCGGGATTAAAACAGTGCTACTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCg                                                                       <  1:813705/150‑1 (MQ=255)
                                                                                gcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTc                                                                     <  1:1233642/150‑1 (MQ=255)
                                                                                gcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTc                                                                     >  1:3890677/1‑150 (MQ=255)
                                                                                  aTCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTcgg                                                                   >  2:4422554/1‑148 (MQ=255)
                                                                                  aTCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTg                                                                   <  1:4217790/150‑1 (MQ=255)
                                                                                       ggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATg                                                              <  1:1037663/150‑1 (MQ=255)
                                                                                       ggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATg                                                              >  2:3891645/1‑150 (MQ=255)
                                                                                         cggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCa                                                            <  2:2771795/150‑1 (MQ=255)
                                                                                         cggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCa                                                            <  2:3431371/150‑1 (MQ=255)
                                                                                          gggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAg                                                           <  1:1240409/150‑1 (MQ=255)
                                                                                          gggATTAAAACAGTGCTAATTCCGGTCGAAAATAAACGCGCTCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGTCATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAg                                                           <  2:623552/150‑1 (MQ=255)
                                                                                           ggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAgg                                                          >  2:2643867/1‑150 (MQ=255)
                                                                                           ggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAgg                                                          >  1:4593160/1‑150 (MQ=255)
                                                                                             aTTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCCGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGtt                                                        >  1:4528682/1‑150 (MQ=255)
                                                                                                aaaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgtg                                                     <  1:579975/150‑1 (MQ=255)
                                                                                                  aaCAGTGCTAATTCCGTTCGAAAATAAACGCGAACTCGAAGATATACCTTACAAAGTAATAGCCGATCTGGACAATCATCCTGTGAAGCGCATTTAGGAAGTTCTGACTCTAGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGAc                                                   >  1:3307537/1‑150 (MQ=255)
                                                                                                    cAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTg                                                 <  1:2996133/150‑1 (MQ=255)
                                                                                                       tGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCaa                                              >  2:545300/1‑150 (MQ=255)
                                                                                                       tGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCaa                                              <  1:1794131/150‑1 (MQ=255)
                                                                                                         cTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCaaaa                                            >  2:1438521/1‑150 (MQ=255)
                                                                                                          tAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAAt                                           >  2:2536496/1‑150 (MQ=255)
                                                                                                           aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                          >  2:2722191/1‑150 (MQ=255)
                                                                                                           aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                          <  1:3483382/150‑1 (MQ=255)
                                                                                                           aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                          <  2:4219023/150‑1 (MQ=255)
                                                                                                            aTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAg                                         <  1:4062789/150‑1 (MQ=255)
                                                                                                             ttCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGt                                        <  1:3151026/150‑1 (MQ=255)
                                                                                                             ttCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCCGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGt                                        >  1:603761/1‑150 (MQ=255)
                                                                                                              tCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTg                                       >  2:3368552/1‑150 (MQ=255)
                                                                                                                cGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGAc                                     <  2:1958521/150‑1 (MQ=255)
                                                                                                                 gTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGAcc                                    <  2:3907166/150‑1 (MQ=255)
                                                                                                                  ttCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCt                                   >  1:885417/1‑150 (MQ=255)
                                                                                                                     gAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcgc                                <  1:4366564/150‑1 (MQ=255)
                                                                                                                         aTAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaa                            >  2:3999793/1‑150 (MQ=255)
                                                                                                                          tAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaa                           <  1:1343822/150‑1 (MQ=255)
                                                                                                                             aCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATg                        <  1:4139153/150‑1 (MQ=255)
                                                                                                                              cgcgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGc                       >  2:1393502/1‑150 (MQ=255)
                                                                                                                              cgcgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGc                       >  1:1064660/1‑150 (MQ=255)
                                                                                                                                cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                     >  1:1471978/1‑150 (MQ=255)
                                                                                                                                     tGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaata                >  1:3737466/1‑150 (MQ=255)
                                                                                                                                      ggAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaataa               >  1:2348614/1‑150 (MQ=255)
                                                                                                                                        aaGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaa             >  2:3848227/1‑150 (MQ=255)
                                                                                                                                         agagATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaaa            >  2:248569/1‑150 (MQ=255)
                                                                                                                                            gaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAg         <  1:399849/150‑1 (MQ=255)
                                                                                                                                            gaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAg         >  1:2447963/1‑150 (MQ=255)
                                                                                                                                                   gACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGc  >  1:2561232/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGC  >  CP009273/456441‑456737

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: