Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,058,131 G→T 100% H61N (CAT→AAT)  cbpM ← modulator of CbpA co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,058,1310GT97.7% 533.0 / 3.8 172H61N (CAT→AAT) cbpMmodulator of CbpA co‑chaperone
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (74/94);  total (76/96)
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.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 172 total aligned reads displayed.

CGCATCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATGCAGGCAAAA  >  CP009273/1057987‑1058275
                                                                                                                                                |                                                                                                                                                
cGCATCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACg                                                                                                                                               <  1:3018038/148‑1 (MQ=255)
cGCATCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACg                                                                                                                                               >  1:270695/1‑148 (MQ=255)
cGCATCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACg                                                                                                                                               <  1:2529576/148‑1 (MQ=255)
   aTCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAg                                                                                                                                            <  1:2350050/148‑1 (MQ=255)
    tCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAgg                                                                                                                                           >  2:1839584/1‑148 (MQ=255)
     ccGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGc                                                                                                                                          >  2:2221603/1‑148 (MQ=255)
     ccGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGc                                                                                                                                          <  1:2697865/148‑1 (MQ=255)
      cGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGTCCAGTCCAGAGCCAGTTCATTACGCAGGCg                                                                                                                                         <  2:974203/148‑1 (MQ=255)
      cGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCg                                                                                                                                         <  1:1619473/148‑1 (MQ=255)
       ggCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCGTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCCGTTCAGTCCGCAGGCGt                                                                                                                                        >  2:2693390/1‑148 (MQ=255)
           aaCACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACcg                                                                                                                                    <  1:603710/148‑1 (MQ=255)
           aaCACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACcg                                                                                                                                    <  1:603762/148‑1 (MQ=255)
            acacacATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACcgc                                                                                                                                   >  2:2549076/1‑148 (MQ=255)
                 caTCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGtt                                                                                                                              <  1:1267632/148‑1 (MQ=255)
                  aTCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTg                                                                                                                             <  1:3349835/148‑1 (MQ=255)
                        gATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGcaccac                                                                                                                       >  2:3318162/1‑148 (MQ=255)
                         aTGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACcaca                                                                                                                      <  1:1865233/148‑1 (MQ=255)
                         aTGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACcaca                                                                                                                      <  2:1063377/148‑1 (MQ=255)
                          tGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACaa                                                                                                                     <  1:324502/148‑1 (MQ=255)
                          tGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACaa                                                                                                                     <  1:2803588/148‑1 (MQ=255)
                           gAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAAt                                                                                                                    <  1:1050749/148‑1 (MQ=255)
                             gCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATgg                                                                                                                  <  1:413573/148‑1 (MQ=255)
                             gCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATgg                                                                                                                  <  1:3450871/148‑1 (MQ=255)
                              cTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATggc                                                                                                                 >  2:1291240/1‑148 (MQ=255)
                                aCAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATggcgg                                                                                                               <  1:2638679/148‑1 (MQ=255)
                                 cAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATggcggc                                                                                                              >  1:1586821/1‑148 (MQ=255)
                                  aaaCCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCa                                                                                                             <  2:474155/148‑1 (MQ=255)
                                    aCCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATg                                                                                                           >  1:554284/1‑148 (MQ=255)
                                    aCCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGGGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATg                                                                                                           >  2:477586/1‑148 (MQ=255)
                                     ccGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATgg                                                                                                          <  1:989308/148‑1 (MQ=255)
                                     ccGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATgg                                                                                                          <  2:1939001/148‑1 (MQ=255)
                                     ccGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATgg                                                                                                          >  2:2856850/1‑148 (MQ=255)
                                         gAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGgtcgt                                                                                                      <  1:2818776/148‑1 (MQ=255)
                                         gAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGgtcgt                                                                                                      >  1:3179039/1‑148 (MQ=255)
                                         gAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGgtcgt                                                                                                      <  1:1398049/148‑1 (MQ=255)
                                            aGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCaa                                                                                                   >  2:647776/1‑148 (MQ=255)
                                            aGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCaa                                                                                                   >  2:2379995/1‑148 (MQ=255)
                                                  tGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATAccc                                                                                             <  1:198792/148‑1 (MQ=255)
                                                  tGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGCCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATAccc                                                                                             <  1:1527516/148‑1 (MQ=255)
                                                      gcagcaGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgt                                                                                         <  2:2961917/148‑1 (MQ=255)
                                                       cagcagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgtt                                                                                        <  2:114396/148‑1 (MQ=255)
                                                       cagcagGCGGTTGTCCTGCTTCAGGTGCGCAATATCATCCATGAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgtt                                                                                        <  2:171677/148‑1 (MQ=255)
                                                         gcagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgttgt                                                                                      <  2:1973109/148‑1 (MQ=255)
                                                         gcagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgttgt                                                                                      >  1:3297671/1‑148 (MQ=255)
                                                         gcagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgttgt                                                                                      <  1:2642474/148‑1 (MQ=255)
                                                          cagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgttgtt                                                                                     >  1:3358139/1‑148 (MQ=255)
                                                          cagGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGgttgtt                                                                                     >  2:2296936/1‑148 (MQ=255)
                                                            ggCGGTTTTCCTGCTTCAGGTGCGCAATATCATGCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTc                                                                                   <  2:40515/148‑1 (MQ=255)
                                                            ggCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTc                                                                                   <  2:40446/148‑1 (MQ=255)
                                                            ggCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTc                                                                                   <  2:1586821/148‑1 (MQ=255)
                                                            ggCGGTTTTCCTGCTACAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCATAGCCAGATCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGACGACAAATACCCACGTTGTTTc                                                                                   >  2:985184/1‑148 (MQ=255)
                                                             gCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGTGGCATGGTCGTCAAATACCCAGGTTGTTTcc                                                                                  <  2:3058357/148‑1 (MQ=255)
                                                             gCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTcc                                                                                  >  1:1489033/1‑148 (MQ=255)
                                                             gCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTcc                                                                                  <  2:1548231/148‑1 (MQ=255)
                                                              cGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCt                                                                                 <  2:70276/148‑1 (MQ=255)
                                                                 ttttCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGa                                                                              >  1:4381/1‑148 (MQ=255)
                                                                 ttttCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGa                                                                              >  2:1356704/1‑148 (MQ=255)
                                                                  tttCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGAt                                                                             >  2:542681/1‑148 (MQ=255)
                                                                    tCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGAtct                                                                           <  1:1696122/148‑1 (MQ=255)
                                                                    tCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGAtct                                                                           <  1:1696143/148‑1 (MQ=255)
                                                                     ccTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGAtctc                                                                          <  1:2796296/148‑1 (MQ=255)
                                                                       tGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAc                                                                        >  2:449761/1‑148 (MQ=255)
                                                                        gCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                                       <  2:3426974/148‑1 (MQ=255)
                                                                        gCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                                       >  1:1546891/1‑148 (MQ=255)
                                                                        gCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                                       <  2:1689880/148‑1 (MQ=255)
                                                                        gCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                                       <  2:1808503/148‑1 (MQ=255)
                                                                         cTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCAcgc                                                                      >  1:1482869/1‑148 (MQ=255)
                                                                              ggTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTc                                                                 >  1:2484763/1‑148 (MQ=255)
                                                                                    cAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAaccac                                                           >  1:2480584/1‑148 (MQ=255)
                                                                                     aaTATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAaccacc                                                          <  1:467710/148‑1 (MQ=255)
                                                                                       tatCATCCATTCACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCAcccc                                                        <  2:1577403/148‑1 (MQ=255)
                                                                                       tatCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCAcccc                                                        <  2:3312225/148‑1 (MQ=255)
                                                                                        atcatcCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCAccccc                                                       <  1:744233/148‑1 (MQ=255)
                                                                                          catcCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCaa                                                     >  1:1810751/1‑148 (MQ=255)
                                                                                            tcCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAAc                                                   >  2:3124684/1‑148 (MQ=255)
                                                                                             ccATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAAcc                                                  >  1:2623960/1‑148 (MQ=255)
                                                                                                 tAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCATTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGGGGCATGGTCGTCAAATACCCTGGTTGTTTCCTGGATCTCACGGGGTTCAACCACCCCCAAACCGACa                                              <  1:983479/148‑1 (MQ=255)
                                                                                                 tAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTAAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACa                                              <  1:2740740/148‑1 (MQ=255)
                                                                                                   aCGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAAt                                            >  1:725768/1‑148 (MQ=255)
                                                                                                    cGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAAtt                                           <  1:1036148/148‑1 (MQ=255)
                                                                                                       cAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATttca                                        >  1:177143/1‑148 (MQ=255)
                                                                                                        aGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATttcat                                       <  2:2932017/148‑1 (MQ=255)
                                                                                                          cgcCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATttcattc                                     <  1:1486954/148‑1 (MQ=255)
                                                                                                          cgcCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATttcattc                                     <  2:2137945/148‑1 (MQ=255)
                                                                                                             cACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAAc                                  >  1:2600417/1‑148 (MQ=255)
                                                                                                              aCCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAAct                                 >  2:147235/1‑148 (MQ=255)
                                                                                                              aCCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAAct                                 <  2:504303/148‑1 (MQ=255)
                                                                                                              aCCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCAGTCAAct                                 >  2:146974/1‑148 (MQ=255)
                                                                                                               ccGCGATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAActc                                >  2:1589458/1‑148 (MQ=255)
                                                                                                                  cgATCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAActcctc                             <  2:1790759/148‑1 (MQ=255)
                                                                                                                    aTCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCtt                           <  2:1394888/148‑1 (MQ=255)
                                                                                                                     tCCCCGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTc                          >  2:1768235/1‑148 (MQ=255)
                                                                                                                      ccccGGCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCa                         >  1:1364574/1‑148 (MQ=255)
                                                                                                                          ggCCAGTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAgaga                     <  1:2928479/148‑1 (MQ=255)
                                                                                                                               gTCCAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCg                <  1:1291240/148‑1 (MQ=255)
                                                                                                                                 ccAGAGCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGt              <  1:1008936/148‑1 (MQ=255)
                                                                                                                                    gatccAGTTCATTACTCAATCGAACCGCGCGTTGCACCACAATGGCTGCATGGACGTAAAATACAAACGTTGTTTCCTTGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATg           >  2:3278982/4‑148 (MQ=255)
                                                                                                                                      gCCAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATGCa         >  2:791301/1‑148 (MQ=255)
                                                                                                                                       ccAGTTCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATGCAg        >  2:390824/1‑148 (MQ=255)
                                                                                                                                            tCATTACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATGCAGGCaaa   <  2:1810751/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGCATCCGGCAAACACACATCACGGATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGGCGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGTCAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAGAGATGCCGGTATGCAGGCAAAA  >  CP009273/1057987‑1058275

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: