Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 706,097 A→C 100% T2T (ACT→ACG
*29G (TGA→GGA) 
fur ←
uof ←
ferric iron uptake regulon transcriptional repressor; autorepressor
ryhB‑regulated fur leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273706,0970AC98.0% 1320.8 / 2.2 339T2T (ACT→ACG
*29G (TGA→GGA) 
fur
uof
ferric iron uptake regulon transcriptional repressor; autorepressor
ryhB‑regulated fur leader peptide
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (173/157);  total (177/162)
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
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 339 total aligned reads displayed.

ACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTG  >  CP009273/705950‑706245
                                                                                                                                                   |                                                                                                                                                    
aCCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGt                                                                                                                                                    >  1:471669/1‑150 (MQ=255)
  cAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCa                                                                                                                                                  <  2:2151764/150‑1 (MQ=255)
    agttcttcACCCATAGCGATCTGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATg                                                                                                                                                <  1:1422987/148‑1 (MQ=255)
           cACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATc                                                                                                                                         >  2:2956703/1‑150 (MQ=255)
            aCCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCt                                                                                                                                        >  1:3219643/1‑150 (MQ=255)
             cccATCTCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTg                                                                                                                                       <  2:3709391/150‑1 (MQ=255)
             cccATATCGATCAGACGTTTGTATAACTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTg                                                                                                                                       <  2:981761/150‑1 (MQ=255)
              ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGt                                                                                                                                      >  2:3872060/1‑150 (MQ=255)
                atatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTcc                                                                                                                                    >  2:5043385/1‑150 (MQ=255)
                atatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTcc                                                                                                                                    <  2:1986361/150‑1 (MQ=255)
                   tCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGt                                                                                                                                 <  2:3822510/150‑1 (MQ=255)
                   tCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGt                                                                                                                                 >  1:2825724/1‑150 (MQ=255)
                    cGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGtt                                                                                                                                >  2:1017274/1‑150 (MQ=255)
                      aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTAc                                                                                                                              >  2:3257049/1‑150 (MQ=255)
                      aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTAc                                                                                                                              >  2:2566991/1‑150 (MQ=255)
                       tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACt                                                                                                                             >  1:1199464/1‑150 (MQ=255)
                         aGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTaa                                                                                                                           >  2:2264945/1‑150 (MQ=255)
                         aGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTaa                                                                                                                           >  1:3614194/1‑150 (MQ=255)
                          gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAg                                                                                                                          >  2:429627/1‑150 (MQ=255)
                          gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAg                                                                                                                          >  2:2810838/1‑150 (MQ=255)
                             gTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTCCTAAGCgg                                                                                                                       >  2:1301282/1‑150 (MQ=255)
                                tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                    >  1:3027355/1‑150 (MQ=255)
                                 gTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCa                                                                                                                   >  2:1748149/1‑150 (MQ=255)
                                 gTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCa                                                                                                                   >  2:2757057/1‑150 (MQ=255)
                                 gTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCa                                                                                                                   >  1:4178427/1‑150 (MQ=255)
                                     aaaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTc                                                                                                               <  2:3732830/150‑1 (MQ=255)
                                        tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                                                                            <  2:5122616/150‑1 (MQ=255)
                                        tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                                                                            <  1:807859/150‑1 (MQ=255)
                                              gCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaaga                                                                                                      <  2:3340844/150‑1 (MQ=255)
                                               cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                                                                                     <  1:2409023/150‑1 (MQ=255)
                                                aCTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGaagaag                                                                                                    <  1:4416642/150‑1 (MQ=255)
                                                 cTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                                                                                   >  1:717717/1‑150 (MQ=255)
                                                  tGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                                                                  <  1:2632427/150‑1 (MQ=255)
                                                  tGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                                                                  <  2:2453971/150‑1 (MQ=255)
                                                   gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                                                                 <  1:3060643/150‑1 (MQ=255)
                                                   gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                                                                 <  1:4577456/150‑1 (MQ=255)
                                                   gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                                                                 >  2:4765018/1‑150 (MQ=255)
                                                    aCGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGAc                                                                                                <  1:1022332/150‑1 (MQ=255)
                                                      gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                                                              >  1:4333388/1‑150 (MQ=255)
                                                      gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                                                              >  1:3690750/1‑150 (MQ=255)
                                                      gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                                                              >  1:2884661/1‑150 (MQ=255)
                                                       tGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                                                             >  1:1681229/1‑150 (MQ=255)
                                                       tGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                                                             >  1:3332083/1‑150 (MQ=255)
                                                        gATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGaa                                                                                            >  1:1954422/1‑150 (MQ=255)
                                                          tGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                                                                          <  1:3205596/150‑1 (MQ=255)
                                                           ggTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAttt                                                                                         <  1:2498493/150‑1 (MQ=255)
                                                           ggTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAttt                                                                                         <  1:2653453/150‑1 (MQ=255)
                                                            gTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTg                                                                                        <  1:202220/150‑1 (MQ=255)
                                                             ttGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                                                                       <  1:3642119/150‑1 (MQ=255)
                                                                  cGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCTCTTCATTAGAAGAAGTGACAGAATTTGCTCTTg                                                                                  >  1:4649030/1‑150 (MQ=255)
                                                                     cTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgaga                                                                               >  2:2922214/1‑150 (MQ=255)
                                                                      tCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGAt                                                                              <  1:5028435/150‑1 (MQ=255)
                                                                       ccTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                                                                             <  2:4470411/150‑1 (MQ=255)
                                                                       ccTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                                                                             <  2:1812011/150‑1 (MQ=255)
                                                                        cTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                                                            >  1:2235159/1‑150 (MQ=255)
                                                                        cTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                                                            >  1:1637167/1‑150 (MQ=255)
                                                                         tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                           <  1:601669/150‑1 (MQ=255)
                                                                           aagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                                                                         <  2:815372/150‑1 (MQ=255)
                                                                               aCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTAt                                                                     >  2:1534790/1‑150 (MQ=255)
                                                                                 ttCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCa                                                                   >  1:4999771/1‑150 (MQ=255)
                                                                                   ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAtt                                                                 <  2:782863/150‑1 (MQ=255)
                                                                                       gATTTTTACACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                                                             <  2:4549488/150‑1 (MQ=255)
                                                                                       gATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                                                             >  2:524411/1‑150 (MQ=255)
                                                                                       gATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                                                             <  1:4165010/150‑1 (MQ=255)
                                                                                           tttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt                                                         >  1:3904641/1‑150 (MQ=255)
                                                                                           tttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt                                                         >  2:1479708/1‑150 (MQ=255)
                                                                                            ttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg                                                        >  1:3467083/1‑150 (MQ=255)
                                                                                             tAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGc                                                       <  2:3550871/150‑1 (MQ=255)
                                                                                                aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc                                                    <  1:3819016/150‑1 (MQ=255)
                                                                                                      aaGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa                                              >  2:4631856/1‑150 (MQ=255)
                                                                                                      aaGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa                                              <  2:3614194/150‑1 (MQ=255)
                                                                                                        gCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAAt                                            <  1:1606044/150‑1 (MQ=255)
                                                                                                          gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                          >  2:4806701/1‑150 (MQ=255)
                                                                                                          gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                          <  1:525950/150‑1 (MQ=255)
                                                                                                             aCTATCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTTACCTGTTACTAATCGGTTCACTACATTAGAAGAAGTGACAGAATTAGCTCTTGAGATAATGCGTATCAATATATAAATGCCACGCCTAAATGaaaa                                       >  2:3852260/1‑150 (MQ=255)
                                                                                                               tttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAAcc                                     >  1:1906996/1‑150 (MQ=255)
                                                                                                                  cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcaca                                  <  1:5181912/150‑1 (MQ=255)
                                                                                                                  cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcaca                                  <  2:647772/150‑1 (MQ=255)
                                                                                                                  cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcaca                                  <  2:4154265/150‑1 (MQ=255)
                                                                                                                  cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcaca                                  >  2:3491101/1‑150 (MQ=255)
                                                                                                                   aGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACaa                                 <  2:1221063/150‑1 (MQ=255)
                                                                                                                     gCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGt                               <  2:1105213/150‑1 (MQ=255)
                                                                                                                        aGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAATCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGATCTTTAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATTAAAACCACAATTccc                            >  1:3494802/1‑150 (MQ=255)
                                                                                                                        aGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTccc                            <  2:4050404/150‑1 (MQ=255)
                                                                                                                        aGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTccc                            >  1:965975/1‑150 (MQ=255)
                                                                                                                         gCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGTTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                           <  2:4316284/150‑1 (MQ=255)
                                                                                                                         gCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                           >  2:4071959/1‑150 (MQ=255)
                                                                                                                           tttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                         <  1:2548358/150‑1 (MQ=255)
                                                                                                                            ttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                        <  1:2073965/150‑1 (MQ=255)
                                                                                                                                  tgggCGGTATTGTTATCCGTCATGGGGTATCTGTCCTGTTACTTAGCGGTTCACTTCATTAGAAGTTGTGACAGATTTTGCTCTTGAGTTAATGCTTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                  <  2:4945413/149‑1 (MQ=255)
                                                                                                                                      cGGTTTTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGATGTGACAGAATTTGCTCTTGAGATAATGCTTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGtt              <  1:2330132/150‑1 (MQ=255)
                                                                                                                                       ggTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta             >  1:4442064/1‑150 (MQ=255)
                                                                                                                                       ggTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta             >  2:2491494/1‑150 (MQ=255)
                                                                                                                                        gTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa            >  2:4853060/1‑150 (MQ=255)
                                                                                                                                        gTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa            >  2:183279/1‑150 (MQ=255)
                                                                                                                                         tATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaat           >  1:1387914/1‑150 (MQ=255)
                                                                                                                                             gTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaaa       >  1:1259181/1‑150 (MQ=255)
                                                                                                                                               tATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAAc     >  1:1346917/1‑150 (MQ=255)
                                                                                                                                                aTCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACg    >  2:5322679/1‑150 (MQ=255)
                                                                                                                                                aTCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACg    >  2:2683451/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTG  >  CP009273/705950‑706245

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: