Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,058,961 C→A 5.4% A91S (GCC→TCC)  cbpA ← DnaK co‑chaperone; curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,058,9610CA 5.4% 1289.9 / 38.0 331A91S (GCC→TCC) cbpADnaK co‑chaperone; curved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (142/171);  new base A (11/7);  total (153/178)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 332 total aligned reads displayed.

TGATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCA  >  CP009273/1058814‑1059109
                                                                                                                                                   |                                                                                                                                                    
tGATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCg                                                                                                                                                     <  2:567940/149‑1 (MQ=255)
tGATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCg                                                                                                                                                     <  2:2351927/149‑1 (MQ=255)
 gATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGt                                                                                                                                                    <  2:2754151/149‑1 (MQ=255)
 gATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGt                                                                                                                                                    <  1:444541/149‑1 (MQ=255)
 gATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGAGt                                                                                                                                                    >  1:785649/1‑149 (MQ=255)
   tGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTa                                                                                                                                                  >  2:2816518/1‑149 (MQ=255)
     gTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTaaa                                                                                                                                                >  2:4722104/1‑149 (MQ=255)
      tACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTaaaa                                                                                                                                               <  1:435466/149‑1 (MQ=255)
      tACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTaaaa                                                                                                                                               >  2:1755017/1‑149 (MQ=255)
       aCGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAAc                                                                                                                                              <  2:1932310/149‑1 (MQ=255)
          cTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTTACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAActct                                                                                                                                           >  2:1440863/1‑149 (MQ=255)
          cTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAActct                                                                                                                                           <  1:1857137/149‑1 (MQ=255)
          cTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAActct                                                                                                                                           >  2:1156541/1‑149 (MQ=255)
          cTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAActct                                                                                                                                           >  1:399708/1‑149 (MQ=255)
           ttATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTg                                                                                                                                          >  1:2107875/1‑149 (MQ=255)
            tATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGa                                                                                                                                         <  2:3846838/149‑1 (MQ=255)
             aTGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGAc                                                                                                                                        <  1:4234784/149‑1 (MQ=255)
                   aGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGc                                                                                                                                  <  1:4655071/149‑1 (MQ=255)
                   aGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGc                                                                                                                                  >  1:3394277/1‑149 (MQ=255)
                     tAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCa                                                                                                                                >  2:206633/1‑149 (MQ=255)
                     tAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCa                                                                                                                                <  1:4554845/149‑1 (MQ=255)
                      aaGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCAt                                                                                                                               <  2:4196220/149‑1 (MQ=255)
                        gCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCAtgg                                                                                                                             <  1:3397652/149‑1 (MQ=255)
                           tttCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCAtggtgg                                                                                                                          >  2:1183360/1‑149 (MQ=255)
                            ttcttcGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGa                                                                                                                         >  2:2608016/1‑149 (MQ=255)
                             tcttcGAGGAATCCCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGaa                                                                                                                        <  1:2624043/149‑1 (MQ=255)
                             tcttcGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGaa                                                                                                                        <  2:4101812/149‑1 (MQ=255)
                               ttcGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACt                                                                                                                      <  2:79055/149‑1 (MQ=255)
                                 cGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGa                                                                                                                    <  1:3570476/149‑1 (MQ=255)
                                   aGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACg                                                                                                                  >  2:3090807/1‑149 (MQ=255)
                                     gAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACggg                                                                                                                >  2:1064100/1‑148 (MQ=255)
                                      aaTACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCACGGTGGAACTGACGGtt                                                                                                               >  2:1489966/1‑149 (MQ=255)
                                       aTACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTa                                                                                                              >  2:2434636/1‑149 (MQ=255)
                                       aTACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTa                                                                                                              <  1:369499/149‑1 (MQ=255)
                                        tACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTaa                                                                                                             <  2:4844748/149‑1 (MQ=255)
                                          ccGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAAt                                                                                                           <  2:4473691/149‑1 (MQ=255)
                                           cGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAAtt                                                                                                          >  2:4656641/1‑149 (MQ=255)
                                            gCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTg                                                                                                         >  1:2388179/1‑149 (MQ=255)
                                               aCTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCACGGTGGAACTGACGGTTAAATTGCgg                                                                                                      >  1:1553802/1‑149 (MQ=255)
                                                cTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGa                                                                                                     <  2:4587356/149‑1 (MQ=255)
                                                cTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGa                                                                                                     >  2:2123985/1‑149 (MQ=255)
                                                  tCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATc                                                                                                   <  1:4506797/149‑1 (MQ=255)
                                                    gATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTTGCCATGGTGGAACTGACGGTTAAATTGCGGATCGt                                                                                                 <  2:1190357/149‑1 (MQ=255)
                                                     aTTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGtt                                                                                                >  2:2164978/1‑149 (MQ=255)
                                                     aTTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGtt                                                                                                >  1:521120/1‑149 (MQ=255)
                                                      tttCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGAGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTg                                                                                               >  1:4419860/1‑149 (MQ=255)
                                                       ttCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGAGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGc                                                                                              >  1:1769067/1‑149 (MQ=255)
                                                          aaTATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGAt                                                                                           <  2:589741/149‑1 (MQ=255)
                                                            tatCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGt                                                                                         <  2:3191958/149‑1 (MQ=255)
                                                             atCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGtt                                                                                        >  1:654884/1‑149 (MQ=255)
                                                                gTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGcc                                                                                     <  1:1030484/149‑1 (MQ=255)
                                                                gTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAAATGCGGATCGTTGCGATGTTGcc                                                                                     >  1:3339787/1‑149 (MQ=255)
                                                                   gCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCcaca                                                                                  >  2:1689078/1‑149 (MQ=255)
                                                                     cgcgTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATc                                                                                <  2:843046/149‑1 (MQ=255)
                                                                      gcgTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCt                                                                               >  1:2584378/1‑149 (MQ=255)
                                                                          gtgGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATc                                                                           <  1:2727416/149‑1 (MQ=255)
                                                                           tgGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCa                                                                          <  2:3932734/149‑1 (MQ=255)
                                                                           tgGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCa                                                                          <  1:3222879/149‑1 (MQ=255)
                                                                            ggCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCAt                                                                         <  2:4902641/149‑1 (MQ=255)
                                                                            ggCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCAt                                                                         <  2:4892404/149‑1 (MQ=255)
                                                                              cGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATAc                                                                       <  2:3776006/149‑1 (MQ=255)
                                                                               ggggCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACt                                                                      <  2:3888749/149‑1 (MQ=255)
                                                                                  gCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAg                                                                   >  1:5183952/1‑149 (MQ=255)
                                                                                       gACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGa                                                              <  2:5091238/149‑1 (MQ=255)
                                                                                            ctctGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACgttg                                                         <  1:1502624/149‑1 (MQ=255)
                                                                                                 gCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTtcat                                                    <  2:1564041/149‑1 (MQ=255)
                                                                                                 gCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTtcat                                                    <  2:4433953/149‑1 (MQ=255)
                                                                                                  cGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTtcatc                                                   <  2:3321114/149‑1 (MQ=255)
                                                                                                      cATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACtt                                               <  1:2231874/149‑1 (MQ=255)
                                                                                                            gACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACt                                         >  1:883850/1‑149 (MQ=255)
                                                                                                              ccGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGTCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTc                                       <  2:3131015/149‑1 (MQ=255)
                                                                                                              ccGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTc                                       <  2:1351544/149‑1 (MQ=255)
                                                                                                                gAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTccc                                     <  2:4297757/149‑1 (MQ=255)
                                                                                                                 aaaaTTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCa                                    <  1:4862402/149‑1 (MQ=255)
                                                                                                                    aTTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGc                                 <  1:1258423/149‑1 (MQ=255)
                                                                                                                      tGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCtt                               <  2:3411465/149‑1 (MQ=255)
                                                                                                                      tGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCtt                               >  2:2246228/1‑149 (MQ=255)
                                                                                                                        aCGGGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCa                             >  1:254823/1‑149 (MQ=255)
                                                                                                                          gagaAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGc                           <  1:3136474/149‑1 (MQ=255)
                                                                                                                           agaagaTATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCg                          <  1:2164978/149‑1 (MQ=255)
                                                                                                                            gaagaTATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGa                         <  2:3371792/149‑1 (MQ=255)
                                                                                                                            gaagaTATCGTCAAAATCTTCGGAGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGa                         >  2:674640/1‑149 (MQ=255)
                                                                                                                              agaTATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGAcc                       >  2:3652938/1‑149 (MQ=255)
                                                                                                                                  atCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCtt                   <  1:1201704/149‑1 (MQ=255)
                                                                                                                                  atCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCtt                   >  1:3663895/1‑149 (MQ=255)
                                                                                                                                  atCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCtt                   <  1:3548951/149‑1 (MQ=255)
                                                                                                                                  atCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCtt                   >  1:4881526/1‑149 (MQ=255)
                                                                                                                                   tCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCttt                  >  1:1961445/1‑149 (MQ=255)
                                                                                                                                      tCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCttgaag               <  2:3626222/149‑5 (MQ=255)
                                                                                                                                      tCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGaa               <  1:4980885/149‑1 (MQ=255)
                                                                                                                                          aaTCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCgg           <  1:3114859/149‑1 (MQ=255)
                                                                                                                                          aaTCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCgg           >  2:386099/1‑149 (MQ=255)
                                                                                                                                            tCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGc         >  2:4083395/1‑149 (MQ=255)
                                                                                                                                            tCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGc         >  1:4702307/1‑149 (MQ=255)
                                                                                                                                            tCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGc         >  1:4157521/1‑149 (MQ=255)
                                                                                                                                              ttCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCtt       <  1:4713963/149‑1 (MQ=255)
                                                                                                                                              ttCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCtt       <  1:1323060/149‑1 (MQ=255)
                                                                                                                                               tCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTc      <  1:1112075/149‑1 (MQ=255)
                                                                                                                                                  gCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTCAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGc   >  1:5000940/1‑149 (MQ=255)
                                                                                                                                                  gCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGc   >  1:19176/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATGGTACGCTTATGCTCAGTAAGCGTTTCTTCGAGGAATACCGCCACTTCGATTTCAATATCGTGGCCGCGTGTGGCGGGGCGTTGACGGCTCTGGCGGGCATGCTGACCGAAAATTGACGAGAAGATATCGTCAAAATCTTCGGCGTTAAAACTCTGACCGTCGCCATGGTGGAACTGACGGTTAAATTGCGGATCGTTGCGATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCA  >  CP009273/1058814‑1059109

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: