Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 165,759 G→T 44.8% E597* (GAG→TAG)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273165,7590GT44.8% 139.2 / 203.5 191E597* (GAG→TAG) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (49/56);  new base T (39/47);  total (88/103)
Fisher's exact test for biased strand distribution p-value = 8.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 195 total aligned reads displayed.

CCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACAC  >  CP009273/165612‑165905
                                                                                                                                                   |                                                                                                                                                  
ccGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCg                                                                                                                                                    >  1:1038233/1‑148 (MQ=255)
 cGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCga                                                                                                                                                   >  2:3338592/1‑148 (MQ=255)
  gTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCtag                                                                                                                                                  <  1:1605462/148‑1 (MQ=255)
     tAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATc                                                                                                                                               <  2:1814979/148‑1 (MQ=255)
     tAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATc                                                                                                                                               <  1:1378385/148‑1 (MQ=255)
         ggTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGCCCCTGAGGGTTCCTTCTTCTCGGTTGACGGTGGCTAGATCcgcg                                                                                                                                           >  2:2480135/1‑148 (MQ=255)
           tAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCa                                                                                                                                         <  1:468751/148‑1 (MQ=255)
            aaGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCAc                                                                                                                                        >  2:1487885/1‑148 (MQ=255)
            aaGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCAc                                                                                                                                        >  1:1463113/1‑148 (MQ=255)
            aaGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATAGTCCGATTGTATTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATTGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCAGAGCAc                                                                                                                                        >  1:3236777/1‑148 (MQ=255)
               cAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTg                                                                                                                                     >  1:2793877/1‑148 (MQ=255)
               cAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAAACTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTg                                                                                                                                     >  2:1682809/1‑148 (MQ=255)
                  tATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGGAGGTGGCTAGATCCGCGCACGTGGCg                                                                                                                                  <  1:221155/148‑1 (MQ=255)
                  tATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCg                                                                                                                                  <  1:3165276/148‑1 (MQ=255)
                   aTGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGt                                                                                                                                 >  2:3117661/1‑148 (MQ=255)
                   aTGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGt                                                                                                                                 <  1:3555503/148‑1 (MQ=255)
                     gAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAg                                                                                                                               >  1:3564980/1‑148 (MQ=255)
                       aGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGaa                                                                                                                             <  1:1794555/148‑1 (MQ=255)
                          cGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATc                                                                                                                          <  1:2340091/148‑1 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGaa                                                                                                                       >  2:3227727/1‑148 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGaa                                                                                                                       >  1:1786372/1‑148 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGaa                                                                                                                       <  1:112655/148‑1 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGaa                                                                                                                       <  1:112613/148‑1 (MQ=255)
                                    tatGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAg                                                                                                                <  1:1938044/148‑1 (MQ=255)
                                     atGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgc                                                                                                               >  1:3071782/1‑148 (MQ=255)
                                     atGTACCGGAAGATCGTCCGAGTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgc                                                                                                               <  2:99062/148‑1 (MQ=255)
                                        tACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggc                                                                                                            >  1:2248924/1‑148 (MQ=255)
                                         aCCGGAAGATCGTCCGATTGTAGTTACTGGTGCCTTGTATAATCTCACTAAAACCAACACCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTATAAATCGAAGCGAAAgcggcg                                                                                                           <  2:2307456/148‑1 (MQ=255)
                                         aCCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggcg                                                                                                           <  2:2307204/148‑1 (MQ=255)
                                         aCCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggcg                                                                                                           >  1:1743657/1‑148 (MQ=255)
                                          ccGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAAAAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGATAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGgcgc                                                                                                          >  1:945798/1‑148 (MQ=255)
                                              aaGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTc                                                                                                      >  2:620643/1‑148 (MQ=255)
                                               aGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCg                                                                                                     <  2:3170295/148‑1 (MQ=255)
                                                  tCGTCCTATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCTAAAGCGGCGCTGTCGGCg                                                                                                  <  2:931883/148‑1 (MQ=255)
                                                      ccGTTTGTTGTTCGTGGTGCCGTTTATAATCTCACTAAAACAAACACCCTGAGGGGGGACCGTGAGGGTTCCTTCTTCTGGGTTGAAGGTGGCGAGAGCCGCGCCCGTGGCGTATAAATCGACGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                              <  1:108853/148‑1 (MQ=255)
                                                      ccGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                              >  2:3092418/1‑148 (MQ=255)
                                                         tttGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACCACCTGAGGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTa                                                                                           <  2:1054568/147‑1 (MQ=255)
                                                         aTTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTa                                                                                           >  2:2284291/1‑148 (MQ=255)
                                                             tAGTTACTGGTGCCGTGTATAATCTCAATAAACCCAACAACCTGATGGCGGACCCGGCGGGTTCCTTCTTCTCGGTTGATGGTGGCTAGATCCGCGCCCGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGt                                                                                       <  1:1907591/148‑1 (MQ=255)
                                                               gTTACTGGTGCCGTGTATAATCTCACTAAAACCAACACCCTGTTGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCAGGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAg                                                                                     <  1:1203034/148‑1 (MQ=255)
                                                               gTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAg                                                                                     >  1:345765/1‑148 (MQ=255)
                                                                ttACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGt                                                                                    >  1:698053/1‑148 (MQ=255)
                                                                   gtGGGGCCGTGTATAATCTCACAAAAACCAACACCCTGTGGGCGGACCCTGCGGGTTCGTTCTTCGCGGTTGACGGGGGCTACACCCGCGCACGTGGCTTACAAATCGAAGCGAACGCGGCGCTGTCGGCGAGTGTTAACGTAGTCgg                                                                                 <  1:453061/147‑1 (MQ=255)
                                                                     ggTGCCTTGTAACATCTCACAAAAACAAACAACCTGATGGCGGACCCTGGGGGTTCATTCTTGTGGGTTGAAGTTGGATACATCCGCGCACGTGGCGTATAAATCCAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGtt                                                                               <  1:2206934/148‑1 (MQ=255)
                                                                     ggTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGtt                                                                               <  2:237946/148‑1 (MQ=255)
                                                                      gTGCCTTGTATAATCTCACTAACACCAACAACCTGATGGGGGACCCTGCGGGTTCCTTCTTCTCGGTTGATGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTc                                                                              <  2:1291735/148‑1 (MQ=255)
                                                                      gTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTc                                                                              >  1:3187270/1‑148 (MQ=255)
                                                                      gTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTATCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAAACGAAACGAAAGCGGCGCTGTCGGCGAGTGTTAACGAAGTCGGGTc                                                                              >  1:2162438/1‑148 (MQ=255)
                                                                       tGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCt                                                                             >  2:1593731/1‑148 (MQ=255)
                                                                        gCCGTGTATAATCTCACTAAAACCAACACCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCtt                                                                            <  2:1550065/148‑1 (MQ=255)
                                                                              tataATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACtt                                                                      <  2:3007232/148‑1 (MQ=255)
                                                                               ataATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTa                                                                     >  2:187505/1‑148 (MQ=255)
                                                                               ataATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTa                                                                     >  1:2204266/1‑148 (MQ=255)
                                                                                taATCTCACTAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTGGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTac                                                                    <  1:2227710/148‑1 (MQ=255)
                                                                                taATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTac                                                                    >  2:1633483/1‑148 (MQ=255)
                                                                                   tctcACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACAcc                                                                 >  2:780978/1‑148 (MQ=255)
                                                                                   tctcACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACAcc                                                                 >  1:901435/1‑148 (MQ=255)
                                                                                      cACTAAAACCAACAACCTTATGGCGGACCATGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTTGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGAt                                                              <  1:2056288/148‑1 (MQ=255)
                                                                                      cACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGAt                                                              >  2:3020688/1‑148 (MQ=255)
                                                                                         tAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCg                                                           >  2:1928594/1‑148 (MQ=255)
                                                                                          aaaaCCAACACCCTGTGGGCGGACCCTGCGGGTTCATTCTTCTGGGTTGATGGTGGCGAGATCCGCGCCGGTGGCGTAGAAATCGAAGCTAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCgg                                                          <  1:1547079/148‑1 (MQ=255)
                                                                                          aaaaCCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCGATATCCGCGCAGGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCgg                                                          <  1:3544071/148‑1 (MQ=255)
                                                                                             aCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                       >  2:774178/1‑148 (MQ=255)
                                                                                             aCCAACAACCTGAGGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGATGGTGGCGAGATCCGCGCACGTGGCGTATAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                       <  1:2490748/148‑1 (MQ=255)
                                                                                              ccAACAACCTGATGGCGGACCCTGTGGGTTCCTTGTTCTGGTTTGAAGGTGGCGAGATCCGCGCACGTGGCTTAGAAATCTATGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATa                                                      <  2:2162438/148‑1 (MQ=255)
                                                                                               caacaaCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATac                                                     >  2:448268/1‑148 (MQ=255)
                                                                                               caacaaCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATac                                                     >  1:2663909/1‑148 (MQ=255)
                                                                                                aacaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATaca                                                    <  2:2980264/148‑1 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 <  2:2601446/148‑1 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 <  2:814683/148‑1 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 <  2:1884725/148‑1 (MQ=255)
                                                                                                   aacCTGAGGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 <  2:1884686/148‑1 (MQ=255)
                                                                                                    acCTGCTGGCGGACCCTGCGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcaccac                                                <  2:2445567/148‑1 (MQ=255)
                                                                                                    acCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcaccac                                                >  1:3008923/1‑148 (MQ=255)
                                                                                                    acCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcaccac                                                <  2:543857/148‑1 (MQ=255)
                                                                                                     cctgtgggCGGACCCTGCGGTTTCCTTCTTCTCGGTTGATGGGGGCTAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                               <  1:416746/142‑1 (MQ=255)
                                                                                                     ccTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                               >  2:2901084/1‑148 (MQ=255)
                                                                                                        gagggCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAt                                            <  2:725549/145‑1 (MQ=255)
                                                                                                        gATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAt                                            >  1:866679/1‑148 (MQ=255)
                                                                                                        gATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAt                                            >  1:2608310/1‑148 (MQ=255)
                                                                                                         aTGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAta                                           >  1:1986221/1‑148 (MQ=255)
                                                                                                           ggCGGACCCTGTGGGTTGCTTCTTCTCGTTTGATGGTGGCTTGATCCGCGCACTTGGGTTATAAATCGAAGCTAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtact                                         <  1:3359655/148‑1 (MQ=255)
                                                                                                             cGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtactac                                       <  2:22508/148‑1 (MQ=255)
                                                                                                             cGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtactac                                       <  1:151608/148‑1 (MQ=255)
                                                                                                                aCCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCta                                    >  2:2659531/1‑148 (MQ=255)
                                                                                                                    tGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATaaa                                <  2:3238644/148‑1 (MQ=255)
                                                                                                                     gAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAg                               >  2:599606/1‑148 (MQ=255)
                                                                                                                        ggTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCa                            >  2:501693/1‑148 (MQ=255)
                                                                                                                          ttCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTATAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAAt                          <  1:3282212/148‑1 (MQ=255)
                                                                                                                          ttCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAAt                          <  2:3452376/148‑1 (MQ=255)
                                                                                                                             cttcttCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACg                       <  2:3536861/148‑1 (MQ=255)
                                                                                                                             cttcttCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACg                       >  1:1115134/1‑148 (MQ=255)
                                                                                                                              ttcttcTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGc                      <  2:1981221/148‑1 (MQ=255)
                                                                                                                                       gTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAgg             >  2:2124660/1‑148 (MQ=255)
                                                                                                                                       gTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAgg             <  2:1038233/148‑1 (MQ=255)
                                                                                                                                         tGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTg           >  1:3193037/1‑148 (MQ=255)
                                                                                                                                             ggTGGCGCGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaa       <  1:3349022/148‑1 (MQ=255)
                                                                                                                                             ggTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaa       >  2:3106285/1‑148 (MQ=255)
                                                                                                                                             ggTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaa       >  1:2902969/1‑148 (MQ=255)
                                                                                                                                              gTGGCTAGATCCGCGCACGTGGCGTATAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaaa      <  1:80716/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACAC  >  CP009273/165612‑165905

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