Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,917,554 C→A 90.5% 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,5540CA90.5% 469.1 / 44.2 186L480L (CTG→CTTmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base A (76/92);  total (84/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 186 total aligned reads displayed.

ACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTG  >  CP009273/3917409‑3917696
                                                                                                                                                 |                                                                                                                                              
aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                              <  2:2750062/148‑1 (MQ=255)
aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                              <  1:1404470/148‑1 (MQ=255)
aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                              >  2:2258948/1‑148 (MQ=255)
aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                              <  2:2749624/148‑1 (MQ=255)
aCGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAg                                                                                                                                              <  1:346297/148‑1 (MQ=255)
 cGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAga                                                                                                                                             >  2:315691/1‑148 (MQ=255)
    cGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgc                                                                                                                                          <  2:1049885/148‑1 (MQ=255)
      tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctg                                                                                                                                        <  1:525480/148‑1 (MQ=255)
      tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctg                                                                                                                                        >  2:1677255/1‑148 (MQ=255)
       cAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAgacgctga                                                                                                                                       >  2:1557008/1‑148 (MQ=255)
            tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctc                                                                                                                                  >  1:1595244/1‑148 (MQ=255)
               ttAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacg                                                                                                                               >  2:2954069/1‑148 (MQ=255)
                tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgc                                                                                                                              <  2:2307876/148‑1 (MQ=255)
                tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgc                                                                                                                              >  1:1227224/1‑148 (MQ=255)
                tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgc                                                                                                                              <  2:2312248/148‑1 (MQ=255)
                tAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgc                                                                                                                              <  2:2440526/148‑1 (MQ=255)
                  aTTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGacgctcacgctc                                                                                                                            >  2:2583492/1‑148 (MQ=255)
                   tttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCg                                                                                                                           >  1:2341286/1‑148 (MQ=255)
                      ttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGata                                                                                                                        <  1:3286750/148‑1 (MQ=255)
                          tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATtct                                                                                                                    <  2:1889969/148‑1 (MQ=255)
                             gTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCaa                                                                                                                 <  2:2885542/148‑1 (MQ=255)
                              tCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAg                                                                                                                >  1:2449429/1‑148 (MQ=255)
                               cATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                               >  1:2328868/1‑148 (MQ=255)
                               cATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                               >  2:340968/1‑148 (MQ=255)
                                  ttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtc                                                                                                            >  2:3243195/1‑148 (MQ=255)
                                   tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtct                                                                                                           <  1:545729/148‑1 (MQ=255)
                                      ggACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGt                                                                                                        <  2:2489220/148‑1 (MQ=255)
                                        acgacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTa                                                                                                      >  2:2684042/1‑148 (MQ=255)
                                          gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTaaa                                                                                                    <  1:2360747/148‑1 (MQ=255)
                                           acgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAg                                                                                                   >  2:2502608/1‑148 (MQ=255)
                                           acgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAg                                                                                                   >  1:2970121/1‑148 (MQ=255)
                                              cagcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgc                                                                                                >  2:2725929/1‑148 (MQ=255)
                                              cagcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgc                                                                                                >  2:1250840/1‑148 (MQ=255)
                                                gcagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgcgc                                                                                              >  2:2808610/1‑148 (MQ=255)
                                                    aTCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAATACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCaa                                                                                          >  1:2835065/1‑148 (MQ=255)
                                                    aTCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCaa                                                                                          >  1:2657290/1‑148 (MQ=255)
                                                        tCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGtt                                                                                      >  1:1091075/1‑148 (MQ=255)
                                                            cacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTc                                                                                  <  1:822438/148‑1 (MQ=255)
                                                            cacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTc                                                                                  >  2:3112071/1‑148 (MQ=255)
                                                             acTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCa                                                                                 <  1:407947/148‑1 (MQ=255)
                                                              cTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                >  1:1916012/1‑148 (MQ=255)
                                                              cTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                <  1:980916/148‑1 (MQ=255)
                                                              cTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                <  1:981357/148‑1 (MQ=255)
                                                                ggCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATcc                                                                              >  2:389804/1‑148 (MQ=255)
                                                                  cTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccac                                                                            <  1:857115/148‑1 (MQ=255)
                                                                    tCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccacca                                                                          <  2:654096/148‑1 (MQ=255)
                                                                    tCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccacca                                                                          <  1:2258948/148‑1 (MQ=255)
                                                                      aCGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAgg                                                                        >  1:3302138/1‑148 (MQ=255)
                                                                       cGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGc                                                                       >  1:387191/1‑148 (MQ=255)
                                                                        gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGcc                                                                      <  1:595116/148‑1 (MQ=255)
                                                                           aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAg                                                                   >  2:1110093/1‑148 (MQ=255)
                                                                              gcggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTc                                                                >  2:482789/1‑148 (MQ=255)
                                                                               cggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCa                                                               >  2:143630/1‑148 (MQ=255)
                                                                                 gcgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACg                                                             >  2:306214/1‑148 (MQ=255)
                                                                                  cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGa                                                            <  2:1848401/148‑1 (MQ=255)
                                                                                     aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCg                                                         <  1:3229144/148‑1 (MQ=255)
                                                                                     aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCg                                                         <  1:1980918/148‑1 (MQ=255)
                                                                                      gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                        <  2:2070828/148‑1 (MQ=255)
                                                                                      gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                        <  2:2070687/148‑1 (MQ=255)
                                                                                       tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                       >  2:152183/1‑148 (MQ=255)
                                                                                       tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                       >  2:3075258/1‑148 (MQ=255)
                                                                                       tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCCGGCCCAGTTCACGACCGAt                                                       >  2:2607888/1‑148 (MQ=255)
                                                                                       tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                       <  2:2191488/148‑1 (MQ=255)
                                                                                       tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                       <  2:3262176/148‑1 (MQ=255)
                                                                                         aGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCCGGCCCAGTTCACGACCGAttt                                                     <  1:3244628/148‑1 (MQ=255)
                                                                                         aGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                     <  1:1793903/148‑1 (MQ=255)
                                                                                         aGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                     <  2:1224008/148‑1 (MQ=255)
                                                                                          ggTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTc                                                    >  2:3200231/1‑148 (MQ=255)
                                                                                           gTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCa                                                   <  2:2844701/148‑1 (MQ=255)
                                                                                           gTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCa                                                   <  2:708040/148‑1 (MQ=255)
                                                                                            tGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAg                                                  <  2:606767/148‑1 (MQ=255)
                                                                                            tGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAg                                                  <  2:856531/148‑1 (MQ=255)
                                                                                             gAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGt                                                 <  1:2440653/148‑1 (MQ=255)
                                                                                                cATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCaa                                              >  1:698027/1‑148 (MQ=255)
                                                                                                 aTTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCaaa                                             <  2:3137087/148‑1 (MQ=255)
                                                                                                  ttCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCAGACGCTCACGCGCGATATTCGCAACTTTCTCCTTAAAGCCCGCCCAACCTTCGTCATCCACCACGCCCACTTCACGAACCAATTCAGTaaaaa                                            >  1:2847840/1‑147 (MQ=255)
                                                                                                    cACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                          >  1:1484146/1‑148 (MQ=255)
                                                                                                    cACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                          <  2:1345388/148‑1 (MQ=255)
                                                                                                    cACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                          >  1:2693450/1‑148 (MQ=255)
                                                                                                     aCTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCa                                         <  1:969628/148‑1 (MQ=255)
                                                                                                       ttCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGa                                       <  1:2725929/148‑1 (MQ=255)
                                                                                                        tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                      >  2:264832/1‑148 (MQ=255)
                                                                                                        tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                      >  1:75304/1‑148 (MQ=255)
                                                                                                             tGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCa                                 <  2:1324358/148‑1 (MQ=255)
                                                                                                               cAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCAtt                               >  2:1973962/1‑148 (MQ=255)
                                                                                                               cAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCGCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCAtt                               >  1:827171/1‑148 (MQ=255)
                                                                                                                  cTTCCGCCGACGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATc                            >  1:548091/1‑148 (MQ=255)
                                                                                                                    tCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCtt                          <  1:2173320/148‑1 (MQ=255)
                                                                                                                            cGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGc                  <  1:428619/148‑1 (MQ=255)
                                                                                                                            cGGGGTTACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGc                  <  1:428822/148‑1 (MQ=255)
                                                                                                                                  tACCCAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGa            <  2:768616/148‑1 (MQ=255)
                                                                                                                                      cAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGat        <  1:1553651/148‑1 (MQ=255)
                                                                                                                                      cAGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGat        <  1:1022491/148‑1 (MQ=255)
                                                                                                                                       aGGTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGata       <  1:1040843/148‑1 (MQ=255)
                                                                                                                                         gTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt     >  2:1061098/1‑148 (MQ=255)
                                                                                                                                         gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt     <  1:2738979/148‑1 (MQ=255)
                                                                                                                                         gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt     <  1:20038/148‑1 (MQ=255)
                                                                                                                                         gTCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt     <  2:2007820/148‑1 (MQ=255)
                                                                                                                                          tCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTc    <  2:1559/148‑1 (MQ=255)
                                                                                                                                          tCGATTTAAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTc    >  1:1578744/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTG  >  CP009273/3917409‑3917696

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: