Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,413,119 G→T 46.8% pseudogene (654/960 nt) yhdW → pseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,413,1190GT46.8% 151.8 / 274.9 224pseudogene (654/960 nt)yhdWpseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (62/57);  new base T (58/47);  total (120/104)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CTGCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCTGAAG  >  CP009273/3412972‑3413266
                                                                                                                                                   |                                                                                                                                                   
  gCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTAt                                                                                                                                                   <  2:113340/148‑1 (MQ=255)
        cGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACg                                                                                                                                             <  1:3009293/148‑1 (MQ=255)
                 cGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCt                                                                                                                                    <  2:285315/148‑1 (MQ=255)
                  ggATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCtt                                                                                                                                   <  1:1015027/148‑1 (MQ=255)
                   gATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCttt                                                                                                                                  >  2:1177107/1‑148 (MQ=255)
                   gATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCttt                                                                                                                                  >  2:2332226/1‑148 (MQ=255)
                   gATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCttt                                                                                                                                  <  2:3054758/148‑1 (MQ=255)
                     tcaatcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTc                                                                                                                                >  2:1358350/1‑148 (MQ=255)
                      caatcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCg                                                                                                                               <  2:3368147/148‑1 (MQ=255)
                       aatcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGc                                                                                                                              >  2:1346393/1‑148 (MQ=255)
                        atcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGcc                                                                                                                             <  2:2934157/148‑1 (MQ=255)
                         tcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCa                                                                                                                            <  1:279465/148‑1 (MQ=255)
                         tcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCa                                                                                                                            <  1:1375788/148‑1 (MQ=255)
                          cacaACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCat                                                                                                                           >  2:502071/1‑148 (MQ=255)
                          cacaACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCat                                                                                                                           >  2:903359/1‑148 (MQ=255)
                           acaACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatg                                                                                                                          >  1:2337030/1‑148 (MQ=255)
                            caACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgc                                                                                                                         <  1:45178/148‑1 (MQ=255)
                             aaCTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgct                                                                                                                        >  1:1962699/1‑148 (MQ=255)
                             aaCTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgct                                                                                                                        >  1:1962769/1‑148 (MQ=255)
                              aCTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctg                                                                                                                       >  1:3033426/1‑148 (MQ=255)
                               cTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgctga                                                                                                                      <  2:1767727/148‑1 (MQ=255)
                                 gTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTAGCCatgctgaat                                                                                                                    >  1:742513/1‑148 (MQ=255)
                                  tttGCCCTGCGCATCAACTTAAGCAACCCCGCTGAATGGATTGTCTTACCGGAAGTTATGTCAAATGAACCGCTTGGTCAGGTAGTTCGTCGTGGCGATGATGATTGGTTCTCGATGGTTCGCTGGACGCTTTTCGCCatgctgaatg                                                                                                                   <  2:3120421/146‑1 (MQ=255)
                                  tATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgctgaatg                                                                                                                   >  2:1053742/1‑148 (MQ=255)
                                  tATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatg                                                                                                                   <  2:1602535/148‑1 (MQ=255)
                                  tATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatg                                                                                                                   <  2:765840/148‑1 (MQ=255)
                                   aTGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgctgaatgc                                                                                                                  >  1:1339094/1‑148 (MQ=255)
                                   aTGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatgc                                                                                                                  <  1:1024618/148‑1 (MQ=255)
                                    tGCCCTGCGCATCAATTTTAGCTACCCAGCTGAATGGTTTGTCTTACCGGAAGTTATCTCAAAAGTACCGCTTGGTGCGGTTGTTCGTCGTGGCGATGATGATTGGTTCTCTATTGTACGCTGGACGCTTTTCGCCttgctgaatgct                                                                                                                 <  2:1135849/148‑1 (MQ=255)
                                     gCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgctgaatgctg                                                                                                                >  2:434083/1‑148 (MQ=255)
                                        cTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAag                                                                                                             <  2:693390/148‑1 (MQ=255)
                                          gcgcATCAAATTAAGCAACCCAGCTGACTGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAagag                                                                                                           <  2:1740854/148‑1 (MQ=255)
                                           cgcATCAAATTCAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAgaga                                                                                                          <  1:2590761/148‑1 (MQ=255)
                                             cATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATg                                                                                                        >  1:1550679/1‑148 (MQ=255)
                                              aTCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATgg                                                                                                       <  2:2001466/148‑1 (MQ=255)
                                                 aaaTTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCa                                                                                                    <  1:3138463/148‑1 (MQ=255)
                                                   aTTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATc                                                                                                  >  2:2843327/1‑148 (MQ=255)
                                                   aTTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATc                                                                                                  <  2:2516928/148‑1 (MQ=255)
                                                   aTTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATc                                                                                                  >  1:3314485/1‑148 (MQ=255)
                                                    ttAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCGTCa                                                                                                 >  1:2203373/1‑148 (MQ=255)
                                                    ttAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCa                                                                                                 >  2:260951/1‑148 (MQ=255)
                                                    ttAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCa                                                                                                 <  1:2480382/148‑1 (MQ=255)
                                                    ttAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCa                                                                                                 <  1:1129524/148‑1 (MQ=255)
                                                     tAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCaa                                                                                                >  2:2329114/1‑148 (MQ=255)
                                                     tAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCaa                                                                                                >  2:2829601/1‑148 (MQ=255)
                                                     tAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCaa                                                                                                <  1:3448785/148‑1 (MQ=255)
                                                        gCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTc                                                                                             >  1:409321/1‑148 (MQ=255)
                                                        gCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTc                                                                                             <  1:2541415/148‑1 (MQ=255)
                                                        gCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTc                                                                                             <  1:2611860/148‑1 (MQ=255)
                                                        gCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTc                                                                                             <  1:2843923/148‑1 (MQ=255)
                                                          aaCCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTccc                                                                                           <  2:2367514/148‑1 (MQ=255)
                                                            cccAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAg                                                                                         >  2:3505994/1‑148 (MQ=255)
                                                            cccAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAg                                                                                         >  2:3505906/1‑148 (MQ=255)
                                                               aGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAAc                                                                                      <  2:2718276/148‑1 (MQ=255)
                                                                 cTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGt                                                                                    <  1:215987/148‑1 (MQ=255)
                                                                 cTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGt                                                                                    <  1:1697835/148‑1 (MQ=255)
                                                                 cTGAAGGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGt                                                                                    <  2:2232274/148‑1 (MQ=255)
                                                                   gAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTTCGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCg                                                                                  <  2:2193779/148‑1 (MQ=255)
                                                                   gAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTTCGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCg                                                                                  <  2:1685251/148‑1 (MQ=255)
                                                                    aaTGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGa                                                                                 >  1:2087643/1‑148 (MQ=255)
                                                                      tGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATg                                                                               >  1:108902/1‑148 (MQ=255)
                                                                      tGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATg                                                                               >  1:3201203/1‑148 (MQ=255)
                                                                           tGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGaaaaa                                                                          >  1:2763339/1‑148 (MQ=255)
                                                                             tCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGc                                                                        >  2:2737736/1‑148 (MQ=255)
                                                                              cTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCg                                                                       >  1:3234742/1‑148 (MQ=255)
                                                                                tACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGc                                                                     >  2:470603/1‑148 (MQ=255)
                                                                                   cGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTaa                                                                  >  2:1834892/1‑148 (MQ=255)
                                                                                    ggAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAAt                                                                 <  2:2335178/148‑1 (MQ=255)
                                                                                     gAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATc                                                                >  1:689076/1‑148 (MQ=255)
                                                                                              tcAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGc                                                       <  1:2706473/148‑1 (MQ=255)
                                                                                                aaaaGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCt                                                     <  1:1996097/148‑1 (MQ=255)
                                                                                                aaaaGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCt                                                     >  2:1478882/1‑148 (MQ=255)
                                                                                                 aaaGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTg                                                    >  1:1798737/1‑148 (MQ=255)
                                                                                                   aGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGat                                                  <  1:3546215/148‑1 (MQ=255)
                                                                                                    gAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGata                                                 <  2:577318/148‑1 (MQ=255)
                                                                                                    gAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGata                                                 <  1:2662867/148‑1 (MQ=255)
                                                                                                     aaCCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGatat                                                >  2:3245661/1‑148 (MQ=255)
                                                                                                       ccGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATgg                                              >  2:2108981/1‑148 (MQ=255)
                                                                                                          cTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGcac                                           <  2:456945/148‑1 (MQ=255)
                                                                                                          cTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGcac                                           >  1:2744625/1‑148 (MQ=255)
                                                                                                            tGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTTAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACAt                                         <  1:903359/148‑1 (MQ=255)
                                                                                                            tGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACAt                                         <  1:502071/148‑1 (MQ=255)
                                                                                                            tGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACAt                                         >  1:1201919/1‑148 (MQ=255)
                                                                                                              gTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGTACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGACGATGAAAAATCGGCAAATCCAGCAACGCCTGATATTGcaaaact                                       >  1:2213095/1‑143 (MQ=255)
                                                                                                              gTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATct                                       >  1:2470654/1‑148 (MQ=255)
                                                                                                                     aGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGt                                >  1:2527552/1‑148 (MQ=255)
                                                                                                                       ttCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTaa                              >  1:1320106/1‑148 (MQ=255)
                                                                                                                           tcgtGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaagaa                          <  2:2891028/148‑1 (MQ=255)
                                                                                                                              tGGCGATGATGAATGGTTCTCTTTTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGc                       <  1:2269845/148‑1 (MQ=255)
                                                                                                                                  gatgatGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGAtt                   <  2:1267413/148‑1 (MQ=255)
                                                                                                                                    tgatgaATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTAc                 >  1:3487914/1‑148 (MQ=255)
                                                                                                                                      atgaATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACgg               >  1:1049823/1‑148 (MQ=255)
                                                                                                                                      atgaATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACgg               >  1:818851/1‑148 (MQ=255)
                                                                                                                                        gaATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCa             >  1:752278/1‑148 (MQ=255)
                                                                                                                                              ttCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATc       <  2:3338859/148‑1 (MQ=255)
                                                                                                                                              ttCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATc       <  2:1518992/148‑1 (MQ=255)
                                                                                                                                              ttCTCGATTGTACGCTGGACGCTTTTCGCCATGCAGAATTCTGAAGAGATGGGCATCAATACCCAGAACGACGATGAAAAAGCGGCTAATCAATCAAAGCCTGATATGGCACATATGCTGGGTAAAGAAGGCGATTCCTGCAAGAATc       >  2:726293/1‑148 (MQ=255)
                                                                                                                                               tctcTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCt      >  2:3455956/1‑148 (MQ=255)
                                                                                                                                               tctcGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCt      <  2:1016038/148‑1 (MQ=255)
                                                                                                                                                   gATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCTGAAg  <  2:3325985/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCTGAAG  >  CP009273/3412972‑3413266

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: