Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,502,777 G→A 5.1% N150N (AAC→AAT fecC ← iron‑dicitrate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,502,7770GA 5.1% 1362.3 / 37.1 373N150N (AAC→AATfecCiron‑dicitrate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (190/164);  new base A (11/8);  total (201/172)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 395 total aligned reads displayed.

AGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGC  >  CP009273/4502628‑4502926
                                                                                                                                                     |                                                                                                                                                     
agCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAgtttg                                                                                                                                                       >  2:524306/1‑150 (MQ=255)
 gCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAgtttgt                                                                                                                                                      <  2:5570625/150‑1 (MQ=255)
 gCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAgtttgt                                                                                                                                                      <  2:1417316/150‑1 (MQ=255)
    gCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTc                                                                                                                                                   <  1:4019142/150‑1 (MQ=255)
     ccAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCt                                                                                                                                                  <  1:4934084/150‑1 (MQ=255)
      cAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTg                                                                                                                                                 >  2:2730247/1‑150 (MQ=255)
        gACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTc                                                                                                                                               >  2:4237238/1‑150 (MQ=255)
          cATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCCGCCAGTAACAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCAt                                                                                                                                             <  2:197945/150‑1 (MQ=255)
          cATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCAt                                                                                                                                             >  2:2045795/1‑150 (MQ=255)
          cATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCAt                                                                                                                                             <  1:3021237/150‑1 (MQ=255)
           aTCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATg                                                                                                                                            <  1:5119070/150‑1 (MQ=255)
               tGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTa                                                                                                                                        <  2:2410432/150‑1 (MQ=255)
                 ccAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATg                                                                                                                                      <  2:4744405/150‑1 (MQ=255)
                    gCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACg                                                                                                                                   <  1:2936707/150‑1 (MQ=255)
                     cGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTGTGTTTCTGTCATGGGTATGACGa                                                                                                                                  >  2:3515671/1‑150 (MQ=255)
                      gggCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGaa                                                                                                                                 <  2:972265/150‑1 (MQ=255)
                      gggCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGaa                                                                                                                                 >  1:1047042/1‑150 (MQ=255)
                        gCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAAt                                                                                                                               <  1:778325/150‑1 (MQ=255)
                            gggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGTATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGc                                                                                                                           >  1:3215505/1‑150 (MQ=255)
                            gggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGAAATCCGc                                                                                                                           >  1:2243503/1‑150 (MQ=255)
                             ggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGcc                                                                                                                          <  2:409605/150‑1 (MQ=255)
                              gACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGcct                                                                                                                         <  2:217586/150‑1 (MQ=255)
                               acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGcctc                                                                                                                        >  2:1942701/1‑150 (MQ=255)
                                cacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTATGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGcctcc                                                                                                                       <  2:3467183/150‑1 (MQ=255)
                                  cTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTg                                                                                                                     >  1:4766167/1‑150 (MQ=255)
                                   tCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGc                                                                                                                    <  1:1049106/150‑1 (MQ=255)
                                      gCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGt                                                                                                                 <  1:2338475/150‑1 (MQ=255)
                                           ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGa                                                                                                            >  1:5578205/1‑150 (MQ=255)
                                             agccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGAcc                                                                                                          >  1:2536030/1‑150 (MQ=255)
                                             agccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATACGCCTCCTGCGGTCATGAcc                                                                                                          >  2:1159938/1‑150 (MQ=255)
                                               ccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACcag                                                                                                        >  2:1335370/1‑150 (MQ=255)
                                               ccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACcag                                                                                                        >  1:1240977/1‑150 (MQ=255)
                                                cagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACcagc                                                                                                       >  1:3360836/1‑150 (MQ=255)
                                                 agTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTCTGACGAAATCCGCCTCCTGCGGTCATGACcagca                                                                                                      >  2:2962806/1‑150 (MQ=255)
                                                   tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCagcagc                                                                                                    >  2:4409053/1‑150 (MQ=255)
                                                    aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGcagcc                                                                                                   >  1:3076098/1‑150 (MQ=255)
                                                    aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGcagcc                                                                                                   >  2:2493711/1‑150 (MQ=255)
                                                       aGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGcagccagc                                                                                                >  2:3190420/1‑150 (MQ=255)
                                                         aTGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTc                                                                                              >  1:626941/1‑150 (MQ=255)
                                                         aTGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTc                                                                                              <  2:5494197/150‑1 (MQ=255)
                                                              gTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCAcgcc                                                                                         >  1:1909895/1‑150 (MQ=255)
                                                               tAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCAcgccg                                                                                        >  2:4468360/1‑150 (MQ=255)
                                                                   cATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGcccc                                                                                    >  1:3303214/1‑150 (MQ=255)
                                                                    aTGTTCTTCGGCCAGCAGCAGGGTGATGAGGGTAAGGCCCATACAAAAGGCCGACAGCTCGATTCCCGCGAGGATCAGTTTGTTTGTGTCATGGGTATGTCGAAATCCGCCTCCTGCGGTCATGACCAGCAGACAGCTCACGCCGccccc                                                                                   <  2:2771510/150‑1 (MQ=255)
                                                                      gATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGc                                                                                 >  2:2372781/1‑150 (MQ=255)
                                                                       aTCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCa                                                                                >  2:3398435/1‑150 (MQ=255)
                                                                       aTCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCa                                                                                >  2:736631/1‑150 (MQ=255)
                                                                        tCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCAt                                                                               <  1:1252752/150‑1 (MQ=255)
                                                                         cTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATg                                                                              >  2:3169518/1‑150 (MQ=255)
                                                                           tCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGTAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCTCGCCGCCCCCGCATGcc                                                                            <  1:1679750/150‑1 (MQ=255)
                                                                           tCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGcc                                                                            <  1:1209264/150‑1 (MQ=255)
                                                                            cGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCcg                                                                           <  2:1635943/150‑1 (MQ=255)
                                                                             ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCcgc                                                                          >  1:4674091/1‑150 (MQ=255)
                                                                             ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCcgc                                                                          >  2:1588640/1‑150 (MQ=255)
                                                                             ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCcgc                                                                          >  1:138573/1‑150 (MQ=255)
                                                                              gCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCcgcg                                                                         <  2:1017263/150‑1 (MQ=255)
                                                                               ccAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGa                                                                        >  2:4937881/1‑150 (MQ=255)
                                                                                 agcagcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatg                                                                      <  1:2714595/150‑1 (MQ=255)
                                                                                 agcagcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatg                                                                      <  2:2680567/150‑1 (MQ=255)
                                                                                   cagcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatgaa                                                                    <  1:4592443/150‑1 (MQ=255)
                                                                                    agcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatgaat                                                                   <  1:4204310/150‑1 (MQ=255)
                                                                                     gcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGatgaatg                                                                  <  1:5665191/150‑1 (MQ=255)
                                                                                       agGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGAc                                                                >  2:2558841/1‑150 (MQ=255)
                                                                                        gggTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACa                                                               >  2:328880/1‑150 (MQ=255)
                                                                                          gTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACagt                                                             <  2:5367098/150‑2 (MQ=255)
                                                                                          gTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACaga                                                             >  1:4391420/1‑150 (MQ=255)
                                                                                             aTGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAgagac                                                          >  2:600816/1‑149 (MQ=255)
                                                                                                cGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAg                                                       <  1:98151/150‑1 (MQ=255)
                                                                                                cGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAg                                                       <  1:1785469/150‑1 (MQ=255)
                                                                                                  ggTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGcc                                                     <  2:1124341/150‑1 (MQ=255)
                                                                                                  ggTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGcc                                                     >  2:5501165/1‑150 (MQ=255)
                                                                                                    tAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTg                                                   <  2:5105884/150‑1 (MQ=255)
                                                                                                     aaGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGc                                                  >  2:1000422/1‑150 (MQ=255)
                                                                                                      aGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCa                                                 <  2:3852483/150‑1 (MQ=255)
                                                                                                        gCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAAt                                               >  2:2274668/1‑150 (MQ=255)
                                                                                                          ccctACCAAACGCCGACAGCGCGATACCCGCGAGGACCAGTTTGTTTCTGTCACGGGTATGACGAAATCCGCCTCCGGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCg                                             <  2:4853782/147‑1 (MQ=255)
                                                                                                            aTACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGc                                           >  1:1804558/1‑150 (MQ=255)
                                                                                                            aTACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGc                                           <  1:366154/150‑1 (MQ=255)
                                                                                                              aCAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGt                                         >  2:1529363/1‑150 (MQ=255)
                                                                                                               cAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTc                                        <  1:4979025/150‑1 (MQ=255)
                                                                                                                 aaaGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCAGGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCgg                                      >  1:1967851/1‑150 (MQ=255)
                                                                                                                  aaGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGa                                     >  1:4268294/1‑150 (MQ=255)
                                                                                                                        gACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAgc                               >  2:5077571/1‑150 (MQ=255)
                                                                                                                        gACAGCGCGATACCCGCGAGGATCAGTTTATTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACCGAGAATAGCCTGCCATCGGCGTCGGACTCAgc                               >  2:3636074/1‑150 (MQ=255)
                                                                                                                          cAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAgcgc                             >  1:372950/1‑150 (MQ=255)
                                                                                                                          cAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAgcgc                             >  2:4498227/1‑150 (MQ=255)
                                                                                                                             cgcgATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCt                          <  2:400419/150‑1 (MQ=255)
                                                                                                                                     ccGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAgcg                  <  1:3819196/150‑1 (MQ=255)
                                                                                                                                      cgcgAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAgcgc                 >  1:881817/1‑150 (MQ=255)
                                                                                                                                       gcgAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGcc                >  2:2115111/1‑150 (MQ=255)
                                                                                                                                        cgAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCa               <  1:5595153/150‑1 (MQ=255)
                                                                                                                                        cgAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGACTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGAAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCa               >  1:2968545/1‑150 (MQ=255)
                                                                                                                                         gAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCAt              >  2:3774823/1‑150 (MQ=255)
                                                                                                                                          aGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATa             >  2:1053936/1‑150 (MQ=255)
                                                                                                                                          aGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAAACGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATa             >  2:3721665/1‑150 (MQ=255)
                                                                                                                                             aTCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGcc          <  1:2634064/150‑1 (MQ=255)
                                                                                                                                               cAGTTTATTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAg        >  2:1764462/1‑150 (MQ=255)
                                                                                                                                                   ttgtttCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAgcgcc    <  2:4528720/150‑1 (MQ=255)
                                                                                                                                                   ttgtttCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAgcgcc    <  2:1249041/150‑1 (MQ=255)
                                                                                                                                                     gtttCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAgcgccgc  <  1:3502605/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGC  >  CP009273/4502628‑4502926

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: