Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,203,148 T→G 12.9% K252T (AAG→ACG)  tsaD ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; glycation binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,203,1480TG12.9% 541.5 / 3.6 204K252T (AAG→ACG) tsaDtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; glycation binding protein
Reads supporting (aligned to +/- strand):  ref base T (50/112);  new base G (30/2);  total (90/114)
Fisher's exact test for biased strand distribution p-value = 1.43e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

TATCAGTACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGACGGTC  >  CP009273/3203002‑3203296
                                                                                                                                                  |                                                                                                                                                    
tATCAGTACAAAATTCCGGACGCGCGTAGAACTCTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTg                                                                                                                                                    <  1:3673794/149‑1 (MQ=255)
tATCAGTACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTg                                                                                                                                                    <  1:2794046/149‑1 (MQ=255)
tATCAGTACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTg                                                                                                                                                    <  1:1448784/149‑1 (MQ=255)
tATCAGTACAAAAGTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAAAGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTg                                                                                                                                                    <  1:3440062/149‑1 (MQ=255)
     gTACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCCGCTGATACAGCGCCCCCATGcccat                                                                                                                                               >  2:4121453/1‑145 (MQ=255)
     gTACAAAATTCCGGACGAGCGTATAACACTTCGCCTCGGCGTTTTTTCCTCAGTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCTTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACtt                                                                                                                                               <  2:2850998/149‑1 (MQ=255)
      tACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATACAGCGCCCGCATGCACATa                                                                                                                                              >  2:1594092/1‑149 (MQ=255)
         aaaaTTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAATCCCGTCTGATCCACCGCCCCCTTGCACTTAata                                                                                                                                           >  2:3860222/1‑148 (MQ=255)
           aaTTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAACCCCCTCTGATACAGCGCCCGCCTGCACCTAataaa                                                                                                                                         >  2:3763984/1‑146 (MQ=255)
             ttCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCGTAAAGCCCCTCAGAACCCGCGCCCCCGTGCACATAAtcatca                                                                                                                                       >  2:4117762/1‑149 (MQ=255)
                cGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCg                                                                                                                                    <  2:3376197/149‑1 (MQ=255)
                cGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCg                                                                                                                                    <  2:1465100/149‑1 (MQ=255)
                cGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCg                                                                                                                                    <  2:1234773/149‑1 (MQ=255)
                  gACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTa                                                                                                                                  >  2:300050/1‑149 (MQ=255)
                  gACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCCTCTGATACAGCGCCCGCGTGCCCATAAACATCAGCGTa                                                                                                                                  >  2:2692524/1‑149 (MQ=255)
                  gACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATAATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGAGACAGCCCCCCCGTGCACTTAATCAGCAGCggc                                                                                                                                  >  2:2691162/1‑147 (MQ=255)
                                acTTCGCCGCGGCGTATTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATc                                                                                                                    <  1:923948/149‑1 (MQ=255)
                                 cTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCt                                                                                                                   <  1:1323577/149‑1 (MQ=255)
                                      ccgcGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaaa                                                                                                              <  2:25318/149‑1 (MQ=255)
                                      ccgcGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaaa                                                                                                              <  1:783122/149‑1 (MQ=255)
                                         cggcgTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCGTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGgc                                                                                                           <  2:2717261/149‑1 (MQ=255)
                                             gTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCgg                                                                                                       <  2:749099/149‑1 (MQ=255)
                                              ttttttCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATACAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCggg                                                                                                      >  2:2023551/1‑149 (MQ=255)
                                               tttttCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGc                                                                                                     <  1:2414441/149‑1 (MQ=255)
                                                ttttCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCg                                                                                                    <  2:1167765/149‑1 (MQ=255)
                                                ttttCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCCTATGATACAGCGCCCGCATGAACAAAAACATAAGAGTATCGAACACCGCATCATCAAAGCCGCGGGCg                                                                                                    >  1:129590/1‑149 (MQ=255)
                                                   tcatcaTTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGata                                                                                                 <  1:4121005/149‑1 (MQ=255)
                                                    catcatTTCAGCCAGCTTCGCCCGTACCGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCGGATCCCGCGCCCGCGTGCACGTACTCATCCCCGTATCGACCCCCGCATCGTCAAAGGCGCGGGCGatat                                                                                                >  2:3069464/1‑149 (MQ=255)
                                                    catcatTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCATGCACTTAATCAACAGCGTATCGACCACCCCATCTTCAAAGGCGCGGGCGatat                                                                                                >  2:3786182/1‑149 (MQ=255)
                                                     atcatTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATc                                                                                               <  1:2858856/149‑1 (MQ=255)
                                                     atcatTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCGTGCACATAATCATCAGCGTATCGACCACCGCCTCTTCCAAGGCGCGGGCGCTATc                                                                                               >  2:1118185/1‑149 (MQ=255)
                                                     atcatTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATACAGCGCCCGCGTGCACGTAATCATCAGCGTATAGACCACCCCATCTTCAAAGGCGCGGGCGATATc                                                                                               >  2:2255609/1‑149 (MQ=255)
                                                        atTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGc                                                                                            <  1:1753862/149‑1 (MQ=255)
                                                        atTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCATGCACTTAATAATCAGCGTATCGACCACCGCAAAGTCAAAGGCGCGGGCGATATCAGc                                                                                            >  1:4233250/1‑149 (MQ=255)
                                                         tttCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCa                                                                                           <  1:761122/149‑1 (MQ=255)
                                                             agccagcTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCGTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCAcgcg                                                                                       >  2:3638856/1‑149 (MQ=255)
                                                              gccagcTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCGTAAAGCCCCTCAGATCCAGCGCCCGCGTGCACGTAATAATAAGCGTATAGACCACCGCAACTTCAAAGGCGCGGGCGATATCAGCACGCGt                                                                                      >  1:3464212/1‑149 (MQ=255)
                                                               ccagcTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAATCCCGTCTTATCCATCGCCCGCTTGCACATAATAATAAGCGAATAGACCACCGCATCTTCAAAGGCGCGGGCGATAACAGCAAGCGTc                                                                                     >  2:697067/1‑149 (MQ=255)
                                                               ccagcTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTc                                                                                     <  1:945101/149‑1 (MQ=255)
                                                               ccagcTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTc                                                                                     >  2:2913570/1‑149 (MQ=255)
                                                                    ttCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGgtc                                                                                <  1:472971/149‑1 (MQ=255)
                                                                        cccGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGtcgtcg                                                                            <  1:2635238/149‑1 (MQ=255)
                                                                         ccGTAACGTACGGTTAGCACTCCCGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCgg                                                                           <  2:1402327/149‑1 (MQ=255)
                                                                         ccGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCgg                                                                           <  2:4125708/149‑1 (MQ=255)
                                                                           gTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAATCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCCCAACTTCAAAGGCGCGGGCGATATCAGCACGCGTCAGGTCGTCGGTg                                                                         >  1:1702670/1‑149 (MQ=255)
                                                                           gTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAATCCCCTCTGATCCAGCGCCCGCTTTCACATAATAATCAGCGTATCGACCACCGCATATTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTg                                                                         >  1:1700926/1‑149 (MQ=255)
                                                                           gTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCGTGCACGTAATCAACAGCGTATCGACCCCCCCATCATCAAAGGCGCGGGCGATATCACCACGCGTCTGGTCGTCgggg                                                                         >  1:3553878/1‑147 (MQ=255)
                                                                             aaCGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGcc                                                                       <  2:1399501/149‑1 (MQ=255)
                                                                               cGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGt                                                                     <  2:3830336/149‑1 (MQ=255)
                                                                               cGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGt                                                                     <  2:3830288/149‑1 (MQ=255)
                                                                                gTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTTAAGCCCGTCAGATCCACCGCCCCCGTGCACATAAACAAAAGCGTATAGACCACCGCAAAATCAAAGGCGCGGGCGCTATAAGCACCCGTCTGGTCGTCGGGGCCGtt                                                                    >  1:4160227/1‑149 (MQ=255)
                                                                                 tACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTg                                                                   <  1:1064083/149‑1 (MQ=255)
                                                                                 tACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTg                                                                   >  1:1810675/1‑149 (MQ=255)
                                                                                  aCGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTGAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGt                                                                  <  2:2726803/149‑1 (MQ=255)
                                                                                   cGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTc                                                                 <  1:4038711/149‑1 (MQ=255)
                                                                                   cGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTc                                                                 <  1:526572/149‑1 (MQ=255)
                                                                                    ggTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCa                                                                >  2:2068007/1‑149 (MQ=255)
                                                                                     gTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCGTGCACATAATCATCAGCGTATCGACCACCCCATCATCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGGTGTCAc                                                               >  2:3687398/1‑149 (MQ=255)
                                                                                      ttAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGGCTGGTCGTCGGTGCCGTTGTCACg                                                              >  2:3873848/1‑149 (MQ=255)
                                                                                         gCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAAt                                                           <  2:531562/149‑1 (MQ=255)
                                                                                           aCTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCGTGCACTTAATAATCAGCGTATAGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATgg                                                         >  2:1817863/1‑149 (MQ=255)
                                                                                            cTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAATGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGt                                                        >  2:4039098/1‑149 (MQ=255)
                                                                                            cTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGt                                                        <  1:1211564/149‑1 (MQ=255)
                                                                                            cTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGt                                                        <  1:2748027/149‑1 (MQ=255)
                                                                                             tCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTa                                                       <  2:293039/149‑1 (MQ=255)
                                                                                              cACGCCGCCCGCCATGACCAGTCGCTTAAATCCCGTCTGATCCAGCGCCCGCGTGCACATAAACATCAGCGTAACGACCACCCCATCGTCAAAGGCGCGGGCGCTATCAGCCCCCGGCGGGGCGTCGGTGCCGTTGTCACGAATGGGAt                                                      >  2:809499/1‑149 (MQ=255)
                                                                                              cACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGTTCCAGCGCCCGCGTGCACTTAATCATCAGCGTATCGACCACCGCATCTTAAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAAGGGTAt                                                      >  2:2277537/1‑149 (MQ=255)
                                                                                                      ccgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCcgcg                                              >  2:1790336/1‑149 (MQ=255)
                                                                                                            tGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGtt                                        <  1:625713/149‑1 (MQ=255)
                                                                                                            tGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGtt                                        <  2:3962651/149‑1 (MQ=255)
                                                                                                            tGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGtt                                        <  1:665025/149‑1 (MQ=255)
                                                                                                            tGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGtt                                        <  2:1742512/149‑1 (MQ=255)
                                                                                                             gACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGttt                                       <  1:3992391/149‑1 (MQ=255)
                                                                                                               ccAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTc                                     <  1:2376144/149‑1 (MQ=255)
                                                                                                               ccAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTc                                     <  1:3095420/149‑1 (MQ=255)
                                                                                                                   tCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGc                                 >  2:2935913/1‑149 (MQ=255)
                                                                                                                   tCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGc                                 <  1:301452/149‑1 (MQ=255)
                                                                                                                   tCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGc                                 <  1:300050/149‑1 (MQ=255)
                                                                                                                       ttAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGga                             <  2:2100923/149‑1 (MQ=255)
                                                                                                                       ttAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGga                             <  2:232419/149‑1 (MQ=255)
                                                                                                                        tAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGgag                            >  1:2943702/1‑149 (MQ=255)
                                                                                                                           aGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAg                         <  1:2986363/149‑1 (MQ=255)
                                                                                                                            gCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGc                        <  1:3018490/149‑1 (MQ=255)
                                                                                                                             cccGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCt                       <  1:2952937/149‑1 (MQ=255)
                                                                                                                             cccGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCt                       <  1:2150324/149‑1 (MQ=255)
                                                                                                                                gTCTGATCCAGCGCCCGCTTGCACTTAATCATAAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCAGGAGAAGCTGaa                    >  2:3445230/1‑149 (MQ=255)
                                                                                                                                 tCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGaaa                   <  1:852652/149‑1 (MQ=255)
                                                                                                                                    gATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATcc                >  1:2582512/1‑149 (MQ=255)
                                                                                                                                    gATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATcc                >  1:3830851/1‑149 (MQ=255)
                                                                                                                                     aTCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCa               <  2:2151530/149‑1 (MQ=255)
                                                                                                                                     aTCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCa               <  2:4372706/149‑1 (MQ=255)
                                                                                                                                          gcgcCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGcccc          >  2:2429295/1‑149 (MQ=255)
                                                                                                                                          gcgcCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGcccc          <  2:451581/149‑1 (MQ=255)
                                                                                                                                           cgcCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCg         <  1:3271018/149‑1 (MQ=255)
                                                                                                                                            gcCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCgg        >  2:1674379/1‑149 (MQ=255)
                                                                                                                                             cccGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGa       >  2:4218134/1‑149 (MQ=255)
                                                                                                                                              ccGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGAc      >  1:810780/1‑149 (MQ=255)
                                                                                                                                                gCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGAcgg    <  1:1043315/149‑1 (MQ=255)
                                                                                                                                                 cTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGAcggt   >  2:1089750/1‑149 (MQ=255)
                                                                                                                                                  ttGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGAcggtc  >  2:2557702/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TATCAGTACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCTGGTCGTCGGTGCCGTTGTCACGAATGGTATTTGCCGCGAAGGTTTTCAGGCCGGAGAAGCTGAAATCCAGCCCCGGACGGTC  >  CP009273/3203002‑3203296

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