Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,505,438 C→A 6.8% D353Y (GAC→TAC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,505,4380CA 6.8% 763.3 / 30.3 219D353Y (GAC→TAC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (102/102);  new base A (6/9);  total (108/111)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 231 total aligned reads displayed.

CAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCA  >  CP009273/4505289‑4505586
                                                                                                                                                     |                                                                                                                                                    
cAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGcgtc                                                                                                                                                      >  1:2391257/1‑150 (MQ=255)
cAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGcgtc                                                                                                                                                      <  2:3886426/150‑1 (MQ=255)
    tAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATGTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTaa                                                                                                                                                  <  1:3166649/150‑1 (MQ=255)
     aGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAAt                                                                                                                                                 >  1:4656464/1‑150 (MQ=255)
     aGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAAt                                                                                                                                                 >  1:3555678/1‑150 (MQ=255)
     aGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAAt                                                                                                                                                 >  2:5192480/1‑150 (MQ=255)
       ccGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATcc                                                                                                                                               >  1:2024973/1‑150 (MQ=255)
        cGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCcg                                                                                                                                              <  1:2268758/150‑1 (MQ=255)
        cGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCcg                                                                                                                                              <  2:4485667/150‑1 (MQ=255)
         gCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGGAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCcgc                                                                                                                                             <  1:4829595/150‑1 (MQ=255)
             cgcAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCga                                                                                                                                         <  1:3372007/150‑1 (MQ=255)
                aGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGAc                                                                                                                                      >  2:3719774/1‑150 (MQ=255)
                 gggTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACa                                                                                                                                     >  2:3809822/1‑150 (MQ=255)
                 gggTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACa                                                                                                                                     >  2:1073818/1‑150 (MQ=255)
                     ttGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGcc                                                                                                                                 >  1:3980464/1‑150 (MQ=255)
                      tGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCa                                                                                                                                >  1:1595303/1‑150 (MQ=255)
                        ggTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCAcc                                                                                                                              >  2:323372/1‑150 (MQ=255)
                         gtgtAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCg                                                                                                                             >  2:2749300/1‑150 (MQ=255)
                            tAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGc                                                                                                                          >  2:781414/1‑150 (MQ=255)
                             aGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCa                                                                                                                         <  2:1791110/150‑1 (MQ=255)
                              gAACCCCTGATTGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCAt                                                                                                                        <  1:4290670/150‑1 (MQ=255)
                              gAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCAt                                                                                                                        <  1:1995199/150‑1 (MQ=255)
                                 ccccTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTc                                                                                                                     <  2:2974142/150‑1 (MQ=255)
                                    cTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGc                                                                                                                  >  2:3490591/1‑150 (MQ=255)
                                    cTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGc                                                                                                                  >  2:2306650/1‑150 (MQ=255)
                                        aTGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCa                                                                                                              <  2:2911997/150‑1 (MQ=255)
                                         tGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAc                                                                                                             >  2:4949056/1‑150 (MQ=255)
                                          gTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAcc                                                                                                            <  2:3490975/150‑1 (MQ=255)
                                             gAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtc                                                                                                         >  2:2867658/1‑150 (MQ=255)
                                               aTTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGGTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtcgt                                                                                                       >  1:4707965/1‑150 (MQ=255)
                                                tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTa                                                                                                      <  2:4841477/150‑1 (MQ=255)
                                                tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTa                                                                                                      >  2:1252500/1‑150 (MQ=255)
                                                   aTGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAata                                                                                                   <  2:2195451/150‑1 (MQ=255)
                                                   aTGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAata                                                                                                   <  1:2130965/150‑1 (MQ=255)
                                                    tGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAatat                                                                                                  <  2:2024973/150‑1 (MQ=255)
                                                     gctggctgTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATAtt                                                                                                 <  1:4394747/150‑1 (MQ=255)
                                                     gctggctgTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATAtt                                                                                                 <  2:3438226/150‑1 (MQ=255)
                                                       tggctgTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTgc                                                                                               <  2:3911182/150‑1 (MQ=255)
                                                             tCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAgg                                                                                         >  1:2483080/1‑150 (MQ=255)
                                                                 gctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGt                                                                                     >  2:1285049/1‑150 (MQ=255)
                                                                 gctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGt                                                                                     <  1:1025289/150‑1 (MQ=255)
                                                                   tggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGACGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTg                                                                                   <  2:5079161/150‑1 (MQ=255)
                                                                    ggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGa                                                                                  <  1:1204822/150‑1 (MQ=255)
                                                                     gAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGaa                                                                                 >  1:4520545/1‑150 (MQ=255)
                                                                        cTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGgt                                                                              >  1:944161/1‑150 (MQ=255)
                                                                         tGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGgtg                                                                             <  2:3129399/150‑1 (MQ=255)
                                                                          ggTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGgtgt                                                                            <  2:2089229/150‑1 (MQ=255)
                                                                              gCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCAc                                                                        <  1:3279589/150‑1 (MQ=255)
                                                                               cccAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCAcc                                                                       >  1:903425/1‑150 (MQ=255)
                                                                                ccAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCt                                                                      <  1:351197/150‑1 (MQ=255)
                                                                                   ggCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCAt                                                                   <  1:4774562/150‑1 (MQ=255)
                                                                                   ggCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCAt                                                                   <  1:3868739/150‑1 (MQ=255)
                                                                                     cTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATcc                                                                 >  1:1298738/1‑150 (MQ=255)
                                                                                     cTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATcc                                                                 <  2:4026789/150‑1 (MQ=255)
                                                                                      tCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCg                                                                >  1:1143838/1‑150 (MQ=255)
                                                                                      tCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCg                                                                >  1:2920642/1‑150 (MQ=255)
                                                                                      tCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCg                                                                >  1:723168/1‑150 (MQ=255)
                                                                                       cGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCCTGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCgg                                                               >  2:3947921/1‑150 (MQ=255)
                                                                                       cGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCgg                                                               >  1:2672722/1‑150 (MQ=255)
                                                                                        gCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGgc                                                              <  1:150503/150‑1 (MQ=255)
                                                                                              ttGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCatat                                                        >  1:4012119/1‑150 (MQ=255)
                                                                                               tGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATAtt                                                       <  2:1455142/150‑1 (MQ=255)
                                                                                               tGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATAtt                                                       <  2:2601058/150‑1 (MQ=255)
                                                                                                gCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATAttt                                                      <  2:410059/150‑1 (MQ=255)
                                                                                                 cgacgaCCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTg                                                     >  2:2924987/1‑150 (MQ=255)
                                                                                                 cgacgaCCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTg                                                     <  2:1145647/150‑1 (MQ=255)
                                                                                                 cgacgaCCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCAGATTTg                                                     >  2:1515145/1‑150 (MQ=255)
                                                                                                     gacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCttt                                                 >  2:3042333/1‑150 (MQ=255)
                                                                                                      acCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCtttt                                                <  1:2614735/150‑1 (MQ=255)
                                                                                                      acCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCtttt                                                <  1:1350340/150‑1 (MQ=255)
                                                                                                        cccAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCa                                              >  2:1713877/1‑150 (MQ=255)
                                                                                                         ccAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAg                                             >  2:4922097/1‑150 (MQ=255)
                                                                                                          cAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGc                                            <  1:1184750/150‑1 (MQ=255)
                                                                                                           aGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCa                                           >  2:4720567/1‑150 (MQ=255)
                                                                                                               gCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAg                                       <  2:1805434/150‑1 (MQ=255)
                                                                                                                cGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAgg                                      >  1:3035234/1‑150 (MQ=255)
                                                                                                                 ggTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgt                                     >  2:3857295/1‑150 (MQ=255)
                                                                                                                   tCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCAGGCGCATATTTGCTTTTCAGCATCAGgtcg                                   >  2:2364937/1‑150 (MQ=255)
                                                                                                                    cATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtcgt                                  <  1:5082825/150‑1 (MQ=255)
                                                                                                                    cATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtcgt                                  >  1:2161612/1‑150 (MQ=255)
                                                                                                                    cATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtcgt                                  <  1:492854/150‑1 (MQ=255)
                                                                                                                    cATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtcgt                                  <  1:3102030/150‑1 (MQ=255)
                                                                                                                        cgggcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                              <  1:1490093/150‑1 (MQ=255)
                                                                                                                        cgggcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                              <  1:1714655/150‑1 (MQ=255)
                                                                                                                            cgggTGGATTGCCAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCgg                          >  1:133788/1‑150 (MQ=255)
                                                                                                                               gTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTg                       <  1:5192480/150‑1 (MQ=255)
                                                                                                                                 ggATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGgc                     >  2:1930854/1‑150 (MQ=255)
                                                                                                                                      gCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTGTTCAGCATCAGGTCGTCGAGGCGGGTGGCGCtgt                >  2:3223470/1‑150 (MQ=255)
                                                                                                                                       ccAGCGATCGGCGTAGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCtgtg               <  2:3773341/150‑1 (MQ=255)
                                                                                                                                        cAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGc              <  1:185735/150‑1 (MQ=255)
                                                                                                                                         aGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCt             <  2:241122/150‑1 (MQ=255)
                                                                                                                                         aGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTGTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCt             >  2:2825513/1‑150 (MQ=255)
                                                                                                                                             aTCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcgcg         >  1:5158410/1‑150 (MQ=255)
                                                                                                                                             aTCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcgcg         >  2:2149504/1‑150 (MQ=255)
                                                                                                                                                ggCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCa      >  2:750122/1‑150 (MQ=255)
                                                                                                                                                ggCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGACGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCa      >  2:4050941/1‑150 (MQ=255)
                                                                                                                                                 gCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAg     <  1:3119368/150‑1 (MQ=255)
                                                                                                                                                  cgtcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGt    >  2:1436091/1‑150 (MQ=255)
                                                                                                                                                    tcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCa  <  1:781414/150‑1 (MQ=255)
                                                                                                                                                    tcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCa  <  1:1664425/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCA  >  CP009273/4505289‑4505586

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