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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,152,2780GT96.1% 841.1 / 22.7 259E16* (GAA→TAA) ycfHputative DNase
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base T (124/125);  total (128/131)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 261 total aligned reads displayed.

CATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGT  >  CP009273/1152132‑1152424
                                                                                                                                                  |                                                                                                                                                  
   cAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCtt                                                                                                                                              >  2:3895610/1‑150 (MQ=255)
     aCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTg                                                                                                                                            <  1:1294445/150‑1 (MQ=255)
     aCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTg                                                                                                                                            <  2:473066/150‑1 (MQ=255)
         cgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATa                                                                                                                                        <  1:2864911/150‑1 (MQ=255)
         cgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATa                                                                                                                                        >  2:277060/1‑150 (MQ=255)
                     caCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATg                                                                                                                            >  1:3515217/1‑150 (MQ=255)
                       ccGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGAc                                                                                                                          <  1:2600708/150‑1 (MQ=255)
                        cGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACg                                                                                                                         <  2:2092119/150‑1 (MQ=255)
                         gATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGt                                                                                                                        >  1:2286350/1‑150 (MQ=255)
                            cTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCt                                                                                                                     <  2:2488876/150‑1 (MQ=255)
                                 aCTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGa                                                                                                                <  2:834711/150‑1 (MQ=255)
                                  cTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGaa                                                                                                               <  2:1984595/150‑1 (MQ=255)
                                    gCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAg                                                                                                             <  2:1330207/150‑1 (MQ=255)
                                    gCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAg                                                                                                             <  2:3653200/150‑1 (MQ=255)
                                      gTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgcc                                                                                                           >  2:2547414/1‑150 (MQ=255)
                                            aGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCa                                                                                                     >  2:2202267/1‑150 (MQ=255)
                                            aGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCa                                                                                                     <  1:3860652/150‑1 (MQ=255)
                                             gCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAc                                                                                                    >  1:4338599/1‑150 (MQ=255)
                                               aTTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgc                                                                                                  >  1:20699/1‑150 (MQ=255)
                                               aTTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgc                                                                                                  >  2:4055056/1‑150 (MQ=255)
                                               aTTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgc                                                                                                  >  2:2344504/1‑150 (MQ=255)
                                                ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgcg                                                                                                 >  1:2693061/1‑150 (MQ=255)
                                                 tACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGa                                                                                                >  2:3796773/1‑150 (MQ=255)
                                                  aCCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAt                                                                                               <  2:3026111/150‑1 (MQ=255)
                                                   ccTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtg                                                                                              >  1:4429744/1‑150 (MQ=255)
                                                    cTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtgt                                                                                             >  1:3572574/1‑150 (MQ=255)
                                                       cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaa                                                                                          >  1:2919686/1‑150 (MQ=255)
                                                       cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaa                                                                                          >  1:1047787/1‑150 (MQ=255)
                                                       cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGAATAAGGACGTGGATGACGTTCTGGCGAAATCCGCCGCACGCGATGTGaa                                                                                          >  2:1847171/1‑150 (MQ=255)
                                                           cGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAttt                                                                                      >  2:541497/1‑150 (MQ=255)
                                                            gggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAtttt                                                                                     <  2:2876682/150‑1 (MQ=255)
                                                             ggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTg                                                                                    <  2:3320267/150‑1 (MQ=255)
                                                             ggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTg                                                                                    <  2:3624628/150‑1 (MQ=255)
                                                             ggCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTg                                                                                    <  2:820326/150‑1 (MQ=255)
                                                                gTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCt                                                                                 >  2:3610230/1‑150 (MQ=255)
                                                                 ttGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTg                                                                                <  2:1488208/150‑1 (MQ=255)
                                                                   gtgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGc                                                                              <  1:3389303/150‑1 (MQ=255)
                                                                   gtgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGc                                                                              >  2:1084480/1‑150 (MQ=255)
                                                                    tgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCa                                                                             <  2:2242959/150‑1 (MQ=255)
                                                                    tgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCa                                                                             <  2:211076/150‑1 (MQ=255)
                                                                    tgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCa                                                                             <  1:749929/150‑1 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          <  1:2760120/150‑1 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          >  1:4019794/1‑150 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          >  1:2614266/1‑150 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          <  2:3552639/150‑1 (MQ=255)
                                                                        aCCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCg                                                                         <  1:1071419/150‑1 (MQ=255)
                                                                        aCCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCg                                                                         >  1:4140691/1‑150 (MQ=255)
                                                                          cGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGcc                                                                       >  1:1415706/1‑150 (MQ=255)
                                                                           ggTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCca                                                                      <  2:206836/150‑1 (MQ=255)
                                                                            gTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCcac                                                                     >  2:2301069/1‑150 (MQ=255)
                                                                            gTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCcac                                                                     <  2:4251778/150‑1 (MQ=255)
                                                                            gTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCcac                                                                     <  2:1825558/150‑1 (MQ=255)
                                                                              tCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaa                                                                   >  1:3247403/1‑150 (MQ=255)
                                                                                 tCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACAt                                                                <  1:3129628/150‑1 (MQ=255)
                                                                                 tCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACAt                                                                <  1:3833867/150‑1 (MQ=255)
                                                                                    tCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTAc                                                             <  2:3059773/150‑1 (MQ=255)
                                                                                      tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           <  1:1484681/150‑1 (MQ=255)
                                                                                      tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           >  1:2916100/1‑150 (MQ=255)
                                                                                      tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           >  1:3673285/1‑150 (MQ=255)
                                                                                       ttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCgg                                                          >  1:2354417/1‑150 (MQ=255)
                                                                                       ttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCgg                                                          <  2:617716/150‑1 (MQ=255)
                                                                                        tAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCggg                                                         >  2:1341297/1‑150 (MQ=255)
                                                                                          agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       <  1:4466507/150‑1 (MQ=255)
                                                                                          agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       <  2:1391087/150‑1 (MQ=255)
                                                                                            agagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAt                                                     >  1:4492713/1‑150 (MQ=255)
                                                                                             gagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAtt                                                    <  1:1777296/150‑1 (MQ=255)
                                                                                             gagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAtt                                                    <  1:872000/150‑1 (MQ=255)
                                                                                              agaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTAttt                                                   <  1:54450/150‑1 (MQ=255)
                                                                                               gaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTa                                                  >  1:1439542/1‑150 (MQ=255)
                                                                                                aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 >  1:2766861/1‑150 (MQ=255)
                                                                                                 catcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACa                                                <  1:2227743/150‑1 (MQ=255)
                                                                                                  atcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACat                                               <  1:3955247/150‑1 (MQ=255)
                                                                                                    catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             <  1:569364/150‑1 (MQ=255)
                                                                                                     atGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCAAAAGCCGCCTCACGCGTTGTGAAATTTTGTCTGTCAGTCACCACAACATTACCGGGTTATTTACATATg                                            <  2:4148371/150‑1 (MQ=255)
                                                                                                      tGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGc                                           >  2:2860894/1‑150 (MQ=255)
                                                                                                       gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          >  2:566762/1‑150 (MQ=255)
                                                                                                       gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          <  2:1000298/150‑1 (MQ=255)
                                                                                                         ttttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCggg                                        <  2:3376864/150‑1 (MQ=255)
                                                                                                           ttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGAt                                      >  1:1110931/1‑150 (MQ=255)
                                                                                                            tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTACTCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATc                                     <  1:1933034/150‑1 (MQ=255)
                                                                                                            tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATc                                     <  1:1279562/150‑1 (MQ=255)
                                                                                                              gTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTg                                   >  1:3348250/1‑150 (MQ=255)
                                                                                                              gTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTg                                   <  2:1654485/150‑1 (MQ=255)
                                                                                                              gTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTg                                   <  2:1654502/150‑1 (MQ=255)
                                                                                                               tCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTgg                                  >  1:1502573/1‑150 (MQ=255)
                                                                                                                cGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGt                                 >  2:3737826/1‑150 (MQ=255)
                                                                                                                  aCTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAg                               <  1:4262276/150‑1 (MQ=255)
                                                                                                                  aCTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAg                               <  2:1290514/150‑1 (MQ=255)
                                                                                                                   cTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAgg                              >  2:4328455/1‑150 (MQ=255)
                                                                                                                     cacaCTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCg                            >  1:3046196/1‑150 (MQ=255)
                                                                                                                     cacaCTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCg                            <  1:2416888/150‑1 (MQ=255)
                                                                                                                          tGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGt                       <  2:2614266/150‑1 (MQ=255)
                                                                                                                           gCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTg                      <  1:777492/150‑1 (MQ=255)
                                                                                                                              aTCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACa                   >  2:2963984/1‑150 (MQ=255)
                                                                                                                              aTCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACa                   <  2:922070/150‑1 (MQ=255)
                                                                                                                                    aTGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcg             >  1:1714507/1‑150 (MQ=255)
                                                                                                                                    aTGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcg             >  2:2805535/1‑150 (MQ=255)
                                                                                                                                     tGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcgt            >  2:3501348/1‑150 (MQ=255)
                                                                                                                                          tGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt       >  2:4465867/1‑150 (MQ=255)
                                                                                                                                            gATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCt     <  1:920111/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGT  >  CP009273/1152132‑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: