Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,917,548 C→A 9.5% S482S (TCG→TCT 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,5480CA 9.5% 1044.4 / 52.0 289S482S (TCG→TCTmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base C (140/121);  new base A (16/12);  total (156/133)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

CAGGGGCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGC  >  CP009273/3917400‑3917697
                                                                                                                                                    |                                                                                                                                                     
cAGGGGCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCg                                                                                                                                                      >  1:795138/1‑150 (MQ=255)
    ggCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAttt                                                                                                                                                  <  1:866571/150‑1 (MQ=255)
     gCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTc                                                                                                                                                 >  1:382030/1‑150 (MQ=255)
     gCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTc                                                                                                                                                 <  1:2495811/150‑1 (MQ=255)
        aaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAga                                                                                                                                              >  1:1877825/1‑150 (MQ=255)
        aaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAga                                                                                                                                              <  1:3470792/150‑1 (MQ=255)
          cGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacg                                                                                                                                            >  2:2053939/1‑150 (MQ=255)
                cAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgacg                                                                                                                                      >  2:639455/1‑150 (MQ=255)
                cAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgacg                                                                                                                                      >  2:2919021/1‑150 (MQ=255)
                    gTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctca                                                                                                                                  >  1:3676749/1‑150 (MQ=255)
                     tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcac                                                                                                                                 >  1:2654665/1‑150 (MQ=255)
                      ggTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcacg                                                                                                                                <  1:2437506/150‑1 (MQ=255)
                       gTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcacgc                                                                                                                               >  1:4057580/1‑150 (MQ=255)
                        ttAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcacgct                                                                                                                              <  2:2519858/150‑1 (MQ=255)
                         tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcacgctc                                                                                                                             <  2:1324402/150‑1 (MQ=255)
                          aaTTTTTCATAATTAATTTCCGGTCGACGCAGCAGATATTCACCACTGGATTCACGGGTAAGCGGGGCAGTCAGGTGTGCATTCACTTCGGCTGCTGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCg                                                                                                                            <  1:1120286/150‑1 (MQ=255)
                             ttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGata                                                                                                                         <  1:3142982/150‑1 (MQ=255)
                             ttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGata                                                                                                                         >  1:4213322/1‑150 (MQ=255)
                             ttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGata                                                                                                                         <  2:431945/150‑1 (MQ=255)
                              tttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGatat                                                                                                                        >  2:1197027/1‑150 (MQ=255)
                              tttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGatat                                                                                                                        <  1:1055802/150‑1 (MQ=255)
                                tCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATtc                                                                                                                      >  1:718533/1‑150 (MQ=255)
                                tCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATtc                                                                                                                      >  1:581974/1‑150 (MQ=255)
                                   tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCa                                                                                                                   <  2:3656615/150‑1 (MQ=255)
                                           ttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGTGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtctc                                                                                                           <  2:1219188/150‑1 (MQ=255)
                                           ttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtctc                                                                                                           <  2:1219359/150‑1 (MQ=255)
                                           ttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtctc                                                                                                           <  2:3117325/150‑1 (MQ=255)
                                            tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCg                                                                                                          <  2:2780960/150‑1 (MQ=255)
                                             ccGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGt                                                                                                         <  1:2959916/150‑1 (MQ=255)
                                              cGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGtt                                                                                                        <  1:2032203/150‑1 (MQ=255)
                                                 acgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTaaa                                                                                                     >  1:2493569/1‑150 (MQ=255)
                                                 acgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTaaa                                                                                                     <  2:2390640/150‑1 (MQ=255)
                                                  cgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAg                                                                                                    >  1:1966713/1‑150 (MQ=255)
                                                  cgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAg                                                                                                    >  1:2325231/1‑150 (MQ=255)
                                                  cgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCAGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAg                                                                                                    >  1:4537837/1‑150 (MQ=255)
                                                   gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                   >  2:3009178/1‑150 (MQ=255)
                                                   gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTACTCGTTAAAgc                                                                                                   >  2:882906/1‑150 (MQ=255)
                                                        agcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGcc                                                                                              >  1:1622932/1‑150 (MQ=255)
                                                         gcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGccc                                                                                             <  2:3235262/150‑1 (MQ=255)
                                                         gcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGccc                                                                                             >  2:727563/1‑150 (MQ=255)
                                                          cagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCa                                                                                            <  2:2225994/150‑1 (MQ=255)
                                                             aTCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACg                                                                                         <  1:2433153/150‑1 (MQ=255)
                                                             aTCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACg                                                                                         >  1:2520461/1‑150 (MQ=255)
                                                               cTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGtt                                                                                       >  1:2838716/1‑150 (MQ=255)
                                                                ttCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTc                                                                                      >  2:4084220/1‑150 (MQ=255)
                                                                ttCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTc                                                                                      <  1:2810008/150‑1 (MQ=255)
                                                                 tCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCg                                                                                     <  1:1867981/150‑1 (MQ=255)
                                                                    ccacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCa                                                                                  <  2:510497/150‑1 (MQ=255)
                                                                       cTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATcc                                                                               <  2:3559123/150‑1 (MQ=255)
                                                                        tGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATcca                                                                              >  2:1580282/1‑150 (MQ=255)
                                                                          gCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccacc                                                                            >  2:4515019/1‑150 (MQ=255)
                                                                            ttCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAg                                                                          <  2:2677581/150‑1 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  1:2651377/1‑150 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  2:415371/1‑150 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     <  1:3805670/150‑1 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  2:2143684/1‑150 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  1:1666127/1‑150 (MQ=255)
                                                                                 gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGCGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  1:3705074/1‑150 (MQ=255)
                                                                                   gAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGt                                                                   <  2:3297509/150‑1 (MQ=255)
                                                                                   gAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGt                                                                   <  2:1021378/150‑1 (MQ=255)
                                                                                    aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGtt                                                                  >  1:1490728/1‑150 (MQ=255)
                                                                                    aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGtt                                                                  >  1:1650521/1‑150 (MQ=255)
                                                                                    aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGtt                                                                  <  2:3749182/150‑1 (MQ=255)
                                                                                      aGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCa                                                                >  1:2934188/1‑150 (MQ=255)
                                                                                       gcggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCAc                                                               <  1:4288614/150‑1 (MQ=255)
                                                                                               gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAtt                                                       >  1:4456541/1‑150 (MQ=255)
                                                                                               gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAtt                                                       >  1:2211639/1‑150 (MQ=255)
                                                                                                tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                      <  1:1805805/150‑1 (MQ=255)
                                                                                                 cAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTc                                                     >  2:1369744/1‑150 (MQ=255)
                                                                                                 cAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTc                                                     >  1:3104531/1‑150 (MQ=255)
                                                                                                    ttGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGt                                                  <  2:3645670/149‑1 (MQ=255)
                                                                                                       aGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCaa                                               >  2:1914570/1‑150 (MQ=255)
                                                                                                          aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACg                                            >  2:1846959/1‑150 (MQ=255)
                                                                                                              aCTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCACGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGa                                        <  1:4149286/150‑1 (MQ=255)
                                                                                                              aCTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGa                                        <  1:2056918/150‑1 (MQ=255)
                                                                                                              aCTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGa                                        <  2:795138/150‑1 (MQ=255)
                                                                                                               cTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                       >  2:1247579/1‑150 (MQ=255)
                                                                                                                 tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                     >  2:4139820/1‑150 (MQ=255)
                                                                                                                 tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                     <  2:1632407/150‑1 (MQ=255)
                                                                                                                   ggCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGc                                   >  1:340401/1‑150 (MQ=255)
                                                                                                                    gCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCa                                  <  2:1776428/150‑1 (MQ=255)
                                                                                                                     cTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCAt                                 >  2:1700051/1‑150 (MQ=255)
                                                                                                                         aGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATc                             >  1:3280449/1‑150 (MQ=255)
                                                                                                                               cgccgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTcgcg                       <  1:3544053/150‑1 (MQ=255)
                                                                                                                                  cgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                    >  1:1983792/1‑150 (MQ=255)
                                                                                                                                  cgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                    <  1:838742/150‑1 (MQ=255)
                                                                                                                                  cgACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                    >  1:2269866/1‑150 (MQ=255)
                                                                                                                                   gACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGc                   >  1:3916441/1‑150 (MQ=255)
                                                                                                                                     cGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCAt                 >  2:835800/1‑150 (MQ=255)
                                                                                                                                      ggggTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATc                >  2:4016648/1‑150 (MQ=255)
                                                                                                                                      ggggTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATc                <  2:1175792/150‑1 (MQ=255)
                                                                                                                                       gggTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCa               >  2:469922/1‑150 (MQ=255)
                                                                                                                                          ttACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGAc            <  1:3204602/150‑1 (MQ=255)
                                                                                                                                           tACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACg           <  1:2289058/150‑1 (MQ=255)
                                                                                                                                              ccAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGata        >  2:4477395/1‑150 (MQ=255)
                                                                                                                                              ccAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGata        >  2:3191480/1‑150 (MQ=255)
                                                                                                                                                aGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt      >  2:685109/1‑150 (MQ=255)
                                                                                                                                                 ggTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCCAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTc     <  1:9364/150‑1 (MQ=255)
                                                                                                                                                  gTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt    <  2:2325231/150‑1 (MQ=255)
                                                                                                                                                    cGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTgc  >  2:1675472/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGGGCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGC  >  CP009273/3917400‑3917697

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: