Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,152,278 G→T 93.9% E16* (GAA→TAA)  ycfH → putative DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,152,2780GT93.9% 839.6 / 37.8 257E16* (GAA→TAA) ycfHputative DNase
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base T (117/124);  total (125/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 260 total aligned reads displayed.

GCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGT  >  CP009273/1152131‑1152424
                                                                                                                                                   |                                                                                                                                                  
 cATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAAt                                                                                                                                                 <  1:3481134/150‑1 (MQ=255)
  aTCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATc                                                                                                                                                >  1:3147285/1‑150 (MQ=255)
   tCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCt                                                                                                                                               <  1:2427584/150‑1 (MQ=255)
   tCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCt                                                                                                                                               <  1:2090749/150‑1 (MQ=255)
   tCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCt                                                                                                                                               >  1:73277/1‑150 (MQ=255)
    cAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCtt                                                                                                                                              <  1:3889621/150‑1 (MQ=255)
     aaCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCttt                                                                                                                                             <  2:1263149/150‑1 (MQ=255)
     aaCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCttt                                                                                                                                             <  1:2869170/150‑1 (MQ=255)
        cgcgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCa                                                                                                                                          <  1:3932886/150‑1 (MQ=255)
         gcgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCAt                                                                                                                                         <  2:1234518/150‑1 (MQ=255)
          cgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATa                                                                                                                                        <  1:3639017/150‑1 (MQ=255)
           gAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATaa                                                                                                                                       <  1:3957983/150‑1 (MQ=255)
           gAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACGCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATaa                                                                                                                                       <  2:3139970/150‑1 (MQ=255)
            aGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATATTTAAGAGAGACATCAAGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAg                                                                                                                                      >  1:3281011/1‑150 (MQ=255)
                 ctcATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGt                                                                                                                                 >  1:2546563/1‑150 (MQ=255)
                  tcatcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTg                                                                                                                                >  2:699425/1‑150 (MQ=255)
                   catcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTgg                                                                                                                               <  1:3861926/150‑1 (MQ=255)
                   catcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTgg                                                                                                                               >  2:4049073/1‑150 (MQ=255)
                   catcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTgg                                                                                                                               >  1:2488312/1‑150 (MQ=255)
                    atcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGa                                                                                                                              <  1:177106/150‑1 (MQ=255)
                    atcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGa                                                                                                                              >  1:1012120/1‑150 (MQ=255)
                       aCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGa                                                                                                                           >  1:400707/1‑150 (MQ=255)
                        ccGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGAc                                                                                                                          >  1:4169667/1‑150 (MQ=255)
                             cTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCt                                                                                                                     <  1:1010522/150‑1 (MQ=255)
                             cTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCt                                                                                                                     >  1:463543/1‑150 (MQ=255)
                                 tACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCg                                                                                                                 >  1:2692029/1‑150 (MQ=255)
                                  aCTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGa                                                                                                                <  2:2487387/150‑1 (MQ=255)
                                  aCTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGa                                                                                                                >  2:3365665/1‑150 (MQ=255)
                                    tGCGTATTGAGCATTACCCGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGaaa                                                                                                              >  1:1593032/1‑150 (MQ=255)
                                      cGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgc                                                                                                            >  1:1690232/1‑150 (MQ=255)
                                       gTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgcc                                                                                                           <  1:3657996/150‑1 (MQ=255)
                                           tGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGccgccg                                                                                                       <  2:1457264/150‑1 (MQ=255)
                                              gCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAc                                                                                                    <  2:4647911/150‑1 (MQ=255)
                                               cATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcg                                                                                                   >  1:4304184/1‑150 (MQ=255)
                                                aTTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgc                                                                                                  >  2:3821191/1‑150 (MQ=255)
                                                 ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgcg                                                                                                 <  2:3906221/150‑1 (MQ=255)
                                                 ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgcg                                                                                                 <  1:394468/150‑1 (MQ=255)
                                                   aCCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAt                                                                                               <  1:1550306/150‑1 (MQ=255)
                                                     cTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAATGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtgt                                                                                             >  2:4203515/1‑150 (MQ=255)
                                                     cTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtgt                                                                                             >  2:4202977/1‑150 (MQ=255)
                                                     cTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtgt                                                                                             >  2:3268588/1‑150 (MQ=255)
                                                      tGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtgtg                                                                                            <  1:1721501/150‑1 (MQ=255)
                                                        cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaa                                                                                          >  2:493824/1‑150 (MQ=255)
                                                        cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaa                                                                                          >  2:2857328/1‑150 (MQ=255)
                                                            cGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAttt                                                                                      >  2:4271833/1‑150 (MQ=255)
                                                             gggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAtttt                                                                                     >  2:1341297/1‑150 (MQ=255)
                                                             gggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAtttt                                                                                     <  1:1288820/150‑1 (MQ=255)
                                                                 gTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCt                                                                                 <  1:190815/150‑1 (MQ=255)
                                                                   tgtgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTgg                                                                               <  2:895861/150‑1 (MQ=255)
                                                                     tgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCa                                                                             <  2:219861/150‑1 (MQ=255)
                                                                       cTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGt                                                                           >  1:1154918/1‑150 (MQ=255)
                                                                            ggTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCca                                                                      >  1:240970/1‑150 (MQ=255)
                                                                            ggTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCca                                                                      <  2:168887/150‑1 (MQ=255)
                                                                            ggTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCca                                                                      <  2:4316133/150‑1 (MQ=255)
                                                                              ttCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCcaca                                                                    <  2:1968080/150‑1 (MQ=255)
                                                                               tCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaa                                                                   >  2:12/1‑150 (MQ=255)
                                                                               tCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaa                                                                   >  2:779320/1‑150 (MQ=255)
                                                                                ccTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaac                                                                  <  1:2909750/150‑1 (MQ=255)
                                                                                       tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           <  1:1761969/150‑1 (MQ=255)
                                                                                       tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           <  1:3838133/150‑1 (MQ=255)
                                                                                       tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           <  1:1761743/150‑1 (MQ=255)
                                                                                         tAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCggg                                                         >  2:3456878/1‑150 (MQ=255)
                                                                                         tAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTATGGCAGTCGCCACAACATTACTggg                                                         >  2:4297523/1‑150 (MQ=255)
                                                                                         tAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGAGGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCggg                                                         <  2:3005621/150‑1 (MQ=255)
                                                                                           agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       >  1:227846/1‑150 (MQ=255)
                                                                                           agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       >  1:227534/1‑150 (MQ=255)
                                                                                           agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       <  2:3063057/150‑1 (MQ=255)
                                                                                             agagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAt                                                     >  1:3580571/1‑150 (MQ=255)
                                                                                             agagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAt                                                     >  1:3580619/1‑150 (MQ=255)
                                                                                             agagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAt                                                     <  1:3819520/150‑1 (MQ=255)
                                                                                              gagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAtt                                                    <  1:1248456/150‑1 (MQ=255)
                                                                                                 acctcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 <  1:3054886/147‑1 (MQ=255)
                                                                                                 aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 >  1:1964081/1‑150 (MQ=255)
                                                                                                 aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 <  2:3630236/150‑1 (MQ=255)
                                                                                                  catcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACa                                                >  2:3869758/1‑150 (MQ=255)
                                                                                                   atcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACat                                               >  2:4506947/1‑150 (MQ=255)
                                                                                                     catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             <  2:2584787/150‑1 (MQ=255)
                                                                                                     catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             <  1:1807471/150‑1 (MQ=255)
                                                                                                     catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             <  1:4224003/150‑1 (MQ=255)
                                                                                                     catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             >  1:2206138/1‑150 (MQ=255)
                                                                                                     catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             >  2:1399107/1‑150 (MQ=255)
                                                                                                        gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          >  2:4553676/1‑150 (MQ=255)
                                                                                                        gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          >  2:3298055/1‑150 (MQ=255)
                                                                                                           tttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGa                                       <  1:3157632/150‑1 (MQ=255)
                                                                                                            ttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGAt                                      <  1:4555723/150‑1 (MQ=255)
                                                                                                            ttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGAt                                      >  2:2136074/1‑150 (MQ=255)
                                                                                                             tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATc                                     >  2:599078/1‑150 (MQ=255)
                                                                                                                tCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTgg                                  <  1:4134133/150‑1 (MQ=255)
                                                                                                                          cTGCCATCGCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACg                        <  2:2655991/150‑1 (MQ=255)
                                                                                                                              cATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGAc                    >  1:3020956/1‑150 (MQ=255)
                                                                                                                               aTCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACa                   >  1:563402/1‑150 (MQ=255)
                                                                                                                               aTCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACa                   >  2:440361/1‑150 (MQ=255)
                                                                                                                                 ctcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAc                 >  1:2304335/1‑150 (MQ=255)
                                                                                                                                  tcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcg                >  1:4084755/1‑150 (MQ=255)
                                                                                                                                      tGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcgt            >  1:2777104/1‑150 (MQ=255)
                                                                                                                                        gTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAt          <  2:1086647/150‑1 (MQ=255)
                                                                                                                                             gATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCt     >  1:1424573/1‑150 (MQ=255)
                                                                                                                                             gATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCt     <  1:1986101/150‑1 (MQ=255)
                                                                                                                                              attattAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCtt    <  2:1887215/150‑1 (MQ=255)
                                                                                                                                                tATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgt  >  2:4253619/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGT  >  CP009273/1152131‑1152424

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: