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.3% 747.8 / 2.8 238L480L (CTG→CTTmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (108/126);  total (109/129)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-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 238 total aligned reads displayed.

AACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGCGC  >  CP009273/3917408‑3917699
                                                                                                                                                  |                                                                                                                                                 
aaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                                 >  2:1414945/1‑149 (MQ=255)
aaCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                                 <  2:285624/149‑1 (MQ=255)
 aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAga                                                                                                                                                <  2:2195081/149‑1 (MQ=255)
 aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAga                                                                                                                                                >  2:794072/1‑149 (MQ=255)
      gTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctg                                                                                                                                           >  1:2662808/1‑149 (MQ=255)
      gTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctg                                                                                                                                           >  2:2631407/1‑149 (MQ=255)
       tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctga                                                                                                                                          >  2:930742/1‑149 (MQ=255)
       tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctga                                                                                                                                          >  2:3300377/1‑149 (MQ=255)
          gCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctgacgc                                                                                                                                       >  2:2973044/1‑149 (MQ=255)
            gTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctc                                                                                                                                     <  2:173024/149‑1 (MQ=255)
            gTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctc                                                                                                                                     <  2:172883/149‑1 (MQ=255)
             tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAATACGCTGacgctca                                                                                                                                    >  1:1263359/1‑149 (MQ=255)
             tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctca                                                                                                                                    >  1:573536/1‑149 (MQ=255)
             tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctca                                                                                                                                    >  1:691404/1‑149 (MQ=255)
              ggTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcac                                                                                                                                   <  1:3219620/149‑1 (MQ=255)
                 tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgct                                                                                                                                <  2:979582/149‑1 (MQ=255)
                    tttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGa                                                                                                                             <  1:1854905/149‑1 (MQ=255)
                    tttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGa                                                                                                                             <  1:3893485/149‑1 (MQ=255)
                      tttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGata                                                                                                                           <  1:1867297/149‑1 (MQ=255)
                       ttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGatat                                                                                                                          >  2:3876553/1‑149 (MQ=255)
                       ttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGatat                                                                                                                          >  2:3876340/1‑149 (MQ=255)
                         cATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtc                                                                                                                        <  1:2312665/149‑1 (MQ=255)
                         cATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtc                                                                                                                        <  2:1830893/149‑1 (MQ=255)
                           tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                      >  2:2455343/1‑149 (MQ=255)
                           tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                      <  1:3424342/149‑1 (MQ=255)
                           tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                      <  1:3347523/149‑1 (MQ=255)
                           tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                      <  1:1643292/149‑1 (MQ=255)
                             aGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCaa                                                                                                                    >  1:2168532/1‑149 (MQ=255)
                               tCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                                  >  2:3937640/1‑149 (MQ=255)
                               tCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                                  <  1:3305647/149‑1 (MQ=255)
                               tCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                                  >  2:712795/1‑149 (MQ=255)
                                cATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGtt                                                                                                                 <  1:3427468/149‑1 (MQ=255)
                                  tttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtc                                                                                                               >  1:2785529/1‑149 (MQ=255)
                                    tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCTAGTTtctc                                                                                                             >  2:2096462/1‑149 (MQ=255)
                                    tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtctc                                                                                                             >  2:139328/1‑149 (MQ=255)
                                      cGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGt                                                                                                           <  1:2260911/149‑1 (MQ=255)
                                          cgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTaaa                                                                                                       <  2:3215583/149‑1 (MQ=255)
                                            acgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                     >  1:2992616/1‑149 (MQ=255)
                                            acgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                     >  2:3450363/1‑149 (MQ=255)
                                              gcagcaGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgc                                                                                                   >  2:2621043/1‑149 (MQ=255)
                                                 gcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGcc                                                                                                >  2:4080594/1‑149 (MQ=255)
                                                  cagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTTGTTAAAGCGCGccc                                                                                               <  1:420340/149‑1 (MQ=255)
                                                  cagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGccc                                                                                               <  2:2814602/149‑1 (MQ=255)
                                                   agATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCa                                                                                              >  1:2242401/1‑149 (MQ=255)
                                                    gATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCaa                                                                                             >  2:3247737/1‑149 (MQ=255)
                                                    gATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCaa                                                                                             <  1:2323228/149‑1 (MQ=255)
                                                     aTCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCAGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAAc                                                                                            <  1:438156/149‑1 (MQ=255)
                                                      tCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACg                                                                                           <  2:4044036/149‑1 (MQ=255)
                                                           accacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGtt                                                                                                                 <  1:3647352/122‑1 (MQ=255)
                                                           accacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGtt                                                                                                                 >  2:3647352/1‑122 (MQ=255)
                                                            ccacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCCTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTc                                                                                     <  1:2102066/149‑1 (MQ=255)
                                                               cTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATc                                                                                  >  2:3141190/1‑149 (MQ=255)
                                                                tGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATcc                                                                                 >  2:13211/1‑149 (MQ=255)
                                                                 ggCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATcca                                                                                >  2:1289398/1‑149 (MQ=255)
                                                                          ggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                       >  2:2656344/1‑149 (MQ=255)
                                                                           gAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCCGGCCCAg                                                                      >  1:2840627/1‑149 (MQ=255)
                                                                            aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGt                                                                     >  2:1491160/1‑149 (MQ=255)
                                                                             aaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGtt                                                                    >  1:1127601/1‑149 (MQ=255)
                                                                             aaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCACGCCCAGtt                                                                    >  1:370943/1‑149 (MQ=255)
                                                                                cggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCAc                                                                 <  1:1388180/149‑1 (MQ=255)
                                                                                  gcgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGa                                                               <  1:1625529/149‑1 (MQ=255)
                                                                                   cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAc                                                              <  1:2370277/149‑1 (MQ=255)
                                                                                      aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGCTCTTATAAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                           <  1:4034023/149‑1 (MQ=255)
                                                                                      aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                           <  2:37726/149‑1 (MQ=255)
                                                                                      aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                           <  1:252664/149‑1 (MQ=255)
                                                                                            gTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAg                                                     <  1:778775/149‑1 (MQ=255)
                                                                                             tGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGt                                                    >  2:3418545/1‑149 (MQ=255)
                                                                                                  aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                               <  2:1734183/149‑1 (MQ=255)
                                                                                                  aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                               >  2:1788632/1‑149 (MQ=255)
                                                                                                    tCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                             <  1:2455343/149‑1 (MQ=255)
                                                                                                     cACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCa                                            <  2:605911/149‑1 (MQ=255)
                                                                                                        ttCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCACGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                         >  2:3813256/1‑149 (MQ=255)
                                                                                                          cGGCTGCAGCGTCCGCCGTCGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAAAGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                       <  1:2942772/149‑1 (MQ=255)
                                                                                                            gCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCCCGCTCGATATTCTCAAGGTTCGCGTGAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGc                                     >  2:234156/1‑149 (MQ=255)
                                                                                                            gCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGc                                     <  1:27027/149‑1 (MQ=255)
                                                                                                             cTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCa                                    <  1:310345/149‑1 (MQ=255)
                                                                                                             cTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCa                                    <  2:1354451/149‑1 (MQ=255)
                                                                                                                 aGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTAt                                <  1:1004138/149‑1 (MQ=255)
                                                                                                                 aGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTAt                                <  1:1003658/149‑1 (MQ=255)
                                                                                                                 aGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTAt                                <  1:2646363/149‑1 (MQ=255)
                                                                                                                  gCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATc                               <  2:1429523/149‑1 (MQ=255)
                                                                                                                    ttCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCtt                             <  1:3141190/149‑1 (MQ=255)
                                                                                                                     tCCGCCGCCGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTc                            <  2:2204717/149‑1 (MQ=255)
                                                                                                                     tCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCGCACGCTCGCTATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCACGTCACGACCGATTTCAGTCAAACCCCGATCCGCATTATCTTc                            >  1:1057129/1‑149 (MQ=255)
                                                                                                                         ccgACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGt                        >  2:1932163/1‑149 (MQ=255)
                                                                                                                               gggTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATc                  <  2:3088994/149‑1 (MQ=255)
                                                                                                                                ggTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCa                 <  1:2973044/149‑1 (MQ=255)
                                                                                                                                ggTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCa                 <  1:2973108/149‑1 (MQ=255)
                                                                                                                                 gTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAg                >  1:4027534/1‑149 (MQ=255)
                                                                                                                                 gTTACCCAGGTCGATTTAAGAAGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAg                >  1:1991512/1‑149 (MQ=255)
                                                                                                                                    aCCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACg             >  2:270393/1‑149 (MQ=255)
                                                                                                                                     cccAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGa            <  2:2975726/149‑1 (MQ=255)
                                                                                                                                        aGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGatat         <  2:1943663/149‑1 (MQ=255)
                                                                                                                                        aGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGatat         <  2:3985305/149‑1 (MQ=255)
                                                                                                                                        aGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGatat         >  1:1805042/1‑149 (MQ=255)
                                                                                                                                          gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTc       <  1:179809/149‑1 (MQ=255)
                                                                                                                                          gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTc       <  1:1115769/149‑1 (MQ=255)
                                                                                                                                            cGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTg     <  1:1845237/149‑1 (MQ=255)
                                                                                                                                              aTTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTgcg   >  2:1160383/1‑149 (MQ=255)
                                                                                                                                              aTTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTgcg   >  2:439588/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGCGC  >  CP009273/3917408‑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: