Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 421,799 C→A 88.6% I69I (ATC→ATA tgt → tRNA‑guanine transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273421,7990CA88.6% 792.7 / 99.7 290I69I (ATC→ATAtgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (13/20);  new base A (126/130);  total (139/151)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 295 total aligned reads displayed.

GATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGA  >  CP009273/421650‑421945
                                                                                                                                                     |                                                                                                                                                  
   cGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATg                                                                                                                                                 <  1:1825162/150‑1 (MQ=255)
        cGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACt                                                                                                                                            <  2:2220000/150‑1 (MQ=255)
         gTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTg                                                                                                                                           >  1:316780/1‑150 (MQ=255)
           aGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCTACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCa                                                                                                                                         >  1:5138501/1‑150 (MQ=255)
           aGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCa                                                                                                                                         <  1:1973456/150‑1 (MQ=255)
            gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCAc                                                                                                                                        <  2:3051085/150‑1 (MQ=255)
             tGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACg                                                                                                                                       >  1:2177477/1‑150 (MQ=255)
                aaaCGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCg                                                                                                                                    <  2:3703936/150‑1 (MQ=255)
                   cGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATc                                                                                                                                 <  2:2445325/150‑1 (MQ=255)
                    gCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCt                                                                                                                                <  2:1515302/150‑1 (MQ=255)
                       ttGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCa                                                                                                                             <  2:3122216/150‑1 (MQ=255)
                        tGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCAc                                                                                                                            <  2:2810573/150‑1 (MQ=255)
                             tATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGAttt                                                                                                                       >  2:3403818/1‑150 (MQ=255)
                              aTGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGAtttt                                                                                                                      >  2:681847/1‑150 (MQ=255)
                              aTGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGAtttt                                                                                                                      >  1:2475867/1‑150 (MQ=255)
                                    gTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAg                                                                                                                <  2:416943/150‑1 (MQ=255)
                                     ttGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGt                                                                                                               >  2:4018311/1‑150 (MQ=255)
                                       ggCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTgg                                                                                                             <  2:2285684/150‑1 (MQ=255)
                                         cACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGaa                                                                                                           >  2:3889424/1‑150 (MQ=255)
                                          aCCTACGGCACCGTAAAAGGTATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAg                                                                                                          >  1:472789/1‑150 (MQ=255)
                                          aCCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAg                                                                                                          >  2:4083542/1‑150 (MQ=255)
                                            cTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAggg                                                                                                        <  1:2100418/150‑1 (MQ=255)
                                            cTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAggg                                                                                                        <  1:662419/150‑1 (MQ=255)
                                             tACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgggg                                                                                                       >  1:2042252/1‑150 (MQ=255)
                                             tACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgggg                                                                                                       >  1:2981055/1‑150 (MQ=255)
                                             tACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgggg                                                                                                       >  1:4191390/1‑150 (MQ=255)
                                              aCGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGc                                                                                                      <  1:2256944/150‑1 (MQ=255)
                                                ggCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCg                                                                                                    >  2:3569926/1‑150 (MQ=255)
                                                  cACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGAt                                                                                                  <  1:61109/150‑1 (MQ=255)
                                                    ccGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATcc                                                                                                <  2:552388/150‑1 (MQ=255)
                                                     cGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCt                                                                                               >  1:138219/1‑150 (MQ=255)
                                                        aaaaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCAc                                                                                            <  2:4853597/150‑1 (MQ=255)
                                                         aaaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAACCTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCAcc                                                                                           <  2:383577/150‑1 (MQ=255)
                                                            gggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGAc                                                                                        <  1:3880234/150‑1 (MQ=255)
                                                            gggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGAc                                                                                        <  1:2136680/150‑1 (MQ=255)
                                                            gggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGAc                                                                                        <  2:1526555/150‑1 (MQ=255)
                                                             ggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACt                                                                                       >  1:3315549/1‑150 (MQ=255)
                                                             ggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACATGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACt                                                                                       >  1:727365/1‑150 (MQ=255)
                                                                      cGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCt                                                                              <  2:2535008/150‑1 (MQ=255)
                                                                       ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                             >  1:4045600/1‑150 (MQ=255)
                                                                       ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                             >  2:3453066/1‑150 (MQ=255)
                                                                       ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                             >  2:3452774/1‑150 (MQ=255)
                                                                         aagaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTcc                                                                           <  1:121784/150‑1 (MQ=255)
                                                                          agaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCa                                                                          <  1:4708711/150‑1 (MQ=255)
                                                                          agaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCa                                                                          <  2:636828/150‑1 (MQ=255)
                                                                          agaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCa                                                                          <  2:2098478/150‑1 (MQ=255)
                                                                            aagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAgg                                                                        >  2:4594372/1‑150 (MQ=255)
                                                                            aagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAgg                                                                        >  2:2039864/1‑150 (MQ=255)
                                                                            aagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAgg                                                                        >  2:1283435/1‑150 (MQ=255)
                                                                              gTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTc                                                                      <  1:4967638/150‑1 (MQ=255)
                                                                                tGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCtt                                                                    <  1:4334457/150‑1 (MQ=255)
                                                                                 gAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTc                                                                   <  1:2309835/150‑1 (MQ=255)
                                                                                  aaGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCa                                                                  <  2:4424276/150‑1 (MQ=255)
                                                                                  aaGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCa                                                                  <  2:519834/150‑1 (MQ=255)
                                                                                    gCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGc                                                                >  2:404028/1‑150 (MQ=255)
                                                                                      cACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCt                                                              >  1:1593537/1‑150 (MQ=255)
                                                                                        cTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTg                                                            >  1:2427965/1‑150 (MQ=255)
                                                                                          ggCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGc                                                          <  2:1093487/150‑1 (MQ=255)
                                                                                            cgcgcAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGa                                                        <  2:2610461/150‑1 (MQ=255)
                                                                                                  aaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCg                                                  >  2:2159414/1‑150 (MQ=255)
                                                                                                  aaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCTATATTCg                                                  >  1:2587923/1‑150 (MQ=255)
                                                                                                   aTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGt                                                 >  2:5014476/1‑150 (MQ=255)
                                                                                                     tATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaa                                               >  2:4291960/1‑150 (MQ=255)
                                                                                                     tATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaa                                               >  1:4546209/1‑150 (MQ=255)
                                                                                                        ccTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAAt                                            <  1:2670676/150‑1 (MQ=255)
                                                                                                        ccTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAAt                                            <  1:2720382/150‑1 (MQ=255)
                                                                                                         cTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATc                                           <  2:4306955/150‑1 (MQ=255)
                                                                                                          tCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCa                                          <  2:3892300/150‑1 (MQ=255)
                                                                                                           cGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCAc                                         >  1:1925272/1‑150 (MQ=255)
                                                                                                             gCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCg                                       >  1:2889056/1‑150 (MQ=255)
                                                                                                                acacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAAc                                    <  1:1833005/150‑1 (MQ=255)
                                                                                                                acacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAAc                                    <  1:1283435/150‑1 (MQ=255)
                                                                                                                 cacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACa                                   >  2:1432332/1‑150 (MQ=255)
                                                                                                                 cacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACa                                   >  2:4486673/1‑150 (MQ=255)
                                                                                                                      tCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGACGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCg                              >  1:1053015/1‑150 (MQ=255)
                                                                                                                      tCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCg                              >  2:1595485/1‑150 (MQ=255)
                                                                                                                      tCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCg                              >  2:3772366/1‑150 (MQ=255)
                                                                                                                        cACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                            >  1:1020052/1‑150 (MQ=255)
                                                                                                                         aCCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGc                           <  1:4923522/150‑1 (MQ=255)
                                                                                                                         aCCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGc                           >  1:2811415/1‑150 (MQ=255)
                                                                                                                          ccTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCa                          >  1:861744/1‑150 (MQ=255)
                                                                                                                            tgtgGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACt                        >  2:831517/1‑150 (MQ=255)
                                                                                                                                 cTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGt                   >  1:3392729/1‑150 (MQ=255)
                                                                                                                                  tGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTa                  >  2:4452391/1‑150 (MQ=255)
                                                                                                                                    cgcCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAc                >  2:3269109/1‑150 (MQ=255)
                                                                                                                                      cccGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAccc              >  1:4829054/1‑150 (MQ=255)
                                                                                                                                       ccGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCg             <  1:3087383/150‑1 (MQ=255)
                                                                                                                                        cGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGa            <  2:3552297/150‑1 (MQ=255)
                                                                                                                                         gggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt           >  2:1031675/1‑150 (MQ=255)
                                                                                                                                          ggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATc          <  1:1004696/150‑1 (MQ=255)
                                                                                                                                           gCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCa         >  1:3810577/1‑150 (MQ=255)
                                                                                                                                            ccAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCaa        >  1:4846563/1‑150 (MQ=255)
                                                                                                                                              aGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGATGGGGCCGATCCTCACCGACTCCGGCGTCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCTCTTCCGTAACCCGATCAACg      <  1:96376/150‑1 (MQ=255)
                                                                                                                                              aGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACg      <  2:3175169/150‑1 (MQ=255)
                                                                                                                                              aGGAAATAATGAAACGGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACg      <  1:889067/150‑1 (MQ=255)
                                                                                                                                               ggAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACgg     >  1:3480104/1‑150 (MQ=255)
                                                                                                                                                gAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc    <  2:4527071/150‑1 (MQ=255)
                                                                                                                                                gAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc    <  1:4637796/150‑1 (MQ=255)
                                                                                                                                                 aaaTAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCg   >  2:1697204/1‑150 (MQ=255)
                                                                                                                                                  aataatGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGa  <  2:2475867/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGA  >  CP009273/421650‑421945

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: