Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,382,218 C→A 94.5% noncoding (42/76 nt) glyV → tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,382,2180CA94.5% 469.8 / 14.9 175noncoding (42/76 nt)glyVtRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base A (64/100);  total (71/104)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

AAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTC  >  CP009273/4382072‑4382362
                                                                                                                                                  |                                                                                                                                                
 aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                >  1:2008995/1‑148 (MQ=255)
 aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                >  1:826030/1‑148 (MQ=255)
 aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                >  1:2658179/1‑148 (MQ=255)
 aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                >  1:2987531/1‑148 (MQ=255)
  aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAggg                                                                                                                                               >  1:1701558/1‑148 (MQ=255)
     aaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGACTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTc                                                                                                                                            >  1:1706920/1‑148 (MQ=255)
      aCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcg                                                                                                                                           <  1:1558894/148‑1 (MQ=255)
       cTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGGGGCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcgc                                                                                                                                          >  1:3389496/1‑148 (MQ=255)
        ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcgcg                                                                                                                                         >  1:2482100/1‑148 (MQ=255)
        ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcgcg                                                                                                                                         <  1:3080397/148‑1 (MQ=255)
        ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcgcg                                                                                                                                         <  2:2781358/148‑1 (MQ=255)
        ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGGAGAGCACGACCTTGCCAAGGTAGGGGTcgcg                                                                                                                                         <  2:1579393/148‑1 (MQ=255)
           tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  2:1380537/148‑1 (MQ=255)
           tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  2:2497361/148‑1 (MQ=255)
           tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  1:35432/148‑1 (MQ=255)
           tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  1:35647/148‑1 (MQ=255)
               gggggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGa                                                                                                                                  <  1:923072/148‑1 (MQ=255)
                  ggTTGCAGAGGGAAAGATTTCTCGTATAATGCTCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGtc                                                                                                                               <  2:734245/148‑1 (MQ=255)
                  ggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGtc                                                                                                                               <  1:2516271/148‑1 (MQ=255)
                     tGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCg                                                                                                                            <  2:3092253/148‑1 (MQ=255)
                       cAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                          >  1:949613/1‑148 (MQ=255)
                       cAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                          <  1:897875/148‑1 (MQ=255)
                         gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                        >  1:3025798/1‑148 (MQ=255)
                         gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                        >  1:3025794/1‑148 (MQ=255)
                          agGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                       >  2:3520374/1‑148 (MQ=255)
                            ggAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCg                                                                                                                     <  2:596369/148‑1 (MQ=255)
                              aaaGATTTCTGGTAGAGTTCGCCTCCCGTTTCGCCGCAGAAATGCGCAAATTACGCAAGCTCAATTAAGTAGGACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                                   <  2:3389496/148‑1 (MQ=18)
                              aaaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                                   <  2:1740186/148‑1 (MQ=255)
                               aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                                  <  1:1951687/148‑1 (MQ=255)
                               aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                                  <  1:3243127/148‑1 (MQ=255)
                                   tttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                              >  2:2655189/1‑148 (MQ=255)
                                    ttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                                                                                             >  1:3233523/1‑148 (MQ=255)
                                      ctcGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                                                                                                           <  2:280767/148‑1 (MQ=255)
                                             aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaagt                                                                                                    >  1:2155204/1‑146 (MQ=255)
                                             aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAg                                                                                                    <  2:78057/148‑1 (MQ=255)
                                             aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAg                                                                                                    <  1:2655189/148‑1 (MQ=255)
                                              aTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGt                                                                                                   <  2:3033823/148‑1 (MQ=255)
                                               tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                                                                                                  <  2:949613/148‑1 (MQ=255)
                                               tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                                                                                                  <  1:1268919/148‑1 (MQ=255)
                                               tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTTaaaacac                                                                                                  >  1:2875596/1‑145 (MQ=255)
                                                    cTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGCTGCa                                                                                             >  2:879909/1‑148 (MQ=255)
                                                     tCCCGTAACGAAGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACTACATTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaa                                                                                            <  2:3178203/148‑1 (MQ=255)
                                                       ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAg                                                                                          <  2:2166855/148‑1 (MQ=255)
                                                           aaCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCTTTTCACGCTCCAAAATTTGAAAAGTGCTGAAAAGcaca                                                                                      >  2:2666571/1‑148 (MQ=255)
                                                            acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAg                                                                                     <  2:3158750/148‑1 (MQ=255)
                                                             cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGac                                                                                   <  2:2155204/148‑1 (MQ=255)
                                                                cgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcacataaaa                                                                                 >  1:1903183/1‑143 (MQ=255)
                                                                    gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCACCAAAATTTGAAAATAGCTGCAAAGCACAGACCACCCa                                                                             >  1:2188970/1‑148 (MQ=255)
                                                                     aaaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAAGTGATGAAAAGCACAGACCACCCaa                                                                            >  1:415662/1‑148 (MQ=255)
                                                                     aaaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCCCTCCAAAAGTTGAAAAGGCCGTCACAGCACAGCCCACCCaa                                                                            >  1:373018/1‑148 (MQ=18)
                                                                       aTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGc                                                                          >  2:1943760/1‑148 (MQ=255)
                                                                         gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg                                                                        <  2:2331710/148‑1 (MQ=255)
                                                                         gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg                                                                        <  2:1957776/148‑1 (MQ=255)
                                                                          cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                                                      <  1:3098113/148‑1 (MQ=255)
                                                                           gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                      <  2:394039/148‑1 (MQ=255)
                                                                           gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                      <  1:3152497/148‑1 (MQ=255)
                                                                            aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATGTGCTGAAAAGCACAGACCACCCACGCgggcc                                                                     >  2:3324770/1‑146 (MQ=255)
                                                                            aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGGTGCTGCACAGCACCGACCCCCCACGCGGGaa                                                                     >  2:2744744/1‑148 (MQ=18)
                                                                             aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAt                                                                    >  1:1406941/1‑148 (MQ=255)
                                                                                 tACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCt                                                                <  1:690366/148‑1 (MQ=255)
                                                                                  aCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTc                                                               <  2:3252515/148‑1 (MQ=255)
                                                                                    gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCCCAGACCACCCAAGCGGGAATAGCTCAg                                                             >  2:203042/1‑148 (MQ=255)
                                                                                    gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATCGCTCAg                                                             >  2:202686/1‑148 (MQ=255)
                                                                                        gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                         <  1:1312047/148‑1 (MQ=255)
                                                                                         cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                        <  2:1065420/148‑1 (MQ=255)
                                                                                         cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                        <  1:241587/148‑1 (MQ=255)
                                                                                             aTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                                    <  1:2151626/148‑1 (MQ=255)
                                                                                               tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCa                                                  <  1:777999/148‑1 (MQ=255)
                                                                                               tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCa                                                  <  1:3268547/148‑1 (MQ=255)
                                                                                                  gtagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGa                                               <  2:2818963/148‑1 (MQ=255)
                                                                                                   tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                                         <  1:3005194/73‑1 (MQ=39)
                                                                                                   tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATGTGGTGCATAGCCCAGCCCACCCAAGCGGGAATAGCGCTGTTGGTAGAGCACGAc                                              >  1:1605063/1‑148 (MQ=14)
                                                                                                             gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAAGTGGTGAAAATCACAGACAACCCAATCGGGAATAGATCAGTTGGTTGAGCACGACCTTGCCAAgg                                    >  2:27609/1‑148 (MQ=11)
                                                                                                             gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                    <  2:3436798/148‑1 (MQ=255)
                                                                                                             gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                    <  1:1457187/148‑1 (MQ=255)
                                                                                                                   cTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                              <  2:2397275/148‑1 (MQ=255)
                                                                                                                     cAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                            >  1:2735472/1‑148 (MQ=255)
                                                                                                                      aGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                           <  2:2218499/148‑1 (MQ=255)
                                                                                                                        ttGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGATGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                         >  1:2910258/1‑148 (MQ=255)
                                                                                                                            tAGAGCACGACCTTGCCAAGGTAGGGGTCACGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                     <  2:3156766/148‑1 (MQ=255)
                                                                                                                              gagCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                   >  1:1966978/1‑148 (MQ=255)
                                                                                                                                   cGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGATGAAAAGCACAGACCACCCAAGCGGGAATAGCTCGGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGCGTTCGAGtct              >  1:1656113/1‑148 (MQ=11)
                                                                                                                                   cGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTGGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct              <  2:802687/148‑1 (MQ=255)
                                                                                                                                     aCCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg            <  2:1966978/148‑1 (MQ=255)
                                                                                                                                     aCCTTGCCAAGGGAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg            <  1:1512592/148‑1 (MQ=255)
                                                                                                                                      ccTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt           >  1:1308027/1‑148 (MQ=255)
                                                                                                                                      ccTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt           <  1:3226420/148‑1 (MQ=255)
                                                                                                                                       ctggCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt          <  2:3359136/145‑1 (MQ=255)
                                                                                                                                        ttGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt         <  2:2987531/148‑1 (MQ=255)
                                                                                                                                        ttGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt         >  1:2985928/1‑148 (MQ=255)
                                                                                                                                          gCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc       <  1:2118591/148‑1 (MQ=255)
                                                                                                                                          gCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc       <  1:2871863/148‑1 (MQ=255)
                                                                                                                                          gCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc       <  2:1580558/148‑1 (MQ=255)
                                                                                                                                          gCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc       <  1:514795/148‑1 (MQ=255)
                                                                                                                                            cAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTGGGGGTCGCGCGTTCGAGTCTCGTTTGCCg     >  1:1207174/1‑148 (MQ=14)
                                                                                                                                            cAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCg     <  1:107610/148‑1 (MQ=255)
                                                                                                                                            cAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCg     <  1:2278494/148‑1 (MQ=255)
                                                                                                                                            cAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCg     <  1:107580/148‑1 (MQ=255)
                                                                                                                                              aGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt   <  1:3043533/148‑1 (MQ=255)
                                                                                                                                               ggTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc  >  1:1857669/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTC  >  CP009273/4382072‑4382362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: