Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,917,233 G→T 10.1% I587I (ATC→ATA mnmG ← 5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,917,2330GT10.1% 1094.4 / 72.7 308I587I (ATC→ATAmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base G (116/161);  new base T (16/15);  total (132/176)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 313 total aligned reads displayed.

TACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCTGTTCCGCCGC  >  CP009273/3917084‑3917382
                                                                                                                                                     |                                                                                                                                                     
tACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCg                                                                                                                                                       >  1:1455911/1‑150 (MQ=255)
     ggCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCAc                                                                                                                                                  >  2:1886/1‑150 (MQ=255)
       cATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACtt                                                                                                                                                >  1:2876294/1‑150 (MQ=255)
         tttttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCg                                                                                                                                              >  2:1035323/1‑150 (MQ=255)
         tttttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCg                                                                                                                                              <  1:1322550/150‑1 (MQ=255)
          ttttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGt                                                                                                                                             >  1:1067606/1‑150 (MQ=255)
           tttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGtt                                                                                                                                            <  2:1088029/150‑1 (MQ=255)
            ttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTa                                                                                                                                           <  2:1614826/150‑1 (MQ=255)
            ttAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTa                                                                                                                                           <  1:2708666/150‑1 (MQ=255)
               aTGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGaa                                                                                                                                        >  1:773560/1‑150 (MQ=255)
                   gTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGAc                                                                                                                                    <  2:2232766/150‑1 (MQ=255)
                    ttATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGAcc                                                                                                                                   >  1:2242242/1‑150 (MQ=255)
                       tGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGa                                                                                                                                >  1:3384320/1‑150 (MQ=255)
                          gcTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATAc                                                                                                                             >  2:1472563/1‑150 (MQ=255)
                          gcTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATAc                                                                                                                             >  2:3860666/1‑150 (MQ=255)
                            tACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCt                                                                                                                           >  1:3006203/1‑150 (MQ=255)
                              cgacgCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTgg                                                                                                                         <  1:371225/150‑1 (MQ=255)
                                  gcagcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGt                                                                                                                     >  2:2029637/1‑150 (MQ=255)
                                  gcagcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGt                                                                                                                     <  2:336992/150‑1 (MQ=255)
                                   cagcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTa                                                                                                                    <  1:2201161/150‑1 (MQ=255)
                                   cagcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTa                                                                                                                    <  1:1111289/150‑1 (MQ=255)
                                    agcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTaa                                                                                                                   >  1:138405/1‑150 (MQ=255)
                                    agcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTaa                                                                                                                   <  2:733287/150‑1 (MQ=255)
                                    agcaTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTATGGCGATCACTTCGTTAGAAACACAGGATACCTGTCGGTaa                                                                                                                   >  1:3344423/1‑150 (MQ=255)
                                       aTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATcc                                                                                                                >  2:84043/1‑150 (MQ=255)
                                       aTACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATcc                                                                                                                >  1:2319440/1‑150 (MQ=255)
                                           ccTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAgtg                                                                                                            <  1:2377445/150‑1 (MQ=255)
                                            cTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAgtgt                                                                                                           >  2:2055202/1‑150 (MQ=255)
                                             tgtttgttCAGCCCCACCAGCAGAAGGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGAGAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCGGGCGGTAATCCAGTGTc                                                                                                          <  2:2802585/150‑1 (MQ=255)
                                             tGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTc                                                                                                          >  1:834316/1‑150 (MQ=255)
                                                   ttCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGt                                                                                                    <  2:2528193/150‑1 (MQ=255)
                                                     cAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAg                                                                                                  <  1:3607236/150‑1 (MQ=255)
                                                       gCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCa                                                                                                <  2:2361008/150‑1 (MQ=255)
                                                       gCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCa                                                                                                <  2:2169444/150‑1 (MQ=255)
                                                        ccACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAg                                                                                               >  2:1211873/1‑150 (MQ=255)
                                                         cacaCCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAgg                                                                                              >  2:4676931/1‑150 (MQ=255)
                                                           cacCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGgt                                                                                            >  1:2652029/1‑150 (MQ=255)
                                                            acCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGgtg                                                                                           <  1:32468/150‑1 (MQ=255)
                                                               agcagAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTtc                                                                                        >  1:2114792/1‑150 (MQ=255)
                                                               agcagAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTtc                                                                                        >  2:822859/1‑150 (MQ=255)
                                                                gcagAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTtct                                                                                       <  1:3307158/150‑1 (MQ=255)
                                                                 cagAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTtctc                                                                                      <  2:3974553/150‑1 (MQ=255)
                                                                 cagAATGGAGATGGCCGCAGGCGTGACACCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTtctc                                                                                      <  1:432160/150‑1 (MQ=255)
                                                                  agAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCg                                                                                     >  1:150132/1‑150 (MQ=255)
                                                                    caTGGAGCTGGCCGCCGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGtt                                                                                   <  1:4412328/149‑1 (MQ=255)
                                                                     aTGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTa                                                                                  >  1:374088/1‑150 (MQ=255)
                                                                     aTGGAGAGGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTa                                                                                  <  2:3417836/150‑1 (MQ=255)
                                                                         agaTGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCa                                                                              <  2:1301951/150‑1 (MQ=255)
                                                                             ggCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAgctg                                                                          <  2:4565702/150‑1 (MQ=255)
                                                                             ggCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAgctg                                                                          <  2:2819447/150‑1 (MQ=255)
                                                                               ccGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACATGGCGGTAATCCAGTGTCTCGGGTAGCAGGGTGTTCTCGTTACGCAgctgct                                                                        >  1:518101/1‑150 (MQ=255)
                                                                                cGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCtt                                                                       >  2:1867990/1‑150 (MQ=255)
                                                                                cGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCtt                                                                       >  2:175664/1‑150 (MQ=255)
                                                                                  cAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCtttt                                                                     <  1:2002785/150‑1 (MQ=255)
                                                                                  cAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCtttt                                                                     >  1:3975054/1‑150 (MQ=255)
                                                                                  cAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGAtttt                                                                     <  2:2220540/150‑1 (MQ=255)
                                                                                   aGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATGCCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTc                                                                    <  1:639046/150‑1 (MQ=255)
                                                                                      cGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGAt                                                                 >  1:4557971/1‑150 (MQ=255)
                                                                                        tGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGAtct                                                               <  2:4430126/150‑1 (MQ=255)
                                                                                         gACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGAtctc                                                              <  1:4012433/150‑1 (MQ=255)
                                                                                          aCGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCa                                                             >  1:4581100/1‑150 (MQ=255)
                                                                                             ccAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCt                                                          <  2:4241879/150‑1 (MQ=255)
                                                                                             ccAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCt                                                          <  2:3210015/150‑1 (MQ=255)
                                                                                             ccAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCt                                                          >  1:1219342/1‑150 (MQ=255)
                                                                                               aGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTg                                                        <  2:2125743/150‑1 (MQ=255)
                                                                                                 aaaTACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCt                                                      >  1:186255/1‑150 (MQ=255)
                                                                                                 aaaTACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCt                                                      <  1:2946907/150‑1 (MQ=255)
                                                                                                   aTACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTgg                                                    <  2:3941458/150‑1 (MQ=255)
                                                                                                   aTACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTgg                                                    <  2:2518026/150‑1 (MQ=255)
                                                                                                       gcgAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGgcgcg                                                <  2:3118557/150‑1 (MQ=255)
                                                                                                         gAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGcgcgcg                                              <  1:2848062/150‑1 (MQ=255)
                                                                                                           aGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGat                                            >  2:1936782/1‑150 (MQ=255)
                                                                                                           aGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGat                                            >  1:261071/1‑150 (MQ=255)
                                                                                                            gCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGata                                           <  1:1780857/150‑1 (MQ=255)
                                                                                                               tGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATaa                                        <  2:3635265/150‑1 (MQ=255)
                                                                                                                ggCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAAc                                       <  1:4643295/150‑1 (MQ=255)
                                                                                                                  ccGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCt                                     <  2:1014756/150‑1 (MQ=255)
                                                                                                                    gATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTc                                   >  2:3577493/1‑150 (MQ=255)
                                                                                                                     aTAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCTGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCg                                  <  2:3384192/150‑1 (MQ=255)
                                                                                                                     aTAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCg                                  >  1:3486254/1‑150 (MQ=255)
                                                                                                                     aTAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCg                                  <  2:511056/150‑1 (MQ=255)
                                                                                                                       agagGCTGGTTTGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTa                                >  2:1638826/1‑150 (MQ=255)
                                                                                                                       agagGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTa                                <  2:482206/150‑1 (MQ=255)
                                                                                                                          ggCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTAttt                             <  1:4023521/150‑1 (MQ=255)
                                                                                                                              ggTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAAcc                         >  2:2923301/1‑150 (MQ=255)
                                                                                                                              ggTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAAcc                         <  2:578454/150‑1 (MQ=255)
                                                                                                                               gTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCt                        <  2:2978490/150‑1 (MQ=255)
                                                                                                                                 ttGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGa                      <  1:534498/150‑1 (MQ=255)
                                                                                                                                 ttGTGATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGa                      <  1:365654/150‑1 (MQ=255)
                                                                                                                                 ttGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGa                      <  2:3796778/150‑1 (MQ=255)
                                                                                                                                  tgtgATCGTTAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGaa                     >  2:207656/1‑150 (MQ=255)
                                                                                                                                     gATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAAtct                  >  1:1570525/1‑150 (MQ=255)
                                                                                                                                     gATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAAtct                  <  1:2793696/150‑1 (MQ=255)
                                                                                                                                     gATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAAtct                  <  2:138478/150‑1 (MQ=255)
                                                                                                                                        cGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCaa               <  2:635003/150‑1 (MQ=255)
                                                                                                                                           tAAGTTTGGCTATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCt            <  2:2876294/150‑1 (MQ=255)
                                                                                                                                           tAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCt            >  1:2509511/1‑150 (MQ=255)
                                                                                                                                              gTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGAGGTAATCCAGTGTCGCGGGTAGCAGTGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCTGtt         <  2:4022954/150‑1 (MQ=255)
                                                                                                                                                 tGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCTGTTccg      >  1:2284856/1‑150 (MQ=255)
                                                                                                                                                     gATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCTGTTCcgccgc  >  1:1016264/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTTTGGCGATCACTTCGTTAGAAAGACCGGATACCTGGCGGTAATCCAGTGTCGCGGGTAGCAGGGTGTTCTCGTTACGCAGCTGCTTTTCGATCTCATCTTGCTGGCGCGCGATATAACCTTCGTATTTAACCTGAATCTCAACCTGTTCCGCCGC  >  CP009273/3917084‑3917382

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: