Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 706,097 A→C 73.4% 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,0970AC73.4% 589.9 / 182.7 270T2T (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 (43/29);  new base C (79/119);  total (122/148)
Fisher's exact test for biased strand distribution p-value = 5.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 271 total aligned reads displayed.

ACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGG  >  CP009273/705950‑706246
                                                                                                                                                   |                                                                                                                                                     
aCCAATTTCTTCACCCATATCGATCAGACGTTTGTATATATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGt                                                                                                                                                     <  2:1789289/150‑1 (MQ=255)
   aaTTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCAt                                                                                                                                                  >  1:4115540/1‑150 (MQ=255)
        cttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGa                                                                                                                                             >  1:923389/1‑150 (MQ=255)
        cttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGa                                                                                                                                             <  2:571383/150‑1 (MQ=255)
         ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGaa                                                                                                                                            <  1:1246120/150‑1 (MQ=255)
         ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGaa                                                                                                                                            >  2:3238463/1‑150 (MQ=255)
         ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGaa                                                                                                                                            >  2:2357708/1‑150 (MQ=255)
          tcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAAt                                                                                                                                           <  1:690127/150‑1 (MQ=255)
          tcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAAt                                                                                                                                           <  2:3940642/150‑1 (MQ=255)
             cccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTg                                                                                                                                        >  1:2558781/1‑150 (MQ=255)
                 tatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCt                                                                                                                                    <  1:1496939/150‑1 (MQ=255)
                     gATCAGACGTTTGTATAACTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTa                                                                                                                                <  2:559397/150‑1 (MQ=255)
                     gATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTa                                                                                                                                <  2:2339645/150‑1 (MQ=255)
                     gATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTa                                                                                                                                >  2:3580953/1‑150 (MQ=255)
                      aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTAc                                                                                                                               >  2:3275709/1‑150 (MQ=255)
                        cAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                                                                                                                             <  1:1007350/150‑1 (MQ=255)
                        cAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                                                                                                                             <  1:190689/150‑1 (MQ=255)
                        cAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                                                                                                                             <  1:3732138/150‑1 (MQ=255)
                        cAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                                                                                                                             <  2:4314703/150‑1 (MQ=255)
                           aCGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGc                                                                                                                          <  1:3152254/150‑1 (MQ=255)
                           aCGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGc                                                                                                                          <  1:3532767/150‑1 (MQ=255)
                            cGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCg                                                                                                                         >  1:2711998/1‑150 (MQ=255)
                             gTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCgg                                                                                                                        >  2:2558598/1‑150 (MQ=255)
                                tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                     <  2:344617/150‑1 (MQ=255)
                                   ataAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACt                                                                                                                  <  1:4514856/150‑1 (MQ=255)
                                     aaaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCATGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTc                                                                                                                <  1:474653/150‑1 (MQ=255)
                                     aaaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTc                                                                                                                >  1:2970714/1‑150 (MQ=255)
                                       aTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAt                                                                                                              >  1:537555/1‑150 (MQ=255)
                                             cGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaag                                                                                                        >  2:1781857/1‑150 (MQ=255)
                                               cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                                                                                      <  2:836547/150‑1 (MQ=255)
                                                 cTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                                                                                    <  2:3880488/150‑1 (MQ=255)
                                                   gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                                                                  >  1:2807272/1‑150 (MQ=255)
                                                     cGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACa                                                                                                <  2:1391465/150‑1 (MQ=255)
                                                      gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                                                               >  2:3028960/1‑150 (MQ=255)
                                                         aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                                                                            >  2:685569/1‑150 (MQ=255)
                                                           ggTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAttt                                                                                          <  2:3720852/150‑1 (MQ=255)
                                                            gTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTg                                                                                         >  1:1994975/1‑150 (MQ=255)
                                                             ttGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                                                                        <  1:3449178/150‑1 (MQ=255)
                                                             ttGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                                                                        >  1:3332342/1‑150 (MQ=255)
                                                             ttGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                                                                        <  1:2312632/150‑1 (MQ=255)
                                                               gTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                                                                                      <  2:3160259/150‑1 (MQ=255)
                                                                 ccGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                                                                                    <  2:3788883/150‑1 (MQ=255)
                                                                   ggCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTga                                                                                  <  2:2154909/150‑1 (MQ=255)
                                                                     cTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgaga                                                                                >  2:3198618/1‑150 (MQ=255)
                                                                       ccTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                                                                              <  1:2003260/150‑1 (MQ=255)
                                                                         tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            <  1:1332463/150‑1 (MQ=255)
                                                                         tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            <  1:3946120/150‑1 (MQ=255)
                                                                         tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            <  1:1686524/150‑1 (MQ=255)
                                                                         tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            <  1:3877681/150‑1 (MQ=255)
                                                                          gaagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATg                                                                           <  2:2636787/150‑1 (MQ=255)
                                                                           aagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                                                                          >  1:2421751/1‑150 (MQ=255)
                                                                           aagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                                                                          >  1:2416249/1‑150 (MQ=255)
                                                                            agaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTAACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCg                                                                         >  2:191538/1‑150 (MQ=255)
                                                                              aaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                                                                       >  1:34537/1‑150 (MQ=255)
                                                                                cTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc                                                                     <  1:4531445/150‑1 (MQ=255)
                                                                                cTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc                                                                     >  2:1602821/1‑150 (MQ=255)
                                                                                 ttCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCATGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCa                                                                    >  2:4401202/1‑150 (MQ=255)
                                                                                   ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAtt                                                                  >  1:2627737/1‑150 (MQ=255)
                                                                                     aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                                                                >  1:3570547/1‑150 (MQ=255)
                                                                                       gATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                                                              <  2:435394/150‑1 (MQ=255)
                                                                                       gATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                                                              <  2:435583/150‑1 (MQ=255)
                                                                                        aTTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGa                                                             >  1:799735/1‑150 (MQ=255)
                                                                                        aTTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGa                                                             >  1:2876489/1‑150 (MQ=255)
                                                                                            ttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg                                                         <  2:4073545/150‑1 (MQ=255)
                                                                                               aaCGGGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCa                                                      >  1:109162/1‑150 (MQ=255)
                                                                                               aaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCa                                                      >  1:3478859/1‑150 (MQ=255)
                                                                                                aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc                                                     <  1:3973605/150‑1 (MQ=255)
                                                                                                 cGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATGAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACg                                                    <  1:2744282/150‑1 (MQ=255)
                                                                                                  gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGc                                                   >  1:2305606/1‑150 (MQ=255)
                                                                                                  gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGc                                                   >  2:1847489/1‑150 (MQ=255)
                                                                                                   aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGcc                                                  <  2:459740/150‑1 (MQ=255)
                                                                                                        gCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATGGAATTGCCACGCCTAAAt                                             >  1:4182816/1‑150 (MQ=255)
                                                                                                           ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaa                                          <  2:1245443/150‑1 (MQ=255)
                                                                                                           ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaa                                          <  2:1611824/150‑1 (MQ=255)
                                                                                                            tACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaa                                         >  2:1366644/1‑150 (MQ=255)
                                                                                                             aCTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaaa                                        >  1:2133103/1‑150 (MQ=255)
                                                                                                             aCTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaaa                                        >  1:1294991/1‑150 (MQ=255)
                                                                                                                ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACca                                     >  1:347006/1‑150 (MQ=255)
                                                                                                                    ggCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAg                                 <  1:2279184/150‑1 (MQ=255)
                                                                                                                       cAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTcc                              <  1:438157/150‑1 (MQ=255)
                                                                                                                       cAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTcc                              <  1:4355792/150‑1 (MQ=255)
                                                                                                                       cAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTcc                              <  1:1730104/150‑1 (MQ=255)
                                                                                                                         gCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                            <  2:3049778/150‑1 (MQ=255)
                                                                                                                         gCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCATGTTACTAATCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATTCGTATCATTATAGAATTGCCAGGACTAAATGAAAACCACAAGTCCCt                            <  2:3656838/150‑1 (MQ=255)
                                                                                                                         gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATGGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                            >  2:1887222/1‑150 (MQ=255)
                                                                                                                           tttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                          <  2:2043372/150‑1 (MQ=255)
                                                                                                                            ttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                         >  1:4128676/1‑150 (MQ=255)
                                                                                                                                ttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAat                     >  2:791723/1‑150 (MQ=255)
                                                                                                                                  aGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                   >  2:1253662/1‑150 (MQ=255)
                                                                                                                                    ggCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACACGTCCCTGGCAAATATTg                 <  2:1476562/150‑1 (MQ=255)
                                                                                                                                    ggCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTg                 <  2:945833/150‑1 (MQ=255)
                                                                                                                                    ggCGGTATTGTTATCCGTCAAGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTg                 <  2:2454389/150‑1 (MQ=255)
                                                                                                                                    ggCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTg                 <  1:2369172/150‑1 (MQ=255)
                                                                                                                                       ggTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTGGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta              >  1:2099623/1‑150 (MQ=255)
                                                                                                                                        gTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa             <  2:1448625/150‑1 (MQ=255)
                                                                                                                                        gTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa             <  2:1634298/150‑1 (MQ=255)
                                                                                                                                         tATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaat            <  2:604041/150‑1 (MQ=255)
                                                                                                                                             gTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaaa        >  2:4097358/1‑150 (MQ=255)
                                                                                                                                               tATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAAc      >  2:2800551/1‑150 (MQ=255)
                                                                                                                                                aTCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACg     >  1:1944603/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGG  >  CP009273/705950‑706246

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: