Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,321,444 T→A 6.5% D60V (GAT→GTT)  yfaP ← DUF2135 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,321,4440TA 6.5% 980.2 / 39.5 285D60V (GAT→GTT) yfaPDUF2135 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base T (154/112);  new base A (11/7);  total (166/119)
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.99e-01

Only 100 of 303 total aligned reads displayed.

CAGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGCCACC  >  CP009273/2321295‑2321593
                                                                                                                                                     |                                                                                                                                                     
cAGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGat                                                                                                                                                       <  2:3852492/150‑1 (MQ=255)
 aGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatc                                                                                                                                                      >  2:5124077/1‑150 (MQ=255)
 aGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCATTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatc                                                                                                                                                      >  2:1654527/1‑150 (MQ=255)
  gTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatca                                                                                                                                                     >  1:3580646/1‑150 (MQ=255)
   tGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGatcat                                                                                                                                                    <  1:228896/150‑1 (MQ=255)
    gAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATCATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGaacatc                                                                                                                                                   >  1:3515952/1‑150 (MQ=255)
     aaaaTATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCt                                                                                                                                                  <  2:217354/150‑1 (MQ=255)
        atatAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGcc                                                                                                                                               <  1:4696717/150‑1 (MQ=255)
        atatAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGcc                                                                                                                                               >  2:2375618/1‑150 (MQ=255)
         tataAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCa                                                                                                                                              >  2:2966690/1‑150 (MQ=255)
          ataAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCAt                                                                                                                                             <  2:1207245/150‑1 (MQ=255)
           taAGGGCGTGCACATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGACAATCTGCCATg                                                                                                                                            >  2:4931789/1‑150 (MQ=255)
           taAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATCTTTTGATCATCTGCCATg                                                                                                                                            >  1:2895352/1‑150 (MQ=255)
           taAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATg                                                                                                                                            <  1:3376858/150‑1 (MQ=255)
             aGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGtt                                                                                                                                          <  1:4489397/150‑1 (MQ=255)
             aGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGtt                                                                                                                                          >  2:411511/1‑150 (MQ=255)
             aGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGtt                                                                                                                                          >  2:1450733/1‑150 (MQ=255)
              gggCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTg                                                                                                                                         <  1:3645061/150‑1 (MQ=255)
                gCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgac                                                                                                                                       >  1:1936837/1‑150 (MQ=255)
                 cGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacg                                                                                                                                      <  2:817574/150‑1 (MQ=255)
                 cGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTgacg                                                                                                                                      <  1:1159070/150‑1 (MQ=255)
                 cGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTgacg                                                                                                                                      <  1:321317/150‑1 (MQ=255)
                  gTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacgg                                                                                                                                     >  2:299319/1‑150 (MQ=255)
                    gaaaaTGAAACATCTGATTCAATTCTCTGCGGCATACTGTCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGACTCATCATAATTATTTATTTTGCCGCGGATCTGCGCTGATATTTTTTGATAATCTGCCATGTTaacggac                                                                                                                                   <  2:543767/148‑1 (MQ=255)
                    gCAAATGAACCATCGGATTCAATACTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTgacggac                                                                                                                                   <  2:4854285/150‑1 (MQ=255)
                        aTGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGacgacg                                                                                                                               >  1:4194924/1‑150 (MQ=255)
                            aCCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCgg                                                                                                                           >  1:2861450/1‑150 (MQ=255)
                             ccATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACAACGCCggg                                                                                                                          <  2:2323652/150‑1 (MQ=255)
                               aTCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTa                                                                                                                        <  2:2562980/150‑1 (MQ=255)
                                 cGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAAt                                                                                                                      >  1:873606/1‑150 (MQ=255)
                                   gATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTg                                                                                                                    >  1:4239786/1‑150 (MQ=255)
                                   gATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTg                                                                                                                    <  1:51203/150‑1 (MQ=255)
                                    aTTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGa                                                                                                                   >  1:3788018/1‑150 (MQ=255)
                                     ttCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGAt                                                                                                                  >  2:3062757/1‑150 (MQ=255)
                                      tCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATc                                                                                                                 <  2:1415410/150‑1 (MQ=255)
                                      tCAAGTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATc                                                                                                                 <  1:2433302/150‑1 (MQ=255)
                                         aTTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGAt                                                                                                              >  1:1056779/1‑150 (MQ=255)
                                          ttCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGAtt                                                                                                             <  2:1272194/150‑1 (MQ=255)
                                          ttCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGAtt                                                                                                             <  2:4094585/150‑1 (MQ=255)
                                           tctcTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATtg                                                                                                            <  2:738548/150‑1 (MQ=255)
                                           tctcTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATtg                                                                                                            >  1:584218/1‑150 (MQ=255)
                                           tctcTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTGACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATtg                                                                                                            <  2:5154994/150‑1 (MQ=255)
                                            ctctGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATtgt                                                                                                           >  1:3121144/1‑150 (MQ=255)
                                              ctgctgCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGa                                                                                                         <  2:2717105/150‑1 (MQ=255)
                                                gCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGaaa                                                                                                       <  1:4629209/150‑1 (MQ=255)
                                                gCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGaaa                                                                                                       <  1:5178098/150‑1 (MQ=255)
                                                 cGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTACTTGATCGATTGTGaaaa                                                                                                      >  2:4751037/1‑150 (MQ=255)
                                                     aTACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTg                                                                                                  <  1:3983505/150‑1 (MQ=255)
                                                       aCTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGc                                                                                                >  1:3867454/1‑150 (MQ=255)
                                                       aCTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGc                                                                                                <  2:340835/150‑1 (MQ=255)
                                                       aCTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGc                                                                                                <  1:2017875/150‑1 (MQ=255)
                                                       aCTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGc                                                                                                <  1:2589899/150‑1 (MQ=255)
                                                           gCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCt                                                                                            <  1:2774854/150‑1 (MQ=255)
                                                           gCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCt                                                                                            >  1:2667652/1‑150 (MQ=255)
                                                             aCCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTc                                                                                          <  2:1789724/150‑1 (MQ=255)
                                                               cGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGc                                                                                        >  2:3482479/1‑150 (MQ=255)
                                                                gTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGcc                                                                                       <  2:431473/150‑1 (MQ=255)
                                                                  tGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCct                                                                                     >  2:2996455/1‑150 (MQ=255)
                                                                       accaGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGcc                                                                                >  1:4867549/1‑150 (MQ=255)
                                                                           gCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGAt                                                                            <  1:3438392/150‑1 (MQ=255)
                                                                              ggCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCGCTGCCGGATGcc                                                                         >  2:3251201/1‑150 (MQ=255)
                                                                                 ccTGCTGAACTTTTACCGCCGCAGCATAATACTTTATTTTGCAGCGGATATGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGTTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGc                                                                      <  2:112926/150‑1 (MQ=255)
                                                                                 ccTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGc                                                                      <  2:4380167/150‑1 (MQ=255)
                                                                                  ctgctgAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGcc                                                                     >  1:125462/1‑150 (MQ=255)
                                                                                   tgctgAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCt                                                                    >  2:626563/1‑150 (MQ=255)
                                                                                      tgAACTTTACCCGCCGCAGCATAATTCCTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTgag                                                                 <  1:2739021/150‑1 (MQ=255)
                                                                                        aaCTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGagag                                                               >  1:4140329/1‑150 (MQ=255)
                                                                                            ttACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGt                                                           <  1:941577/150‑1 (MQ=255)
                                                                                             tcccggccgcAGCAATCTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTg                                                          <  2:2686653/145‑1 (MQ=255)
                                                                                                ccgccgCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCAt                                                       <  1:3487589/150‑1 (MQ=255)
                                                                                                 cgccgcAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATc                                                      >  1:4013118/1‑150 (MQ=255)
                                                                                                    cgcAGCATAATTCTTTATTTTGCCGCGGATCTGCGCAGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCAtactg                                                   >  2:3705363/1‑146 (MQ=255)
                                                                                                        gcaTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCg                                               <  2:1056779/150‑1 (MQ=255)
                                                                                                          aTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGAc                                             >  1:2982927/1‑150 (MQ=255)
                                                                                                          aTAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGAc                                             >  1:1778603/1‑150 (MQ=255)
                                                                                                           tAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACt                                            <  1:5076646/150‑1 (MQ=255)
                                                                                                             aTTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTc                                          >  1:4249932/1‑150 (MQ=255)
                                                                                                               tCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCt                                        >  2:4512816/1‑150 (MQ=255)
                                                                                                                    aTTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGcaca                                   >  1:779118/1‑150 (MQ=255)
                                                                                                                      tttGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACAcc                                 >  1:4417829/1‑150 (MQ=255)
                                                                                                                        tGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCtt                               >  1:1088584/1‑150 (MQ=255)
                                                                                                                        tGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCtt                               <  1:534592/150‑1 (MQ=255)
                                                                                                                           cgcgGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCAc                            >  1:1644520/1‑150 (MQ=255)
                                                                                                                             cgGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtg                          >  2:1337653/1‑150 (MQ=255)
                                                                                                                              ggATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgt                         >  2:1036942/1‑150 (MQ=255)
                                                                                                                               gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACtgtg                        >  1:3054502/1‑150 (MQ=255)
                                                                                                                               gATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCACGACTTCCTGCACACCTTCACtgtg                        >  1:1478181/1‑150 (MQ=255)
                                                                                                                                 tCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAg                      <  1:247970/150‑1 (MQ=255)
                                                                                                                                 tCTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTTAg                      >  2:2481280/1‑150 (MQ=255)
                                                                                                                                 ncTGCGCTGAGATATTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAg                      >  2:2261603/2‑150 (MQ=255)
                                                                                                                                  cTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGc                     >  2:4616653/1‑150 (MQ=255)
                                                                                                                                    gcgcTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGa                   <  2:933770/150‑1 (MQ=255)
                                                                                                                                     cgcTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGAc                  >  1:3342647/1‑150 (MQ=255)
                                                                                                                                        tGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAgg               >  2:2107609/1‑150 (MQ=255)
                                                                                                                                         gagaTATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGc              <  1:1303811/150‑1 (MQ=255)
                                                                                                                                            atatTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGaa           >  1:419973/1‑150 (MQ=255)
                                                                                                                                            atatTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGaa           >  1:1150691/1‑150 (MQ=255)
                                                                                                                                            atatTTTGAACATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGaa           >  2:5015951/1‑150 (MQ=255)
                                                                                                                                               ttttGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAgc        >  1:3174525/1‑150 (MQ=255)
                                                                                                                                                   gATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGcca    <  2:1741741/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTGAAAATATAAGGGCGTGCAAATGAACCATCGGATTCAATTCTCTGCGGCATACTGGCACCGTTGACCACCAGCCGGCCCTGCTGAACTTTACCCGCCGCAGCATAATTCTTTATTTTGCCGCGGATCTGCGCTGAGATATTTTGATCATCTGCCATGTTGACGGACGACGCCGGGTAATTGATCGATTGTGAAAAACTGGCATCTTCGCCCTCTGCCGGATGCCAGCCTGAGAGCGGTGCATCAATCTCGACTTCCTGCACACCTTCACTGTGAGCGACAGGCGAAAGCGCCACC  >  CP009273/2321295‑2321593

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: