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,2780GT99.5% 709.0 / ‑3.8 200E16* (GAA→TAA) ycfHputative DNase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (110/89);  total (110/90)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

CATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGT  >  CP009273/1152132‑1152424
                                                                                                                                                  |                                                                                                                                                  
cATCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAAt                                                                                                                                                 >  1:2366047/1‑150 (MQ=255)
 aTCAACCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACCTCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATc                                                                                                                                                <  2:2143067/150‑1 (MQ=255)
    aaCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCttt                                                                                                                                             >  2:2126773/1‑150 (MQ=255)
    aaCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCttt                                                                                                                                             >  2:1946541/1‑150 (MQ=255)
     aCCGCGAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCGCGATGGTCTGGATTATTAATCTTTg                                                                                                                                            >  1:600879/1‑150 (MQ=255)
         cgAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATa                                                                                                                                        <  2:724762/150‑1 (MQ=255)
          gAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATaa                                                                                                                                       >  1:3549945/1‑150 (MQ=255)
          gAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATaa                                                                                                                                       >  1:3455157/1‑150 (MQ=255)
          gAGCTTCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATaa                                                                                                                                       >  2:2341570/1‑150 (MQ=255)
              ttCTCATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGAc                                                                                                                                   <  2:3494748/150‑1 (MQ=255)
               tctcATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGTGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACg                                                                                                                                  >  2:3923684/1‑150 (MQ=255)
                ctcATCACCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTGTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGt                                                                                                                                 >  2:3981296/1‑150 (MQ=255)
                 tcatcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGTGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTg                                                                                                                                >  2:243055/1‑150 (MQ=255)
                 tcatcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTg                                                                                                                                <  2:2735815/150‑1 (MQ=255)
                  catcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTgg                                                                                                                               <  1:2070790/150‑1 (MQ=255)
                   atcaCCGATCTTTTACTGCGTATTTAGCATTACCTGCAACCGGTAGATGAGCTACAGTATCCTCATCAGTAAGACAGACATCATGTTTTTATTCGACTCACAATGACAACTCGATGGTCTGGCTTATTAATAATTGCATAAGGACGTGGa                                                                                                                              >  1:560852/1‑150 (MQ=255)
                    tcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGAt                                                                                                                             >  1:3749396/1‑150 (MQ=255)
                    tcaCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGAt                                                                                                                             >  1:1651549/1‑150 (MQ=255)
                      aCCGATCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGa                                                                                                                           >  2:1454888/1‑150 (MQ=255)
                          aTCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGtt                                                                                                                       <  1:2159305/150‑1 (MQ=255)
                          aTCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGtt                                                                                                                       <  1:1293373/150‑1 (MQ=255)
                           tCTTTTACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTc                                                                                                                      <  1:4067211/150‑1 (MQ=255)
                             ttttACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTg                                                                                                                    <  1:4073433/150‑1 (MQ=255)
                             ttttACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTg                                                                                                                    >  1:67011/1‑150 (MQ=255)
                             ttttACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTg                                                                                                                    >  1:1347162/1‑150 (MQ=255)
                              tttACTGCGTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTgg                                                                                                                   >  2:237923/1‑150 (MQ=255)
                                      gTATTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgcc                                                                                                           >  2:1223319/1‑150 (MQ=255)
                                        aTTGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgccgc                                                                                                         >  2:3004839/1‑150 (MQ=255)
                                         ttGAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAgccgcc                                                                                                        <  2:45665/150‑1 (MQ=255)
                                           gAGCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCcgccgc                                                                                                      >  2:628809/1‑150 (MQ=255)
                                             gCATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAc                                                                                                    >  2:1651579/1‑150 (MQ=255)
                                              cATTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcg                                                                                                   <  2:174137/150‑1 (MQ=255)
                                               aTTACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgc                                                                                                  >  1:2847413/1‑150 (MQ=255)
                                                ttACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCAcgcg                                                                                                 <  2:3900540/150‑1 (MQ=255)
                                                 tACCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGa                                                                                                >  2:3750778/1‑150 (MQ=255)
                                                  aCCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAt                                                                                               >  1:931259/1‑150 (MQ=255)
                                                  aCCTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAt                                                                                               <  1:777887/150‑1 (MQ=255)
                                                   ccTGCAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGAtg                                                                                              >  2:96000/1‑150 (MQ=255)
                                                       cAACCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaa                                                                                          >  1:1479641/1‑150 (MQ=255)
                                                        aaCCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGaaa                                                                                         >  2:2883581/1‑150 (MQ=255)
                                                         aCCGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAt                                                                                        <  2:175376/150‑1 (MQ=255)
                                                          ccGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAtt                                                                                       >  2:3064325/1‑150 (MQ=255)
                                                          ccGGGCGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAAtt                                                                                       >  2:3734396/1‑150 (MQ=255)
                                                               cGTTGTGCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTc                                                                                  >  2:1016326/1‑150 (MQ=255)
                                                                   gtgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGc                                                                              >  2:1546308/1‑150 (MQ=255)
                                                                    tgCTACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCa                                                                             <  2:1558258/150‑1 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGATTCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          <  1:2641401/150‑1 (MQ=255)
                                                                       tACCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTc                                                                          >  1:955466/1‑150 (MQ=255)
                                                                        aCCGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCg                                                                         <  2:3605059/150‑1 (MQ=255)
                                                                          cGGTTCCTCATCTTTCAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCNCGCGATGTGAAATTTTGTCTGGCAGTCGcc                                                                       <  1:3855385/150‑1 (MQ=255)
                                                                          cGGTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGcc                                                                       >  1:3569795/1‑150 (MQ=255)
                                                                           ggTTCCTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCca                                                                      >  1:989525/1‑150 (MQ=255)
                                                                               ccTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaac                                                                  >  2:2775541/1‑150 (MQ=255)
                                                                                cTCATCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCacaaca                                                                 <  1:3621897/150‑1 (MQ=255)
                                                                                   aTCTTTCAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTa                                                              <  2:2479774/150‑1 (MQ=255)
                                                                                    tCTTTAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTAc                                                             <  2:3885238/150‑1 (MQ=255)
                                                                                      tttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCg                                                           >  1:708675/1‑150 (MQ=255)
                                                                                       ttAAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCgg                                                          >  1:4166907/1‑150 (MQ=255)
                                                                                          agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       >  2:2809663/1‑150 (MQ=255)
                                                                                          agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGtt                                                       <  1:3272452/150‑1 (MQ=255)
                                                                                           gagagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTa                                                      >  1:1982332/1‑150 (MQ=255)
                                                                                           gagagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTa                                                      >  2:516049/1‑150 (MQ=255)
                                                                                                aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 >  2:3718176/1‑150 (MQ=255)
                                                                                                aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 <  1:561836/150‑1 (MQ=255)
                                                                                                aCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTAc                                                 <  2:146707/150‑1 (MQ=255)
                                                                                                 catcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACa                                                >  1:846145/1‑150 (MQ=255)
                                                                                                    catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACatat                                             <  1:1727464/150‑1 (MQ=255)
                                                                                                     atGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATg                                            <  1:650619/150‑1 (MQ=255)
                                                                                                      tGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGc                                           <  1:827682/150‑1 (MQ=255)
                                                                                                       gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          <  1:4156365/150‑1 (MQ=255)
                                                                                                       gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCg                                          <  2:3110783/150‑1 (MQ=255)
                                                                                                        tttttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCgg                                         <  1:1964325/150‑1 (MQ=255)
                                                                                                         ttttAGTCGACTCACACTGCCATCTCGATGTTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTTTCTGGCAGTCGCCACAACATTACCGGGTTATTTAAATATGAggg                                        >  2:1026698/1‑150 (MQ=255)
                                                                                                               tCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTgg                                  >  2:2871608/1‑150 (MQ=255)
                                                                                                               tCGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTgg                                  >  2:544421/1‑150 (MQ=255)
                                                                                                                cGACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGt                                 <  2:1330130/150‑1 (MQ=255)
                                                                                                                 gACTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTa                                >  2:2126576/1‑150 (MQ=255)
                                                                                                                  aCTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAg                               <  2:2461633/150‑1 (MQ=255)
                                                                                                                   cTCACACTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAgg                              >  1:3374221/1‑150 (MQ=255)
                                                                                                                       cacTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGaa                          >  1:3710635/1‑150 (MQ=255)
                                                                                                                         cTGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACg                        <  2:3391747/150‑1 (MQ=255)
                                                                                                                          tGCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCAGGTAGGCGAACGt                       <  2:824577/150‑1 (MQ=255)
                                                                                                                           gCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGTCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTg                      <  2:1798960/150‑1 (MQ=255)
                                                                                                                           gCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTg                      <  2:90551/150‑1 (MQ=255)
                                                                                                                           gCCATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTg                      >  2:1959712/1‑150 (MQ=255)
                                                                                                                             cATCTCGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGAc                    >  1:103251/1‑150 (MQ=255)
                                                                                                                                ctcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAc                 <  2:4050156/150‑1 (MQ=255)
                                                                                                                                ctcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAc                 <  2:1178900/150‑1 (MQ=255)
                                                                                                                                 tcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcg                >  1:1916826/1‑150 (MQ=255)
                                                                                                                                 tcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcg                >  2:1701182/1‑150 (MQ=255)
                                                                                                                                 tcGATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcg                <  2:3549945/150‑1 (MQ=255)
                                                                                                                                   gATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtc              >  1:4013246/1‑150 (MQ=255)
                                                                                                                                   gATGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtc              <  1:4102701/150‑1 (MQ=255)
                                                                                                                                    aTGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcg             >  1:3000485/1‑150 (MQ=255)
                                                                                                                                     tGGTCTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcgt            <  2:1349780/150‑1 (MQ=255)
                                                                                                                                         cTGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAttt        <  1:2511596/150‑1 (MQ=255)
                                                                                                                                          tGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt       <  2:1994009/150‑1 (MQ=255)
                                                                                                                                          tGGATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt       <  2:4096738/150‑1 (MQ=255)
                                                                                                                                           ggATTATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTc      >  1:1540656/1‑150 (MQ=255)
                                                                                                                                               tATTAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACACCATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgt  <  1:27823/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: ‑