Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,398,321 G→T 5.3% S617S (TCG→TCT rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,398,3210GT 5.3% 1053.7 / 18.3 272S617S (TCG→TCTrnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base G (124/133);  new base T (7/8);  total (131/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

CGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTG  >  CP009273/4398173‑4398468
                                                                                                                                                    |                                                                                                                                                   
 gCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGa                                                                                                                                                    >  1:601260/1‑149 (MQ=255)
  ccGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGAt                                                                                                                                                   <  2:2691677/149‑1 (MQ=255)
    gATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATgg                                                                                                                                                 <  1:225801/149‑1 (MQ=255)
    gATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATgg                                                                                                                                                 >  1:1703033/1‑149 (MQ=255)
     aTTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGc                                                                                                                                                >  2:1301014/1‑149 (MQ=255)
         gtcgtTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGaa                                                                                                                                            >  1:430442/1‑149 (MQ=255)
         gtcgtTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGaa                                                                                                                                            >  1:1983100/1‑149 (MQ=255)
          tcgtTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAAc                                                                                                                                           >  2:2840962/1‑149 (MQ=255)
            gtTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAAcgt                                                                                                                                         >  2:1525419/1‑149 (MQ=255)
             ttATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAAcgtc                                                                                                                                        >  1:3786862/1‑149 (MQ=255)
                tCCAGTCCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTg                                                                                                                                     >  1:1422552/1‑149 (MQ=255)
                    gACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCcga                                                                                                                                 >  2:3018397/1‑149 (MQ=255)
                     aCCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCcgac                                                                                                                                >  1:355665/1‑149 (MQ=255)
                     aCCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCcgac                                                                                                                                >  2:2127319/1‑149 (MQ=255)
                          aCGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGc                                                                                                                           <  2:1417859/149‑1 (MQ=255)
                            gCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCaa                                                                                                                         <  1:3891314/149‑1 (MQ=255)
                               gCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAAcgc                                                                                                                      >  2:3155597/1‑149 (MQ=255)
                               gCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAAcgc                                                                                                                      >  2:2918038/1‑149 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               >  2:1653104/1‑149 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               >  1:2680651/1‑149 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               >  1:777455/1‑149 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               >  2:3702356/1‑149 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               <  1:2165833/149‑1 (MQ=255)
                                      gcCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGAtgt                                                                                                               <  2:1770580/149‑1 (MQ=255)
                                         aTTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCACCGCGCGATGTGGc                                                                                                            >  1:804857/1‑149 (MQ=255)
                                          ttAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCt                                                                                                           >  2:781563/1‑149 (MQ=255)
                                             aaTATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGAc                                                                                                        <  1:2604182/149‑1 (MQ=255)
                                               tatCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTg                                                                                                      >  1:2214550/1‑149 (MQ=255)
                                                 tCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGc                                                                                                    <  1:357276/149‑1 (MQ=255)
                                                  ctgctgGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAAc                                                                                                                        <  2:3119495/128‑1 (MQ=255)
                                                      tgGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAg                                                                                               <  2:1161552/149‑1 (MQ=255)
                                                          gAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGtgtg                                                                                           >  1:508610/1‑149 (MQ=255)
                                                          gAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGtgtg                                                                                           >  2:322883/1‑149 (MQ=255)
                                                            aaGAGCAGGGGCATCAGGGCAACACCACTGAAACCGCCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGAc                                                                                         <  2:2848968/149‑1 (MQ=255)
                                                             agagCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGACGCACCGCGCGATGTGGCTGACTGGCTGAAGTGTGACt                                                                                        >  2:1566020/1‑149 (MQ=255)
                                                             agagCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACt                                                                                        >  1:1068265/1‑149 (MQ=255)
                                                             agagCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACt                                                                                        >  2:1937783/1‑149 (MQ=255)
                                                                gCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCa                                                                                     >  1:976048/1‑149 (MQ=255)
                                                                 cAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCAt                                                                                    <  1:1771673/149‑1 (MQ=255)
                                                                 cAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCAt                                                                                    <  1:1569809/149‑1 (MQ=255)
                                                                  aGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCATCTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATg                                                                                   <  2:3834325/149‑1 (MQ=255)
                                                                   ggggCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGc                                                                                  <  2:334742/149‑1 (MQ=255)
                                                                   ggggCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGc                                                                                  >  1:1598655/1‑149 (MQ=255)
                                                                         tCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGAcc                                                                            <  1:2719964/149‑1 (MQ=255)
                                                                              gCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCaggta                                                                       <  1:893216/149‑1 (MQ=255)
                                                                                 acacCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCaggtaggt                                                                    <  1:72487/149‑1 (MQ=255)
                                                                                  caccacTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAggtaggta                                                                   <  1:2691877/149‑1 (MQ=255)
                                                                                   accacTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTaa                                                                  <  1:1525419/149‑1 (MQ=255)
                                                                                    ccacTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAAc                                                                 <  1:628433/149‑1 (MQ=255)
                                                                                      acTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGt                                                               <  2:1993413/149‑1 (MQ=255)
                                                                                      acTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGt                                                               <  2:4030727/149‑1 (MQ=255)
                                                                                        tGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCt                                                             <  2:242047/149‑1 (MQ=255)
                                                                                        tGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCt                                                             <  1:2528770/149‑1 (MQ=255)
                                                                                         gAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCtt                                                            >  1:3700654/1‑149 (MQ=255)
                                                                                         gAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCtt                                                            <  2:2883973/149‑1 (MQ=255)
                                                                                          aaaCCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCttt                                                           <  1:118357/149‑1 (MQ=255)
                                                                                            aCCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCCACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTaa                                                         >  2:69894/1‑149 (MQ=255)
                                                                                              cggcggCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAg                                                       <  1:1908990/149‑1 (MQ=255)
                                                                                                    cTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTaa                                                 <  1:1747796/149‑1 (MQ=255)
                                                                                                    cTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTaa                                                 <  1:1011354/149‑1 (MQ=255)
                                                                                                    cTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTaa                                                 <  2:2574016/149‑1 (MQ=255)
                                                                                                     tACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAAt                                                <  1:2918038/149‑1 (MQ=255)
                                                                                                      aCCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAAtt                                               <  2:2486191/149‑1 (MQ=255)
                                                                                                       ccATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAAttt                                              <  1:2127319/149‑1 (MQ=255)
                                                                                                        cATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTc                                             >  1:3000187/1‑149 (MQ=255)
                                                                                                               cGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTc                                      >  1:1868132/1‑149 (MQ=255)
                                                                                                                 aTGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCAc                                    <  2:3662232/149‑1 (MQ=255)
                                                                                                                 aTGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCAc                                    <  1:564074/149‑1 (MQ=255)
                                                                                                                 aTGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCAc                                    <  1:1194632/149‑1 (MQ=255)
                                                                                                                    gAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtgg                                 >  1:3262239/1‑149 (MQ=255)
                                                                                                                    gAAGAGATGCTGCAACTGGGTCAGCACTGTACGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtgg                                 >  1:3262196/1‑149 (MQ=255)
                                                                                                                      agagATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggct                               <  1:2454720/149‑1 (MQ=255)
                                                                                                                       gagaTGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggctt                              >  2:1358919/1‑149 (MQ=255)
                                                                                                                        agaTGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggcttt                             <  2:1422552/149‑1 (MQ=255)
                                                                                                                          aTGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggctttgg                           <  2:3926655/149‑1 (MQ=255)
                                                                                                                          aTGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggctttgg                           >  1:3299679/1‑149 (MQ=255)
                                                                                                                            gctgcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACtggctttggct                         >  1:365039/1‑149 (MQ=255)
                                                                                                                              tgcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttc                       <  1:1303351/149‑1 (MQ=255)
                                                                                                                               gcAACTGGGTCAGCACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttct                      <  1:2792365/149‑1 (MQ=255)
                                                                                                                               gcAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttct                      <  1:2019825/149‑1 (MQ=255)
                                                                                                                                cAACTGGGTCAGCACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttctt                     <  2:882306/149‑1 (MQ=255)
                                                                                                                                cAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttctt                     <  2:3786862/149‑1 (MQ=255)
                                                                                                                                cAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttctt                     <  1:3527310/149‑1 (MQ=255)
                                                                                                                                cAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttctt                     <  1:1126265/149‑1 (MQ=255)
                                                                                                                                cAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGcttctt                     <  1:542869/149‑1 (MQ=255)
                                                                                                                                 aaCTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCttcttc                    >  2:2225435/1‑149 (MQ=255)
                                                                                                                                     gggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTcgtcc                >  2:1081020/1‑149 (MQ=255)
                                                                                                                                     gggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTcgtcc                >  1:39620/1‑149 (MQ=255)
                                                                                                                                     gggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTcgtcc                <  1:401424/149‑1 (MQ=255)
                                                                                                                                      ggTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCACCGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTcgtccg               >  2:509664/1‑149 (MQ=255)
                                                                                                                                            cACTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGga         <  2:785231/149‑1 (MQ=255)
                                                                                                                                            cACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGga         <  1:2874656/149‑1 (MQ=255)
                                                                                                                                             aCTGTTCTATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGgac        >  2:2743465/1‑149 (MQ=255)
                                                                                                                                              cTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGgacg       >  1:2380866/1‑149 (MQ=255)
                                                                                                                                              cTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGgacg       >  2:1430251/1‑149 (MQ=255)
                                                                                                                                               tGTTCGATTGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGgacga      <  2:1268748/149‑1 (MQ=255)
                                                                                                                                                 ttCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACt    <  1:3155618/149‑1 (MQ=255)
                                                                                                                                                  tCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACtt   >  2:935393/1‑149 (MQ=255)
                                                                                                                                                   cGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTg  <  1:3739700/149‑1 (MQ=255)
                                                                                                                                                   cGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTg  <  1:445276/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAAAGAGCAGGGGCATCAGGGCAACACCACTGAAACCGGCGGCTACCATTATTCGATGGAAGAGATGCTGCAACTGGGTCAGCACTGTTCGATGGCGGAACGTCGTGCCGACGAAGCAACGCGCGATGTGGCTGACTGGCTGAAGTGTGACTTCATGCTCGACCAGGTAGGTAACGTCTTTAAAGGCGTAATTTCCAGCGTCACTGGCTTTGGCTTCTTCGTCCGTCTGGACGACTTG  >  CP009273/4398173‑4398468

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: