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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273421,7990CA29.2% 446.4 / 189.6 229I69I (ATC→ATAtgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (85/77);  new base A (33/34);  total (118/111)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

CGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCC  >  CP009273/421653‑421948
                                                                                                                                                  |                                                                                                                                                     
cGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATg                                                                                                                                                    <  1:4438903/150‑1 (MQ=255)
     cGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACt                                                                                                                                               <  2:1580912/150‑1 (MQ=255)
        aGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCa                                                                                                                                            <  2:3675474/150‑1 (MQ=255)
         gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCAc                                                                                                                                           >  1:3817484/1‑150 (MQ=255)
          tGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACg                                                                                                                                          <  1:1687069/150‑1 (MQ=255)
          tGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACg                                                                                                                                          <  1:4532723/150‑1 (MQ=255)
            gAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGc                                                                                                                                        <  2:370339/150‑1 (MQ=255)
            gAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGc                                                                                                                                        >  1:1360454/1‑150 (MQ=255)
            gAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGc                                                                                                                                        <  2:4897382/150‑1 (MQ=255)
            gAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGc                                                                                                                                        <  1:1319430/150‑1 (MQ=255)
                 gCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGACATCATGAAACTGCACGGCGATCt                                                                                                                                   >  1:2804834/1‑150 (MQ=255)
                  ccTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTg                                                                                                                                  <  1:3230364/150‑1 (MQ=255)
                  ccTTGTTTTATGCCTGTTGGCAACTACGGCACCGTAAAATGGATGACGCCGGAAGAATTTGAAGCCACTGGCGCTCAAATTATCCTCGGCAACAACTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTg                                                                                                                                  <  2:464293/150‑1 (MQ=255)
                   cTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGc                                                                                                                                 <  2:856259/150‑1 (MQ=255)
                    ttGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCa                                                                                                                                >  2:4543441/1‑150 (MQ=255)
                       ttttATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGa                                                                                                                             >  1:4795467/1‑150 (MQ=255)
                           aTGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGAtttt                                                                                                                         >  1:3029825/1‑150 (MQ=255)
                             gCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTAt                                                                                                                       <  2:4466396/150‑1 (MQ=255)
                              ccTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATg                                                                                                                      <  1:3167482/150‑1 (MQ=255)
                               cTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGc                                                                                                                     >  2:5119291/1‑150 (MQ=255)
                                 gTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAg                                                                                                                   >  2:3066463/1‑150 (MQ=255)
                                 gTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAg                                                                                                                   >  1:825522/1‑150 (MQ=255)
                                   tGGCACCTACGGCACCGTAAATGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTg                                                                                                                 <  1:4865861/150‑1 (MQ=255)
                                    ggCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTgg                                                                                                                <  2:3704348/150‑1 (MQ=255)
                                      cACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGaa                                                                                                              <  2:2640223/150‑1 (MQ=255)
                                        ccTCCGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgg                                                                                                            <  2:1412556/150‑1 (MQ=255)
                                        ccTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgg                                                                                                            >  2:816185/1‑150 (MQ=255)
                                        ccTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgg                                                                                                            >  2:4204299/1‑150 (MQ=255)
                                             ggCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCg                                                                                                       >  1:2088186/1‑150 (MQ=255)
                                                aCCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATc                                                                                                    >  1:1747705/1‑150 (MQ=255)
                                                 ccGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGTAGGGGCCGATcc                                                                                                   <  2:3050781/150‑1 (MQ=255)
                                                   gTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTc                                                                                                 >  2:714353/1‑150 (MQ=255)
                                                      aaaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCAcc                                                                                              <  1:3566314/150‑1 (MQ=255)
                                                       aaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCg                                                                                             >  1:1456120/1‑150 (MQ=255)
                                                       aaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCg                                                                                             <  2:808819/150‑1 (MQ=255)
                                                          ggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACt                                                                                          <  1:4906390/150‑1 (MQ=255)
                                                              gACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcgg                                                                                      <  2:1134529/150‑1 (MQ=255)
                                                                cGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcggcg                                                                                    <  2:3225317/150‑1 (MQ=255)
                                                                cGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcggcg                                                                                    <  2:1689867/150‑1 (MQ=255)
                                                                cGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcggcg                                                                                    >  1:3767068/1‑150 (MQ=255)
                                                                 gCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcggcgg                                                                                   <  1:1808823/150‑1 (MQ=255)
                                                                  ccGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCggcggc                                                                                  >  2:5167889/1‑150 (MQ=255)
                                                                    ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                                >  2:2316807/1‑150 (MQ=255)
                                                                    ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                                >  1:1103829/1‑150 (MQ=255)
                                                                       agaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCa                                                                             >  2:1661743/1‑150 (MQ=255)
                                                                        gaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAg                                                                            <  1:3944649/150‑1 (MQ=255)
                                                                           gTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTc                                                                         <  1:1175394/150‑1 (MQ=255)
                                                                            ttGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCt                                                                        <  2:4950999/150‑1 (MQ=255)
                                                                               aaGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCa                                                                     <  2:335025/150‑1 (MQ=255)
                                                                                  ccACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGcc                                                                  <  2:4062728/150‑1 (MQ=255)
                                                                                  ccACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGcc                                                                  >  1:5021747/1‑150 (MQ=255)
                                                                                   cACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCt                                                                 >  2:3989402/1‑150 (MQ=255)
                                                                                      tGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTgg                                                              >  2:2768832/1‑150 (MQ=255)
                                                                                             cAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATAtt                                                       >  1:1313831/1‑150 (MQ=255)
                                                                                              aaaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTc                                                      >  2:847910/1‑150 (MQ=255)
                                                                                              aaaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTc                                                      >  2:1288288/1‑150 (MQ=255)
                                                                                                  tATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaa                                                  <  1:2503407/150‑1 (MQ=255)
                                                                                                   aTCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaaa                                                 >  1:2257771/1‑150 (MQ=255)
                                                                                                      cTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACGGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATCTTCGTAAAATc                                              >  2:123536/1‑150 (MQ=255)
                                                                                                         ggCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCAcc                                           <  1:2167394/150‑1 (MQ=255)
                                                                                                          gCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCg                                          <  2:4559981/150‑1 (MQ=255)
                                                                                                           cAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGa                                         >  1:992799/1‑150 (MQ=255)
                                                                                                            aaCACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGaa                                        >  1:1415793/1‑150 (MQ=255)
                                                                                                             acacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAAc                                       >  1:2337678/1‑150 (MQ=255)
                                                                                                              cacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACa                                      >  2:4178716/1‑150 (MQ=255)
                                                                                                              cacCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACa                                      <  1:3091480/150‑1 (MQ=255)
                                                                                                                ccTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAgg                                    >  2:4073658/1‑150 (MQ=255)
                                                                                                                    ccACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCTCGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGt                                <  1:4226695/150‑1 (MQ=255)
                                                                                                                    ccACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGt                                >  2:3307955/1‑150 (MQ=255)
                                                                                                                     cACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                               <  1:4948376/150‑1 (MQ=255)
                                                                                                                     cACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                               <  2:1170025/150‑1 (MQ=255)
                                                                                                                     cACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                               <  2:4201933/150‑1 (MQ=255)
                                                                                                                       ccTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCa                             >  2:4371225/1‑150 (MQ=255)
                                                                                                                        cTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCAc                            >  1:2942763/1‑150 (MQ=255)
                                                                                                                        cTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCAc                            >  1:2813610/1‑150 (MQ=255)
                                                                                                                        cTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCAc                            >  1:3598633/1‑150 (MQ=255)
                                                                                                                         tgtgGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACt                           >  1:4173807/1‑150 (MQ=255)
                                                                                                                         tgtgGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACt                           >  1:433035/1‑150 (MQ=255)
                                                                                                                           tgGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTc                         >  2:3720637/1‑150 (MQ=255)
                                                                                                                            ggCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTcc                        <  1:1554031/150‑1 (MQ=255)
                                                                                                                            ggCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTcc                        <  2:4694503/150‑1 (MQ=255)
                                                                                                                             gCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCg                       >  2:2608205/1‑150 (MQ=255)
                                                                                                                              cTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGt                      >  2:1391655/1‑150 (MQ=255)
                                                                                                                               tGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTa                     <  1:4689820/150‑1 (MQ=255)
                                                                                                                                  gcCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAcc                  <  2:3925919/150‑1 (MQ=255)
                                                                                                                                    ccGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCg                <  1:552830/150‑1 (MQ=255)
                                                                                                                                     cGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGa               >  1:4303122/1‑150 (MQ=255)
                                                                                                                                        gCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCa            <  1:2721023/150‑1 (MQ=255)
                                                                                                                                        gCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCa            <  1:2586859/150‑1 (MQ=255)
                                                                                                                                          cAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAAc          >  1:1163311/1‑150 (MQ=255)
                                                                                                                                           aGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACg         >  2:1954554/1‑150 (MQ=255)
                                                                                                                                           aGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACg         >  1:1219716/1‑150 (MQ=255)
                                                                                                                                           aGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCAGTAACCCGATCAACg         >  2:1954656/1‑150 (MQ=255)
                                                                                                                                             gAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc       >  2:2353021/1‑150 (MQ=255)
                                                                                                                                               aataatGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGa     <  2:2098673/150‑1 (MQ=255)
                                                                                                                                                aTCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGAt    <  2:5064774/150‑1 (MQ=255)
                                                                                                                                                 tCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATc   >  2:395815/1‑150 (MQ=255)
                                                                                                                                                 tCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATc   >  2:387830/1‑150 (MQ=255)
                                                                                                                                                  cATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATcc  >  1:708706/1‑150 (MQ=255)
                                                                                                                                                  cATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATcc  >  2:1808128/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCC  >  CP009273/421653‑421948

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: