Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 619,048 C→A 100% A307S (GCC→TCC)  fepB ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273619,0480CA98.6% 930.9 / 2.5 289A307S (GCC→TCC) fepBiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (117/168);  total (117/172)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 290 total aligned reads displayed.

GAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAT  >  CP009273/618901‑619194
                                                                                                                                                   |                                                                                                                                                  
gAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGgagc                                                                                                                                                  <  1:4675059/150‑1 (MQ=255)
 aGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCt                                                                                                                                                 >  2:2901978/1‑150 (MQ=255)
 aGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCt                                                                                                                                                 <  1:1769809/150‑1 (MQ=255)
  gCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGAACTTGCATGGAGCTg                                                                                                                                                <  1:1500802/150‑1 (MQ=255)
   cGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgt                                                                                                                                               >  2:2052816/1‑150 (MQ=255)
    ggTTTTGGTTTGTTGATTCTCGGCGTGTTGTTAGTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                              <  1:4479706/150‑1 (MQ=255)
    ggTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                              <  2:4360150/150‑1 (MQ=255)
    ggTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                              <  1:4892846/150‑1 (MQ=255)
    ggTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                              <  1:937870/150‑1 (MQ=255)
      ttttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtagt                                                                                                                                            >  2:3418321/1‑150 (MQ=255)
        ttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTaa                                                                                                                                          <  1:2899250/150‑1 (MQ=255)
         tGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAAt                                                                                                                                         <  2:1606679/150‑1 (MQ=255)
         tGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAAt                                                                                                                                         <  1:2234123/150‑1 (MQ=255)
         tGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATACAGCACTTGCATGGAGCTGTAGTAAt                                                                                                                                         >  1:1136523/1‑150 (MQ=255)
          ggTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATc                                                                                                                                        >  1:3912035/1‑150 (MQ=255)
             ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAg                                                                                                                                     <  2:5205687/150‑1 (MQ=255)
              tgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGa                                                                                                                                    <  2:1162702/150‑1 (MQ=255)
                 tgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACgg                                                                                                                                 >  2:3271739/1‑150 (MQ=255)
                 tgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACgg                                                                                                                                 >  2:907091/1‑150 (MQ=255)
                    ttATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACTTTTTCGCCAGTCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCTCTTGGATGGAGCTGTAGTAATCCAGACGGAAc                                                                                                                              <  1:2938817/150‑1 (MQ=255)
                    ttATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAAc                                                                                                                              >  1:4812651/1‑150 (MQ=255)
                      aTCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                            >  1:1316446/1‑150 (MQ=255)
                        cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtct                                                                                                                          >  2:3069629/1‑150 (MQ=255)
                         ggCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtctc                                                                                                                         <  1:2986044/150‑1 (MQ=255)
                           tgtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCgg                                                                                                                       >  2:3181628/2‑150 (MQ=255)
                            gtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGACCGACAGATAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGt                                                                                                                      >  2:2679667/1‑150 (MQ=255)
                               ttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTcc                                                                                                                   >  1:3813706/1‑150 (MQ=255)
                                tgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTccc                                                                                                                  <  1:2586109/150‑1 (MQ=255)
                                 gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                                 <  1:2304780/150‑1 (MQ=255)
                                 gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                                 >  2:3981334/1‑150 (MQ=255)
                                    aTTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAgcg                                                                                                              >  2:677145/1‑150 (MQ=255)
                                      tgctggtgctggTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCa                                                                                                            <  1:4213342/150‑1 (MQ=255)
                                       gctggtgctggtGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCAt                                                                                                           <  2:4279552/150‑1 (MQ=255)
                                        ctggagcaggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTAACCAGCACTTGAATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATa                                                                                                          >  2:271567/1‑150 (MQ=255)
                                         tggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                         <  1:1851242/150‑1 (MQ=255)
                                           gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAAc                                                                                                       <  2:2499303/150‑1 (MQ=255)
                                           gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAAc                                                                                                       <  2:1960851/150‑1 (MQ=255)
                                              ctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTg                                                                                                    >  2:178215/1‑150 (MQ=255)
                                                 gtggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCtt                                                                                                 <  2:2328292/150‑1 (MQ=255)
                                                      gTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGtttt                                                                                            <  1:3112901/150‑1 (MQ=255)
                                                        tGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGt                                                                                          <  1:906914/150‑1 (MQ=255)
                                                          cGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTAc                                                                                        <  2:3502143/150‑1 (MQ=255)
                                                           gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACt                                                                                       <  2:3816367/150‑1 (MQ=255)
                                                               ttttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgcag                                                                                   >  1:4702223/1‑150 (MQ=255)
                                                                 ttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgcaggc                                                                                 <  2:4794043/150‑1 (MQ=255)
                                                                     ccAGACGCCGCCGCAGGTGACAGCGTCCGACATTTAATGCTTAAAACAGCGACTTAAGCCTATCCAGCACTTGCATTGAGATGTAGTAATCCAGACGGAACGTCTCGGTTCCCATCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                             <  1:280756/150‑1 (MQ=255)
                                                                     ccAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                             <  2:1253547/150‑1 (MQ=255)
                                                                        gACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgcg                                                                          <  2:4885792/150‑1 (MQ=255)
                                                                            ccgccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGag                                                                      <  2:5060909/150‑1 (MQ=255)
                                                                              gccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagca                                                                    <  2:1531241/150‑1 (MQ=255)
                                                                               ccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagcag                                                                   <  2:2492676/150‑1 (MQ=255)
                                                                                agcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCATCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCTagctgc                                                                  <  1:116514/149‑1 (MQ=255)
                                                                                   aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                               >  2:885537/1‑150 (MQ=255)
                                                                                   aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                               <  2:842848/150‑1 (MQ=255)
                                                                                    ggTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAt                                                              <  2:4991081/150‑1 (MQ=255)
                                                                                      tGACAGCTTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                            <  2:1473032/150‑1 (MQ=255)
                                                                                      tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                            <  2:2253945/150‑1 (MQ=255)
                                                                                      tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                            >  1:4771012/1‑150 (MQ=255)
                                                                                      tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                            <  1:2188009/150‑1 (MQ=255)
                                                                                       gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAg                                                           >  1:334011/1‑150 (MQ=255)
                                                                                       gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAg                                                           >  1:5277845/1‑150 (MQ=255)
                                                                                        aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGc                                                          <  2:3765045/150‑1 (MQ=255)
                                                                                          aGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCAt                                                        <  1:27565/150‑1 (MQ=255)
                                                                                           gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATa                                                       <  1:3611573/150‑1 (MQ=255)
                                                                                             gTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTGGCATaaa                                                     >  2:423499/1‑150 (MQ=255)
                                                                                              tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAt                                                    <  2:2090383/150‑1 (MQ=255)
                                                                                              tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAt                                                    <  2:2209139/150‑1 (MQ=255)
                                                                                                cGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAg                                                  <  1:4948256/150‑1 (MQ=255)
                                                                                                    aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                              <  2:3100995/150‑1 (MQ=255)
                                                                                                    aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                              <  1:844595/150‑1 (MQ=255)
                                                                                                    aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                              >  2:1102927/1‑150 (MQ=255)
                                                                                                      ttAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTGTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCgg                                            >  2:4753117/1‑150 (MQ=255)
                                                                                                       tAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGc                                           >  2:550705/1‑150 (MQ=255)
                                                                                                        aaTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                          <  2:2523299/150‑1 (MQ=255)
                                                                                                          tGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATc                                        <  2:4125096/150‑1 (MQ=255)
                                                                                                           gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                       >  2:4658192/1‑150 (MQ=255)
                                                                                                           gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                       <  1:3467597/150‑1 (MQ=255)
                                                                                                               aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                   >  1:4121375/1‑150 (MQ=255)
                                                                                                                aaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTg                                  <  1:1191755/150‑1 (MQ=255)
                                                                                                                  aCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGAt                                <  1:2564567/150‑1 (MQ=255)
                                                                                                                    aGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCa                              >  1:1006631/1‑150 (MQ=255)
                                                                                                                           tAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGa                       >  2:3867231/1‑150 (MQ=255)
                                                                                                                           tAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGa                       <  1:4363511/150‑1 (MQ=255)
                                                                                                                            aaGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGaa                      <  2:135265/150‑1 (MQ=255)
                                                                                                                               ccTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAg                   >  1:3537680/1‑150 (MQ=255)
                                                                                                                                cTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAgg                  >  1:109773/1‑150 (MQ=255)
                                                                                                                                 tATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGAGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGa                 >  1:1969531/1‑150 (MQ=255)
                                                                                                                                     cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAg             <  1:5100450/150‑1 (MQ=255)
                                                                                                                                     cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAg             <  1:73175/150‑1 (MQ=255)
                                                                                                                                     cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAg             <  1:3181628/150‑1 (MQ=255)
                                                                                                                                      aGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGt            >  2:4549770/1‑150 (MQ=255)
                                                                                                                                        cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa          >  1:1671872/1‑150 (MQ=255)
                                                                                                                                         aCTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAc         >  1:237560/1‑150 (MQ=255)
                                                                                                                                          cTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACt        >  1:1887890/1‑150 (MQ=255)
                                                                                                                                           ttGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTc       <  2:4182817/150‑1 (MQ=255)
                                                                                                                                           ttGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTc       <  1:4027832/150‑1 (MQ=255)
                                                                                                                                            tGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCg      <  2:1745619/150‑1 (MQ=255)
                                                                                                                                              cATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc    >  2:1295816/1‑150 (MQ=255)
                                                                                                                                               aTGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCa   <  2:3912035/150‑1 (MQ=255)
                                                                                                                                                tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAt  <  1:3271739/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAT  >  CP009273/618901‑619194

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: