Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,917,554 C→A 100% L480L (CTG→CTT)  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,5540CA98.7% 944.5 / 2.3 304L480L (CTG→CTT) mnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (124/176);  total (124/180)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 304 total aligned reads displayed.

AAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGCGC  >  CP009273/3917407‑3917699
                                                                                                                                                   |                                                                                                                                                 
aaaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                                 <  2:2851041/150‑1 (MQ=255)
aaaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                                 >  1:4969893/1‑150 (MQ=255)
 aaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAga                                                                                                                                                >  1:3590064/1‑150 (MQ=255)
  aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgac                                                                                                                                               <  2:1149599/150‑1 (MQ=255)
  aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgac                                                                                                                                               <  1:1897900/150‑1 (MQ=255)
     gCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgct                                                                                                                                            <  2:3219703/150‑1 (MQ=255)
      cGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctg                                                                                                                                           <  2:619651/150‑1 (MQ=255)
       gTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctga                                                                                                                                          >  1:1480305/1‑150 (MQ=255)
       gTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTAACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctga                                                                                                                                          >  1:955856/1‑150 (MQ=255)
        tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctgac                                                                                                                                         <  2:5202005/150‑1 (MQ=255)
        tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctgac                                                                                                                                         >  2:2439894/1‑150 (MQ=255)
        tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctgac                                                                                                                                         <  1:2411480/150‑1 (MQ=255)
          aGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACAGGCATCACGGGACAGCGGCGCAGTCAGGAGAGCATACACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCACGTCGATTTAAgacgctgacgc                                                                                                                                       >  2:5065654/1‑150 (MQ=255)
            cGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCATCTTCCGCCGACGGGGTTACCCAGGTCGATTTAATACGCTGtctctc                                                                                                                                     >  1:2397193/1‑150 (MQ=255)
               ggTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacg                                                                                                                                  >  1:109036/1‑150 (MQ=255)
                gTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgc                                                                                                                                 >  2:2656561/1‑150 (MQ=255)
                 ttAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgct                                                                                                                                >  1:1247222/1‑150 (MQ=255)
                        ttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATAtt                                                                                                                         <  2:1087693/150‑1 (MQ=255)
                        ttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATAtt                                                                                                                         <  2:3537909/150‑1 (MQ=255)
                           aTAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                      <  1:2926903/150‑1 (MQ=255)
                            tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCa                                                                                                                     <  2:1731591/150‑1 (MQ=255)
                            tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCa                                                                                                                     >  1:1121421/1‑150 (MQ=255)
                              aGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAg                                                                                                                   >  1:2232631/1‑150 (MQ=255)
                              aGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAg                                                                                                                   >  1:4801642/1‑150 (MQ=255)
                               gTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                                  >  1:4592147/1‑150 (MQ=255)
                                  aTTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtc                                                                                                               <  1:1005953/150‑1 (MQ=255)
                                   tttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtct                                                                                                              <  1:3019564/150‑1 (MQ=255)
                                     tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGTTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCg                                                                                                            <  2:3923599/150‑1 (MQ=255)
                                     tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCg                                                                                                            <  2:3626334/150‑1 (MQ=255)
                                      ccGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGt                                                                                                           >  1:2099374/1‑150 (MQ=255)
                                       cGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGtt                                                                                                          <  2:5129000/150‑1 (MQ=255)
                                        ggACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGAAGCTCACGCTCGATATTCTCAAGTTTATCGTTa                                                                                                         >  2:3408563/1‑150 (MQ=255)
                                            gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                     >  1:3349822/1‑150 (MQ=255)
                                               gcagcaGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgtgcg                                                                                                  >  2:1284385/1‑150 (MQ=255)
                                               gcagcaGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgcg                                                                                                  >  1:3751260/1‑150 (MQ=255)
                                                 agcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAGTTAAGACGCTGACGCTCACGCGCGATATTCTCAAGTTTATCGGTAAAGCGCGcc                                                                                                >  1:1813594/1‑150 (MQ=255)
                                                       tCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGt                                                                                          >  2:3735715/1‑150 (MQ=255)
                                                       tCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGt                                                                                          >  2:4676470/1‑150 (MQ=255)
                                                         ttCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTc                                                                                        <  1:553585/150‑1 (MQ=255)
                                                          tCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCg                                                                                       >  2:2419987/1‑150 (MQ=255)
                                                               acTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATc                                                                                  <  1:1671143/150‑1 (MQ=255)
                                                                   gcctCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccacc                                                                              <  2:3250556/147‑1 (MQ=255)
                                                                     ttCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAg                                                                            <  2:4779426/150‑1 (MQ=255)
                                                                         cGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGccc                                                                        <  1:4733460/150‑1 (MQ=255)
                                                                           ggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAg                                                                      >  2:4365027/1‑150 (MQ=255)
                                                                           ggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAg                                                                      >  1:3842077/1‑150 (MQ=255)
                                                                              aaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTc                                                                   <  2:4205613/150‑1 (MQ=255)
                                                                                    cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAcc                                                             <  2:5300381/150‑1 (MQ=255)
                                                                                       aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATAGTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                          <  2:1569527/150‑1 (MQ=255)
                                                                                         tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                        <  2:2123724/150‑1 (MQ=255)
                                                                                         tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                        <  1:3457766/150‑1 (MQ=255)
                                                                                         tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                        <  2:2167434/150‑1 (MQ=255)
                                                                                           aGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCa                                                      >  1:216223/1‑150 (MQ=255)
                                                                                            ggTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCCTCCACCAGGCCCAGTTCACGACCGATTTAAg                                                     >  2:2297770/1‑150 (MQ=255)
                                                                                            ggTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAg                                                     <  2:2388887/150‑1 (MQ=255)
                                                                                             gTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGt                                                    <  1:4470993/150‑1 (MQ=255)
                                                                                              tGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTc                                                   >  2:2022945/1‑150 (MQ=255)
                                                                                                 gCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGTTTTCAGTCaaa                                                <  1:2735844/150‑1 (MQ=255)
                                                                                                 gCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCaaa                                                >  1:1975351/1‑150 (MQ=255)
                                                                                                  cATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                               <  2:3472393/150‑1 (MQ=255)
                                                                                                  cATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                               <  1:2847768/150‑1 (MQ=255)
                                                                                                  cATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                               >  2:2708360/1‑150 (MQ=255)
                                                                                                   aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACg                                              <  1:3541136/150‑1 (MQ=255)
                                                                                                   aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACg                                              <  1:2547423/150‑1 (MQ=255)
                                                                                                    ttCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                             >  2:3360546/1‑150 (MQ=255)
                                                                                                        cTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                         <  2:1480305/150‑1 (MQ=255)
                                                                                                        cTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                         <  2:3768648/150‑1 (MQ=255)
                                                                                                          tCGGCTGCAGCTTCCGCCGACGGTGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                       <  1:5313662/150‑1 (MQ=255)
                                                                                                          tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                       <  1:1512781/150‑1 (MQ=255)
                                                                                                            ggCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGc                                     <  2:1637604/150‑1 (MQ=255)
                                                                                                             gCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCa                                    <  2:1247222/150‑1 (MQ=255)
                                                                                                               tGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCAtt                                  <  1:3096142/150‑1 (MQ=255)
                                                                                                                   gCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCt                              >  1:4680766/1‑150 (MQ=255)
                                                                                                                    cTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCCTCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCtt                             >  2:1412382/1‑150 (MQ=255)
                                                                                                                     ttCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAATGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTttcttc                            <  1:679222/150‑1 (MQ=255)
                                                                                                                     ttCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTc                            <  1:1420043/150‑1 (MQ=255)
                                                                                                                     ttCCGCCGACGGGGTTACACAGGTCTTTTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCAAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTc                            <  1:1434639/150‑1 (MQ=255)
                                                                                                                      tCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCAGTATCTTcg                           >  1:3255248/1‑150 (MQ=255)
                                                                                                                          ccgAGGGGGTTACCCTGGTCGATTTATGACGCTGACGCTCACGCTCGTTATTCTCATGTTTCTCGTTATTGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGTTTTCAGTCAAACGCAGATCCGCATTATCTTCGCTTa                       <  2:3895880/150‑1 (MQ=255)
                                                                                                                          ccgACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTa                       <  1:775058/150‑1 (MQ=255)
                                                                                                                           cgACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                      <  2:4957018/150‑1 (MQ=255)
                                                                                                                           cgACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                      <  1:3264282/150‑1 (MQ=255)
                                                                                                                                gggTTACACAGGTCGATTTAAGACGATGACGCTCACGCCCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTACACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCa                 <  2:1805877/150‑1 (MQ=255)
                                                                                                                                 ggTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAg                >  1:2041471/1‑150 (MQ=255)
                                                                                                                                 ggTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAg                <  2:259504/150‑1 (MQ=255)
                                                                                                                                   ttACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGAc              <  1:523542/150‑1 (MQ=255)
                                                                                                                                   ttACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGAc              <  1:4382964/150‑1 (MQ=255)
                                                                                                                                   ttACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCACAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCACTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGAc              <  2:3101556/150‑1 (MQ=255)
                                                                                                                                     aCCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGa            >  2:2958613/1‑150 (MQ=255)
                                                                                                                                      cccAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGat           <  2:722546/150‑1 (MQ=255)
                                                                                                                                      cccAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGat           <  2:773203/150‑1 (MQ=255)
                                                                                                                                      cccAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGat           <  2:1900995/150‑1 (MQ=255)
                                                                                                                                        cAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGatat         <  2:1750420/150‑1 (MQ=255)
                                                                                                                                         aGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt        >  1:3931029/1‑150 (MQ=255)
                                                                                                                                           gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt      <  2:5032477/150‑1 (MQ=255)
                                                                                                                                           gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt      <  1:4265449/150‑1 (MQ=255)
                                                                                                                                           gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt      >  1:2318339/1‑150 (MQ=255)
                                                                                                                                            tCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTg     >  2:4568524/1‑150 (MQ=255)
                                                                                                                                             cGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTgc    >  1:4167188/1‑150 (MQ=255)
                                                                                                                                              gATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTgcg   <  2:3085869/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGCGC  >  CP009273/3917407‑3917699

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: ‑