Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,505,471 C→T 6.6% D342N (GAC→AAC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,505,4710CT 6.6% 1080.2 / 41.8 302D342N (GAC→AAC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (127/154);  new base T (9/11);  total (137/165)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

CCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGT  >  CP009273/4505324‑4505618
                                                                                                                                                   |                                                                                                                                                   
ccTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCg                                                                                                                                                    >  2:372686/1‑149 (MQ=255)
 cTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCgg                                                                                                                                                   >  2:1874186/1‑149 (MQ=255)
    aaTGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCtt                                                                                                                                                >  1:2064656/1‑149 (MQ=255)
      tGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCa                                                                                                                                              >  1:1015528/1‑149 (MQ=255)
       gTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAc                                                                                                                                             >  1:3628321/1‑149 (MQ=255)
        ttGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAcc                                                                                                                                            >  2:2240415/1‑149 (MQ=255)
        ttGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAcc                                                                                                                                            >  1:496603/1‑149 (MQ=255)
        ttGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAcc                                                                                                                                            >  1:5053000/1‑149 (MQ=255)
         tgtaTTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCTGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAcct                                                                                                                                           <  1:2651509/146‑2 (MQ=255)
          gAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgt                                                                                                                                          >  1:1623157/1‑149 (MQ=255)
            aTTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtcg                                                                                                                                        <  1:198866/149‑1 (MQ=255)
             tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACcgtcgt                                                                                                                                       <  2:260816/149‑1 (MQ=255)
             tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtcgt                                                                                                                                       >  2:905691/1‑149 (MQ=255)
             tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtcgt                                                                                                                                       <  2:310903/149‑1 (MQ=255)
             tttATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcgtcgt                                                                                                                                       >  1:4786183/1‑149 (MQ=255)
               tATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTaa                                                                                                                                     <  2:3475021/149‑1 (MQ=255)
               tATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTaa                                                                                                                                     <  2:4192835/149‑1 (MQ=255)
                aTGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAat                                                                                                                                    >  1:942988/1‑149 (MQ=255)
                 tGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAata                                                                                                                                   <  1:4826467/149‑1 (MQ=255)
                          tCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGcagcag                                                                                                                          >  2:2497976/1‑149 (MQ=255)
                          tCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGcagcag                                                                                                                          <  2:2413748/149‑1 (MQ=255)
                          tCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGcagcag                                                                                                                          >  2:1107518/1‑149 (MQ=255)
                          tCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGcagcag                                                                                                                          >  2:2056180/1‑149 (MQ=255)
                            tggctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGc                                                                                                                        <  1:685401/149‑1 (MQ=255)
                            tggctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGc                                                                                                                        <  2:1646744/149‑1 (MQ=255)
                              gctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTg                                                                                                                      >  2:4148535/1‑149 (MQ=255)
                               ctggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGt                                                                                                                     <  2:2371626/149‑1 (MQ=255)
                                 ggAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTg                                                                                                                   >  2:1012634/1‑149 (MQ=255)
                                  gAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGa                                                                                                                  >  2:94492/1‑149 (MQ=255)
                                   aaCTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGaa                                                                                                                 <  2:5007779/149‑1 (MQ=255)
                                      tGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGgt                                                                                                              <  2:4034420/149‑1 (MQ=255)
                                       ggTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGgtg                                                                                                             >  2:1379824/1‑149 (MQ=255)
                                          aGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGc                                                                                                          <  1:2806903/149‑1 (MQ=255)
                                           gCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCa                                                                                                         >  1:4313333/1‑149 (MQ=255)
                                           gCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCa                                                                                                         >  1:2378942/1‑149 (MQ=255)
                                             ccAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCAcc                                                                                                       >  1:4403155/1‑149 (MQ=255)
                                              cAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCt                                                                                                      >  1:287386/1‑149 (MQ=255)
                                               aGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTc                                                                                                     >  2:3249031/1‑149 (MQ=255)
                                               aGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTc                                                                                                     <  2:4093314/149‑1 (MQ=255)
                                                ggCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCa                                                                                                    <  2:1903911/149‑1 (MQ=255)
                                                ggCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCa                                                                                                    <  2:4387315/149‑1 (MQ=255)
                                                 gCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCAt                                                                                                   <  1:3881565/149‑1 (MQ=255)
                                                 gCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCAt                                                                                                   <  2:3927891/149‑1 (MQ=255)
                                                  cTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATc                                                                                                  <  1:2607738/149‑1 (MQ=255)
                                                  cTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATc                                                                                                  <  1:2460274/149‑1 (MQ=255)
                                                  cTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATc                                                                                                  <  1:3905919/149‑1 (MQ=255)
                                                       cAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGgcg                                                                                             >  2:688719/1‑149 (MQ=255)
                                                         gCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCa                                                                                           <  2:90928/149‑1 (MQ=255)
                                                          cTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCat                                                                                          <  1:306335/149‑1 (MQ=255)
                                                             gCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATAtt                                                                                       <  2:4937662/149‑1 (MQ=255)
                                                             gCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATAtt                                                                                       <  2:2085558/149‑1 (MQ=255)
                                                               gacgacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTg                                                                                     >  1:638207/1‑149 (MQ=255)
                                                                acgacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGc                                                                                    <  1:1483168/149‑1 (MQ=255)
                                                                 cgacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCt                                                                                   <  2:3861609/149‑1 (MQ=255)
                                                                  gacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCTTCGTAATCCGCGAGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCtt                                                                                  <  2:1480191/149‑1 (MQ=255)
                                                                  gacCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCtt                                                                                  <  1:2204963/149‑1 (MQ=255)
                                                                   acCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCttt                                                                                 <  2:326948/149‑1 (MQ=255)
                                                                    ccccAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCtttt                                                                                <  2:4600770/149‑1 (MQ=255)
                                                                     cccAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTc                                                                               <  1:2249755/149‑1 (MQ=255)
                                                                       cAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAg                                                                             >  1:4301934/1‑149 (MQ=255)
                                                                        aGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGc                                                                            <  1:302514/149‑1 (MQ=255)
                                                                          aaGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCAt                                                                          <  2:1090443/149‑1 (MQ=255)
                                                                          aaGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCAt                                                                          <  2:2855926/149‑1 (MQ=255)
                                                                           aGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATc                                                                         >  1:1107138/1‑149 (MQ=255)
                                                                              ggTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAgg                                                                      <  2:2770135/149‑1 (MQ=255)
                                                                              ggTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAgg                                                                      >  1:3101913/1‑149 (MQ=255)
                                                                               gTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgt                                                                     >  2:1260872/1‑149 (MQ=255)
                                                                                tCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtc                                                                    <  2:2347086/149‑1 (MQ=255)
                                                                                   tACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgtcgtc                                                                 <  1:5143454/149‑1 (MQ=255)
                                                                                      gggcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                                                              <  2:5105669/149‑1 (MQ=255)
                                                                                      gggcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                                                              <  2:3371410/149‑1 (MQ=255)
                                                                                        gcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGc                                                            <  2:4786183/149‑1 (MQ=255)
                                                                                        gcgggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGc                                                            >  1:3791700/1‑149 (MQ=255)
                                                                                          gggTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGTTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCgg                                                          >  2:2031269/1‑149 (MQ=255)
                                                                                            gTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGt                                                        <  1:372686/149‑1 (MQ=255)
                                                                                              ggATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTgg                                                      <  1:4224602/149‑1 (MQ=255)
                                                                                                    ccAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCtgt                                                >  1:1343272/1‑149 (MQ=255)
                                                                                                     cAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCtgtg                                               >  1:1694345/1‑149 (MQ=255)
                                                                                                     cAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCtgtg                                               >  2:3608953/1‑149 (MQ=255)
                                                                                                      aGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGc                                              >  2:1315789/1‑149 (MQ=255)
                                                                                                       gCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCt                                             >  1:4016883/1‑149 (MQ=255)
                                                                                                         gATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcg                                           >  2:2578634/1‑149 (MQ=255)
                                                                                                         gATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcg                                           >  2:4921127/1‑149 (MQ=255)
                                                                                                          aTCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcgc                                          >  1:861498/1‑149 (MQ=255)
                                                                                                             ggCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGcc                                       <  1:3423947/149‑1 (MQ=255)
                                                                                                             ggCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGcc                                       <  2:2918048/149‑1 (MQ=255)
                                                                                                             ggCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGcc                                       >  1:1925530/1‑149 (MQ=255)
                                                                                                                gtcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGt                                    >  2:1534730/1‑149 (MQ=255)
                                                                                                                 tcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTc                                   <  1:1240211/149‑1 (MQ=255)
                                                                                                                 tcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTc                                   <  2:297916/149‑1 (MQ=255)
                                                                                                                 tcgtAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTc                                   >  2:4866129/1‑149 (MQ=255)
                                                                                                                      aTCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGc                              >  2:2018796/1‑149 (MQ=255)
                                                                                                                      aTCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGc                              <  2:1347917/149‑1 (MQ=255)
                                                                                                                       tCCGCGCGAGACAGGCCCCCGGGCATGTTGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCAGATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGcc                             >  2:386302/1‑149 (MQ=255)
                                                                                                                       tCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGcc                             <  2:3233695/149‑1 (MQ=255)
                                                                                                                                    ggCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAg                <  2:775622/149‑1 (MQ=255)
                                                                                                                                      ccACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTa              <  1:4070779/149‑1 (MQ=255)
                                                                                                                                          cGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGagc          <  2:1813630/149‑1 (MQ=255)
                                                                                                                                           gggCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGagca         >  1:1802430/1‑149 (MQ=255)
                                                                                                                                            ggCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGagcag        >  2:4672788/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGGT  >  CP009273/4505324‑4505618

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: