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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273421,7990CA91.8% 744.4 / 57.5 265I69I (ATC→ATAtgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base A (111/131);  total (124/141)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 265 total aligned reads displayed.

GTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCG  >  CP009273/421654‑421944
                                                                                                                                                 |                                                                                                                                                 
gTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGa                                                                                                                                               >  2:1837028/1‑150 (MQ=255)
gTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGa                                                                                                                                               >  1:1431619/1‑150 (MQ=255)
 tGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGaa                                                                                                                                              <  1:4308400/150‑1 (MQ=255)
   gCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAAc                                                                                                                                            <  2:2804911/150‑1 (MQ=255)
   gCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAAc                                                                                                                                            <  2:2750736/150‑1 (MQ=255)
    cGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACt                                                                                                                                           >  1:2114108/1‑150 (MQ=255)
     gTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTg                                                                                                                                          <  1:1506807/150‑1 (MQ=255)
      tAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGc                                                                                                                                         >  2:4283727/1‑150 (MQ=255)
      tAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGc                                                                                                                                         >  2:4284346/1‑150 (MQ=255)
       aGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCa                                                                                                                                        <  2:3166273/150‑1 (MQ=255)
        gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCAc                                                                                                                                       <  1:3216369/150‑1 (MQ=255)
        gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCAc                                                                                                                                       <  1:132929/150‑1 (MQ=255)
        gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCAc                                                                                                                                       <  1:1754619/150‑1 (MQ=255)
        gTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCAc                                                                                                                                       <  1:2165998/150‑1 (MQ=255)
         tGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACg                                                                                                                                      <  1:754850/150‑1 (MQ=255)
          ggAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACgg                                                                                                                                     >  1:3997865/1‑150 (MQ=255)
            taaCGCCTTGTTTTATGCCTGTTGGCACCTACGGCAACGTAAAAGGGATGACGCAGGAAGAAGTTGAAGCCAATGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGAACGGCg                                                                                                                                   <  1:1028464/149‑1 (MQ=255)
            aaaCGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCg                                                                                                                                   >  1:627547/1‑150 (MQ=255)
            aaaCGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCg                                                                                                                                   >  1:643480/1‑150 (MQ=255)
             aaCGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGa                                                                                                                                  <  1:4014612/150‑1 (MQ=255)
              aCGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGAt                                                                                                                                 <  1:4113337/150‑1 (MQ=255)
                  cTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCCCTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGc                                                                                                                             >  2:3287552/1‑150 (MQ=255)
                      ttttATGCCTGTTGGCACCTACGGCACCGTAAAAGTGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACATTCCACCAGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGa                                                                                                                         >  1:2313094/1‑150 (MQ=255)
                      ttttATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGa                                                                                                                         >  1:2313114/1‑150 (MQ=255)
                         tATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGAttt                                                                                                                      >  1:5377/1‑150 (MQ=255)
                                gTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAg                                                                                                               <  2:3580928/150‑1 (MQ=255)
                                   ggCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTgg                                                                                                            >  1:2882227/1‑150 (MQ=255)
                                     cACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGaa                                                                                                          >  1:3247078/1‑150 (MQ=255)
                                       ccTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgg                                                                                                        >  1:3312205/1‑150 (MQ=255)
                                       ccTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAgg                                                                                                        >  1:581885/1‑150 (MQ=255)
                                        cTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAggg                                                                                                       >  1:1403187/1‑150 (MQ=255)
                                          aCGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGc                                                                                                     >  1:157129/1‑150 (MQ=255)
                                            ggCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCg                                                                                                   <  1:2876418/150‑1 (MQ=255)
                                              cACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGAt                                                                                                 >  2:3974210/1‑150 (MQ=255)
                                              cACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGAt                                                                                                 >  2:3065218/1‑150 (MQ=255)
                                                ccGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATcc                                                                                               <  2:920312/150‑1 (MQ=255)
                                                 cGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCt                                                                                              <  1:2488933/150‑1 (MQ=255)
                                                  gTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTc                                                                                             <  1:3805503/150‑1 (MQ=255)
                                                   tAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCa                                                                                            <  1:254579/150‑1 (MQ=255)
                                                   tAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCa                                                                                            >  2:757804/1‑150 (MQ=255)
                                                     aaaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCAcc                                                                                          <  1:2827652/150‑1 (MQ=255)
                                                     aaaGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGGGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCAcc                                                                                          <  1:468323/150‑1 (MQ=255)
                                                       aGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGa                                                                                        <  2:3338969/150‑1 (MQ=255)
                                                         ggATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACt                                                                                      <  2:1165551/150‑1 (MQ=255)
                                                           aTGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTcc                                                                                    <  1:1075465/150‑1 (MQ=255)
                                                            tGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcg                                                                                   <  1:3079575/150‑1 (MQ=255)
                                                             gACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCCTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCcgg                                                                                  >  1:3884733/1‑150 (MQ=255)
                                                               cGCCGTAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGAACGATATTATGCAGAGGAAGGGGCAGATCCTCACCGAATCaggag                                                                                >  2:2519729/1‑148 (MQ=255)
                                                                 ccGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCggcggc                                                                              >  2:1705928/1‑150 (MQ=255)
                                                                   ggAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCtt                                                                            >  1:1857276/1‑150 (MQ=255)
                                                                    gaagaaGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTc                                                                           <  1:4503793/150‑1 (MQ=255)
                                                                      agaagTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCa                                                                         <  2:337498/150‑1 (MQ=255)
                                                                             gAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGACACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTc                                                                  <  2:236687/150‑1 (MQ=255)
                                                                               aGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAg                                                                <  2:1413040/150‑1 (MQ=255)
                                                                                 ccACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGcc                                                              >  2:3352707/1‑150 (MQ=255)
                                                                                 ccACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGcc                                                              >  2:3352733/1‑150 (MQ=255)
                                                                                 ccACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGcc                                                              <  2:4542427/150‑1 (MQ=255)
                                                                                  cACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCt                                                             <  2:2752047/150‑1 (MQ=255)
                                                                                   aCTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCtt                                                            >  1:150335/1‑150 (MQ=255)
                                                                                   aCTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCtt                                                            >  1:2034219/1‑150 (MQ=255)
                                                                                    cTGGCTCTCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGTGACAGTAAATAATGAAACTTCACGGCTATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCATCCttt                                                           <  2:2313094/150‑2 (MQ=255)
                                                                                    cTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTg                                                           <  1:4289751/150‑1 (MQ=255)
                                                                                      ggCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGc                                                         >  2:4446613/1‑150 (MQ=255)
                                                                                        cgcgcAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGa                                                       >  1:2108745/1‑150 (MQ=255)
                                                                                         gcgcAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATTATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTTGCGat                                                      <  2:2382183/150‑1 (MQ=255)
                                                                                            cAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATAtt                                                   <  1:1022647/150‑1 (MQ=255)
                                                                                             aaaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTc                                                  >  2:848205/1‑150 (MQ=255)
                                                                                              aaTTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGTCCATGAAATAATGAAACTGCACGTCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCg                                                 <  2:2879522/150‑1 (MQ=255)
                                                                                                 tATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaa                                              >  1:4091419/1‑150 (MQ=255)
                                                                                                  aTCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaaa                                             >  2:3705351/1‑150 (MQ=255)
                                                                                                   tCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTaaaa                                            >  2:1114104/1‑150 (MQ=255)
                                                                                                    ccTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAAt                                           <  2:2730419/150‑1 (MQ=255)
                                                                                                      tCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCa                                         >  1:226256/1‑150 (MQ=255)
                                                                                                        ggCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCAcc                                       <  2:2261993/150‑1 (MQ=255)
                                                                                                         gCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCg                                      <  2:3858965/150‑1 (MQ=255)
                                                                                                          cAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGa                                     >  2:2692372/1‑150 (MQ=255)
                                                                                                           aaCACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGaa                                    >  1:116661/1‑150 (MQ=255)
                                                                                                              acCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAg                                 <  2:2278625/150‑1 (MQ=255)
                                                                                                                cTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGATACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAggg                               <  2:354132/150‑1 (MQ=255)
                                                                                                                cTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAggg                               >  1:1636939/1‑150 (MQ=255)
                                                                                                                   ccACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCTCGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGt                            <  2:1666158/150‑1 (MQ=255)
                                                                                                                   ccACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGt                            <  2:2315139/150‑1 (MQ=255)
                                                                                                                    cACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                           >  2:159237/1‑150 (MQ=255)
                                                                                                                    cACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                           >  1:3445105/1‑150 (MQ=255)
                                                                                                                    cACCTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTg                           >  1:3445081/1‑150 (MQ=255)
                                                                                                                       cTGTGGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCAc                        <  1:2243896/150‑1 (MQ=255)
                                                                                                                        tgtgGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACt                       >  1:973176/1‑150 (MQ=255)
                                                                                                                          tgGCTGCGCCCGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTc                     >  1:4128901/1‑150 (MQ=255)
                                                                                                                                  cccGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAccc             >  2:2402465/1‑150 (MQ=255)
                                                                                                                                  cccGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAccc             >  2:2402376/1‑150 (MQ=255)
                                                                                                                                  cccGGGCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAAccc             >  1:3968455/1‑150 (MQ=255)
                                                                                                                                     gggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt          <  2:420587/150‑1 (MQ=255)
                                                                                                                                     gggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGAt          <  1:3588729/150‑1 (MQ=255)
                                                                                                                                     gggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGGAACCCGAt          >  2:2783293/1‑150 (MQ=255)
                                                                                                                                      ggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGagc         >  2:1898769/1‑148 (MQ=255)
                                                                                                                                      ggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATc         <  2:4451317/150‑1 (MQ=255)
                                                                                                                                      ggCCAGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTAAGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGAGCACTTCCGTAACCCGATc         >  2:1898809/1‑150 (MQ=255)
                                                                                                                                          aGGAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACg     >  1:1238544/1‑150 (MQ=255)
                                                                                                                                           ggAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACgg    <  2:83587/150‑1 (MQ=255)
                                                                                                                                            gAAATAATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGc   <  2:1368489/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGAAATCATGAAACTGCACGGCGATCTGCACGATTTTATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCG  >  CP009273/421654‑421944

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: