Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,672,918 A→C 8.8% S79R (AGC→CGC)  ydgH → DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,672,9180AC 8.8% 498.7 / 39.1 189S79R (AGC→CGC) ydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (66/106);  new base C (5/12);  total (71/118)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

CAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGG  >  CP009273/1672771‑1673065
                                                                                                                                                   |                                                                                                                                                   
cAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACa                                                                                                                                                     <  2:3423859/148‑1 (MQ=255)
  aGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGc                                                                                                                                                   <  1:1716319/148‑1 (MQ=255)
   gCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCg                                                                                                                                                  <  2:3439567/148‑1 (MQ=255)
     ggCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAGGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGt                                                                                                                                                <  1:3118336/148‑1 (MQ=255)
      gCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTa                                                                                                                                               <  2:2518698/148‑1 (MQ=255)
        aGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAAc                                                                                                                                             <  2:3203080/148‑1 (MQ=255)
         gCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACt                                                                                                                                            <  2:2106870/148‑1 (MQ=255)
         gCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACt                                                                                                                                            <  2:1637827/148‑1 (MQ=255)
             ttAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCg                                                                                                                                        <  1:2764134/148‑1 (MQ=255)
                aaCCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATTCTATTTGCGAAGCGGTGTAAGCCGTTTCTAGTCGCGCACATAAAGAATGTGCCGCCGCTTTTTATTTTTTCGACACTTCTGATTTTGTTAAAAGCGGTAAATGGCgtgt                                                                                                                                     >  2:2479483/1‑148 (MQ=255)
                 aCCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGtgtg                                                                                                                                    <  1:359422/148‑1 (MQ=255)
                  ccTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATTCTATTGGCGAAGCGGTGAAATCCGGTTGTCGTCGCGCCGATAAAGAAGGTGCCGCCTCTTTTTTTTGTGTAGACACATCTGATTTTGGTAACAGCGGTAACAGGCGGGTgg                                                                                                                                   >  2:1236732/1‑148 (MQ=255)
                  ccTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTgg                                                                                                                                   <  2:980558/148‑1 (MQ=255)
                     tttGACCGTGTAGTGGTTACCGGTCGTTTTTAATAAATTTTGGAATCCGTGAAACACGTTTTTAGTCGCGCAGATAAAAAAGGTGCCGCCTCTTTTTATTTTTTTGAAACATATTTTTTTTTTAACAACGGGAAATCGCGTGTggggg                                                                                                                                >  2:2897769/1‑145 (MQ=255)
                     tttGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCg                                                                                                                                <  2:51149/148‑1 (MQ=255)
                      ttGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGc                                                                                                                               <  1:2583023/148‑1 (MQ=255)
                      ttGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGc                                                                                                                               <  1:2073346/148‑1 (MQ=255)
                         aCCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTTGGGTAGCGGTGAAAGCCCGTTATCGGGGCGCAGATAAAGAATGTGCCGCCACTTTTTATGTTTTCGACACATCTGATTTTTGTAAAAGCGGTTAAATGCGTGTGGTCGCTGa                                                                                                                            >  2:3709156/1‑148 (MQ=255)
                             tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACctc                                                                                                                        <  1:3041060/148‑1 (MQ=255)
                              gtAGTGGTTACCGGTCGTTTTAATTCTATTTGGGAAGCGGGGAAAACCGGTTATCGTCGCGCCCAAAAAGAAGGTGCCGCCACTTTTTATGTTGTAGACACATCTGGTTTTTGTAACACCGGTAACAGGCGTGTGGTCGCTGTCcact                                                                                                                       >  2:2140623/1‑148 (MQ=255)
                                  tGGTTACCGGTCGTTTTAATTCTAATGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAATAAGGTTCCGCCTCTTTTTATTTTGTCGACACTTCTGATTTTTGTAAAAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaa                                                                                                                   >  2:3445399/1‑148 (MQ=255)
                                  tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaa                                                                                                                   <  1:687029/148‑1 (MQ=255)
                                  tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaa                                                                                                                   <  2:3544712/148‑1 (MQ=255)
                                  tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaa                                                                                                                   <  2:1377672/148‑1 (MQ=255)
                                    gTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAg                                                                                                                 <  1:1393734/148‑1 (MQ=255)
                                     ttACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGc                                                                                                                <  1:2585238/148‑1 (MQ=255)
                                       aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCg                                                                                                              <  2:809030/148‑1 (MQ=255)
                                        ccGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                                                                                             <  2:3071009/148‑1 (MQ=255)
                                        ccGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                                                                                             <  2:2514833/148‑1 (MQ=255)
                                        ccGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                                                                                             <  1:336689/148‑1 (MQ=255)
                                          ggTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATg                                                                                                           <  2:1149724/148‑1 (MQ=255)
                                          ggTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATg                                                                                                           <  2:1099801/148‑1 (MQ=255)
                                               ttttAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaa                                                                                                      >  1:262472/1‑148 (MQ=255)
                                                   aaTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAg                                                                                                  <  1:614313/148‑1 (MQ=255)
                                                   aaTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAg                                                                                                  <  1:2938134/148‑1 (MQ=255)
                                                    aTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGc                                                                                                 <  2:74977/148‑1 (MQ=255)
                                                    aTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGc                                                                                                 >  1:1908970/1‑148 (MQ=255)
                                                    aTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGc                                                                                                 <  1:3462776/148‑1 (MQ=255)
                                                      gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCag                                                                                               <  2:210846/148‑1 (MQ=255)
                                                       cTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCaga                                                                                              <  1:498867/148‑1 (MQ=255)
                                                            ggCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGaaa                                                                                         <  1:1900168/148‑1 (MQ=255)
                                                              cGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACa                                                                                       <  1:235925/148‑1 (MQ=255)
                                                              cGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACa                                                                                       <  2:2824307/148‑1 (MQ=255)
                                                                aaGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAATCAGAAGAAACAAg                                                                                     >  1:1309912/1‑148 (MQ=255)
                                                                 aGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGt                                                                                    >  2:3092063/1‑148 (MQ=255)
                                                                   cGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTaa                                                                                  <  1:2510966/148‑1 (MQ=255)
                                                                   cGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTaa                                                                                  <  1:2232760/148‑1 (MQ=255)
                                                                    ggTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAAt                                                                                 <  1:832647/148‑1 (MQ=255)
                                                                        aaaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcgcg                                                                             >  2:556664/1‑148 (MQ=255)
                                                                          aGCCGTTTCTCGTCGCGCAGATAAAGAATGTGCCGCCTCTTTTTATTTTGTCGACACTTATGATTTTTGTAACAGCGGTAAATGGCGTGTGGTCGCTGACCTCTAAAAAACCGAATCTGAAAAATCAGAAGAAACAACTAATAGCGTa                                                                           >  1:3211040/1‑148 (MQ=255)
                                                                          aGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTa                                                                           <  1:2675566/148‑1 (MQ=255)
                                                                          aGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTa                                                                           <  2:2179192/148‑1 (MQ=255)
                                                                          aGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTa                                                                           >  2:2056491/1‑148 (MQ=255)
                                                                          aGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTa                                                                           <  1:1124523/148‑1 (MQ=255)
                                                                           gCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAATCAGAAGAAACAAGTAATCGCGTaa                                                                          <  2:1616244/148‑1 (MQ=255)
                                                                           gCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTaa                                                                          >  1:1166086/1‑148 (MQ=255)
                                                                           gCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTaa                                                                          >  2:3153695/1‑148 (MQ=255)
                                                                               tttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa                                                                      >  1:56577/1‑148 (MQ=255)
                                                                               tttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa                                                                      >  1:3780206/1‑148 (MQ=255)
                                                                               tttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa                                                                      <  2:2303364/148‑1 (MQ=255)
                                                                                ttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAAc                                                                     >  1:2920450/1‑148 (MQ=255)
                                                                                 tctcGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGTCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACg                                                                    <  1:3284218/148‑1 (MQ=255)
                                                                                 tctcGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACg                                                                    <  1:3442057/148‑1 (MQ=255)
                                                                                  ctcGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACgg                                                                   <  2:830499/148‑1 (MQ=255)
                                                                                   tcgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgt                                                                  <  2:2004213/148‑1 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                 <  1:1797638/148‑1 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGAGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                 <  1:2083548/148‑1 (MQ=255)
                                                                                       cgcgCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGtgttg                                                              <  1:1624632/148‑1 (MQ=255)
                                                                                       cgcgCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGtgttg                                                              <  1:2847940/148‑1 (MQ=255)
                                                                                            aGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGGAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACACGTAATCGCGTAATTAACGGTGTTGTCGac                                                         <  1:148148/148‑2 (MQ=255)
                                                                                                aaaGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGc                                                     >  2:1151182/1‑148 (MQ=255)
                                                                                                aaaGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGc                                                     >  1:925056/1‑148 (MQ=255)
                                                                                                   gaagGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCGCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGa                                                  >  1:1408310/1‑148 (MQ=255)
                                                                                                      ggTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAg                                               >  2:221538/1‑148 (MQ=255)
                                                                                                       gTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAATCAGAAGAAACAAGTAATAGCGTAATTAAAGGTGTTGTCGAACTGCCGAAAGa                                              >  1:1695954/1‑148 (MQ=255)
                                                                                                       gTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGa                                              <  1:2799584/148‑1 (MQ=255)
                                                                                                         gccgccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATc                                            <  2:3178070/148‑1 (MQ=255)
                                                                                                         gccgccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATc                                            <  1:189279/148‑1 (MQ=255)
                                                                                                         gccgccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTATATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATc                                            <  1:201197/148‑1 (MQ=255)
                                                                                                            gccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAgg                                         <  2:2349485/148‑1 (MQ=255)
                                                                                                             ccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGc                                        <  2:3584554/148‑1 (MQ=255)
                                                                                                                ctTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGt                                     >  2:1734878/1‑148 (MQ=255)
                                                                                                                          tgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTACTCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGaa                           <  2:378567/148‑1 (MQ=255)
                                                                                                                          tgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGaa                           <  1:793905/148‑1 (MQ=255)
                                                                                                                              gACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGt                       <  1:401400/148‑1 (MQ=255)
                                                                                                                              gACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGt                       <  1:398951/148‑1 (MQ=255)
                                                                                                                              gACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGt                       <  2:3302673/148‑1 (MQ=255)
                                                                                                                               acacTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGtt                      >  2:665188/1‑148 (MQ=255)
                                                                                                                               acacTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGtt                      >  1:217293/1‑148 (MQ=255)
                                                                                                                               acacTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGtt                      >  1:2761134/1‑148 (MQ=255)
                                                                                                                                   ttCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGac                  >  1:2408651/1‑148 (MQ=255)
                                                                                                                                     cTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCATAAGAAACAAGTAATCCCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGacac                >  1:2985075/1‑148 (MQ=255)
                                                                                                                                      tGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACg               >  2:3372056/1‑148 (MQ=255)
                                                                                                                                       gATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACTGTTTGACACgg              <  1:433340/148‑1 (MQ=255)
                                                                                                                                       gATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACTGTTTGACACgg              <  1:3617120/148‑1 (MQ=255)
                                                                                                                                        aTTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGGAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGt             >  1:695279/1‑148 (MQ=255)
                                                                                                                                         ttttGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTc            >  2:2837404/1‑148 (MQ=255)
                                                                                                                                         ttttGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTc            >  2:1565357/1‑148 (MQ=255)
                                                                                                                                          tttGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCa           <  1:1828022/148‑1 (MQ=255)
                                                                                                                                             ggTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCg        <  2:1842709/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGG  >  CP009273/1672771‑1673065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: