Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,672,918 A→C 5.1% 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 5.1% 579.9 / 21.0 196S79R (AGC→CGC) ydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (72/114);  new base C (5/5);  total (77/119)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 197 total aligned reads displayed.

CAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAG  >  CP009273/1672771‑1673064
                                                                                                                                                   |                                                                                                                                                  
 aaGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAg                                                                                                                                                   <  1:2495790/148‑1 (MQ=255)
 aaGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAg                                                                                                                                                   <  2:3279401/148‑1 (MQ=255)
     ggCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGt                                                                                                                                               <  2:2534283/148‑1 (MQ=255)
     ggCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGt                                                                                                                                               <  2:2536933/148‑1 (MQ=255)
      gCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTa                                                                                                                                              <  2:3236795/148‑1 (MQ=255)
       cAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATTCTATTGGCGAAGCGGTGAAATCCGGTTCTCGTCGCGCCGATAAAGAAGGTGCCCCCTCTTTTTATGTTGTCGACACAACTGATTTTGGTAACACCGGTaa                                                                                                                                             >  2:2853955/1‑148 (MQ=255)
         gCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACt                                                                                                                                           <  2:3177612/148‑1 (MQ=255)
          cGGTTAATCCTTTTGACCGTGTAGTGGTTACCGGTAGTTTTAATTCTATTGGCGAATCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCTACACTTCTGATTTTGGTAACAGCGGTAACTg                                                                                                                                          <  1:1216386/148‑1 (MQ=255)
          cGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTg                                                                                                                                          <  1:193910/148‑1 (MQ=255)
              tAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCgt                                                                                                                                      <  2:2111439/148‑1 (MQ=255)
                 aCCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGtgtg                                                                                                                                   <  2:1920650/148‑1 (MQ=255)
                  ccTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTgg                                                                                                                                  <  1:1627446/148‑1 (MQ=255)
                        gACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTGGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTg                                                                                                                            <  1:2372948/148‑1 (MQ=255)
                          ccGTGTAGTGGTTACCGGTCGTTTTAATTCTATTTGCGAACCGGTGTAAGCCGTTTCTCGGCGCGCAGATAAATAAGGTGCCGCCACTTTTTATTTTGTGGACACTTCGGATTTTTGTAAAAGCGGTAACAGGCGGGTGGTCGCTGAc                                                                                                                          >  2:565659/1‑148 (MQ=255)
                             tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACctc                                                                                                                       <  2:3189201/148‑1 (MQ=255)
                              gtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACctct                                                                                                                      <  1:2295336/148‑1 (MQ=255)
                               ttgTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCta                                                                                                                     <  2:2416559/146‑1 (MQ=255)
                                aGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCtat                                                                                                                    <  1:2289954/148‑1 (MQ=255)
                                  tGGTTACCGGTCGTTTTAATTCTATTGGCGAAGCGGTGAAACCCCGTTATCGGCGCGCACATAAAAAACGTGCCCCCTCTTTTTATGTTGTAGACACATATGATTTTGGTAACACCCGTAAAATGCGGGTGGTCGCTGACCACTATaa                                                                                                                  >  2:2493570/1‑148 (MQ=255)
                                     ttACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGc                                                                                                               <  1:791345/148‑1 (MQ=255)
                                       aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCg                                                                                                             <  2:3063338/148‑1 (MQ=255)
                                        ccGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                                                                                            <  1:2613346/148‑1 (MQ=255)
                                           gTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTCACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGc                                                                                                         <  1:2596896/148‑1 (MQ=255)
                                            tCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCt                                                                                                        <  2:200233/148‑1 (MQ=255)
                                               ttttAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGACGATGCTGaa                                                                                                     >  1:946954/1‑148 (MQ=255)
                                                   aaTGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAg                                                                                                 >  1:1442201/1‑148 (MQ=255)
                                                     tGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCGCTATAAAGCCGATGCTGAAAAAGCa                                                                                               >  1:3330625/1‑148 (MQ=255)
                                                     tGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCa                                                                                               <  2:1812769/148‑1 (MQ=255)
                                                         aTTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCagaag                                                                                           <  1:1671907/148‑1 (MQ=255)
                                                          ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCagaaga                                                                                          <  1:805513/148‑1 (MQ=255)
                                                          ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCagaaga                                                                                          <  2:859320/148‑1 (MQ=255)
                                                            ggCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGTCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGaaa                                                                                        <  1:3446581/148‑1 (MQ=255)
                                                            ggCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGaaa                                                                                        <  2:331620/148‑1 (MQ=255)
                                                             gCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAAc                                                                                       >  2:1102348/1‑148 (MQ=255)
                                                              cGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACa                                                                                      <  2:1388970/148‑1 (MQ=255)
                                                                  gCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTa                                                                                  >  1:1041569/1‑148 (MQ=255)
                                                                  gCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTa                                                                                  <  1:3006754/148‑1 (MQ=255)
                                                                   cGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTaa                                                                                 <  1:1151286/148‑1 (MQ=255)
                                                                    ggTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAAt                                                                                <  1:3548225/148‑1 (MQ=255)
                                                                      tGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcg                                                                              <  1:2708706/148‑1 (MQ=255)
                                                                       gAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcgc                                                                             >  2:2670617/1‑148 (MQ=255)
                                                                        aaaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcgcg                                                                            >  2:2331357/1‑148 (MQ=255)
                                                                        aaaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcgcg                                                                            >  1:3018805/1‑148 (MQ=255)
                                                                         aaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGt                                                                           <  1:3439089/148‑1 (MQ=255)
                                                                         aaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATAGCGt                                                                           >  1:681528/1‑148 (MQ=255)
                                                                         aaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGAGGCTGAAAAAGCAGAAGAAACAAGTAATCGCGt                                                                           >  1:3147438/1‑148 (MQ=255)
                                                                         aaGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGACGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGt                                                                           >  2:1447429/1‑148 (MQ=255)
                                                                           gCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTaa                                                                         <  2:1593933/148‑1 (MQ=255)
                                                                           gCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTaa                                                                         >  2:405929/1‑148 (MQ=255)
                                                                            ccGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAAt                                                                        >  2:2893952/1‑148 (MQ=255)
                                                                            ccGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAAt                                                                        <  1:1338620/148‑1 (MQ=255)
                                                                              gtgtctcGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTa                                                                      <  2:1951196/145‑1 (MQ=255)
                                                                              gTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTa                                                                      <  2:156723/148‑1 (MQ=255)
                                                                                ttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAAc                                                                    <  1:1505332/148‑1 (MQ=255)
                                                                                ttCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAAc                                                                    >  1:713552/1‑148 (MQ=255)
                                                                                 tctcGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACg                                                                   <  1:1583495/148‑1 (MQ=255)
                                                                                   tcgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgt                                                                 >  2:1023511/1‑148 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                <  1:1580341/148‑1 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                >  1:1904392/1‑148 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                <  1:3259065/148‑1 (MQ=255)
                                                                                    cgtcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGgtg                                                                <  1:1114317/148‑1 (MQ=255)
                                                                                      tcgCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGtgtt                                                              >  2:1303305/1‑148 (MQ=255)
                                                                                       cgcgCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGtgttg                                                             <  1:2062602/148‑1 (MQ=255)
                                                                                           cAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGa                                                         <  2:2160744/148‑1 (MQ=255)
                                                                                               tAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTg                                                     <  2:899750/148‑1 (MQ=255)
                                                                                               tAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTg                                                     >  1:588743/1‑148 (MQ=255)
                                                                                                  agaagGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGGAATTAACGGTGTTGTCGAACTGCCg                                                  >  1:3051954/1‑148 (MQ=255)
                                                                                                   gaagGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGa                                                 >  2:1311318/1‑148 (MQ=255)
                                                                                                    aagGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGaa                                                >  1:830728/1‑148 (MQ=255)
                                                                                                     agGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGaaa                                               >  1:1602279/1‑148 (MQ=255)
                                                                                                      ggTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAg                                              >  2:94579/1‑148 (MQ=255)
                                                                                                      ggTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAg                                              >  2:1892474/1‑148 (MQ=255)
                                                                                                             ccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGc                                       <  1:3419740/148‑1 (MQ=255)
                                                                                                             ccTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGc                                       <  1:2103414/148‑1 (MQ=255)
                                                                                                              ctctTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCt                                      >  2:1086216/1‑148 (MQ=255)
                                                                                                                ctTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGt                                    <  1:1023511/148‑1 (MQ=255)
                                                                                                                 tttttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAATCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGtt                                   >  1:1295166/1‑148 (MQ=255)
                                                                                                                 tttttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGtt                                   <  2:1934138/148‑1 (MQ=255)
                                                                                                                 tttttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGtt                                   <  2:1779931/148‑1 (MQ=255)
                                                                                                                  ttttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTc                                  <  1:2493574/148‑1 (MQ=255)
                                                                                                                    ttATGTTGTCGACACTTCTGATTTTGGTAACCGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTg                                <  1:898399/148‑1 (MQ=255)
                                                                                                                    ttATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTg                                >  1:2180706/1‑148 (MQ=255)
                                                                                                                      aTGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGAt                              >  2:1654907/1‑148 (MQ=255)
                                                                                                                        gttgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTg                            >  1:1840038/1‑148 (MQ=255)
                                                                                                                        gttgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTg                            >  1:791910/1‑148 (MQ=255)
                                                                                                                        gttgtCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTg                            >  1:1349224/1‑148 (MQ=255)
                                                                                                                              gACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGt                      <  2:2772082/148‑1 (MQ=255)
                                                                                                                                cacTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGttt                    >  1:397387/1‑148 (MQ=255)
                                                                                                                                 acTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTg                   >  1:1762428/1‑148 (MQ=255)
                                                                                                                                    tCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATAAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCCTTTGaca                >  2:812686/1‑148 (MQ=255)
                                                                                                                                     cTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGacac               >  1:871258/1‑148 (MQ=255)
                                                                                                                                          tttGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCa          >  1:3102792/1‑148 (MQ=255)
                                                                                                                                             ggTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCg       <  1:2853955/148‑1 (MQ=255)
                                                                                                                                               tAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTc     <  2:2622601/148‑1 (MQ=255)
                                                                                                                                                aaCAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTcc    <  1:1557439/148‑1 (MQ=255)
                                                                                                                                                aaCAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTcc    >  2:181364/1‑148 (MQ=255)
                                                                                                                                                aaCAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTcc    >  2:2085596/1‑148 (MQ=255)
                                                                                                                                                aaCAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTcc    <  1:3321046/148‑1 (MQ=255)
                                                                                                                                                  cAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAg  <  2:949620/148‑1 (MQ=255)
                                                                                                                                                  cAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAg  >  2:2053106/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAACGGTGTTGTCGAACTGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAG  >  CP009273/1672771‑1673064

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: