Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 706,097 A→C 91.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,0970AC91.4% 874.2 / 36.5 233T2T (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 (12/8);  new base C (93/120);  total (105/128)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 237 total aligned reads displayed.

AGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGG  >  CP009273/705948‑706246
                                                                                                                                                     |                                                                                                                                                     
   ccAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTc                                                                                                                                                    <  1:1597471/150‑1 (MQ=255)
       tttCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGc                                                                                                                                                <  1:3651637/150‑1 (MQ=255)
           ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGaa                                                                                                                                            <  1:2999874/150‑1 (MQ=255)
           ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGaa                                                                                                                                            >  2:321040/1‑150 (MQ=255)
            tcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAAt                                                                                                                                           <  2:4563374/150‑1 (MQ=255)
             cACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATc                                                                                                                                          >  1:4172316/1‑150 (MQ=255)
             cACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATc                                                                                                                                          >  2:2067038/1‑150 (MQ=255)
               cccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTg                                                                                                                                        <  2:3132777/150‑1 (MQ=255)
                ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGt                                                                                                                                       <  1:1095842/150‑1 (MQ=255)
                ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGt                                                                                                                                       <  1:980667/150‑1 (MQ=255)
                 cATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTc                                                                                                                                      <  2:4250585/150‑1 (MQ=255)
                  atatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTcc                                                                                                                                     <  1:284303/150‑1 (MQ=255)
                    ctCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTg                                                                                                                                   <  2:4451647/149‑1 (MQ=255)
                     tCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGt                                                                                                                                  >  2:4071116/1‑150 (MQ=255)
                        aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTAc                                                                                                                               >  2:3040233/1‑150 (MQ=255)
                         tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACt                                                                                                                              <  1:765726/150‑1 (MQ=255)
                         tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACt                                                                                                                              >  1:1778777/1‑150 (MQ=255)
                         tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACt                                                                                                                              >  1:1452833/1‑150 (MQ=255)
                         tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACt                                                                                                                              >  2:2165600/1‑150 (MQ=255)
                          cAGACGTTTGTATAAATCTTCCGCACTGACATGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTa                                                                                                                             <  1:3279117/150‑1 (MQ=255)
                            gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAg                                                                                                                           <  2:100744/150‑1 (MQ=255)
                             aCGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGc                                                                                                                          <  1:4349929/150‑1 (MQ=255)
                              cGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCg                                                                                                                         >  1:3153312/1‑150 (MQ=255)
                              cGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCg                                                                                                                         <  2:1963463/150‑1 (MQ=255)
                                 ttGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGtt                                                                                                                      <  1:1523256/150‑1 (MQ=255)
                                  tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                     <  1:4204254/150‑1 (MQ=255)
                                  tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                     >  2:618252/1‑150 (MQ=255)
                                  tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                     >  2:823437/1‑150 (MQ=255)
                                  tGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                                                                                     <  1:1156970/150‑1 (MQ=255)
                                    tataAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTGTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCAc                                                                                                                   >  2:1294210/1‑150 (MQ=255)
                                         aTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAt                                                                                                              <  2:2779142/150‑1 (MQ=255)
                                          tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                                                                             >  2:1930089/1‑150 (MQ=255)
                                          tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                                                                             <  1:639014/150‑1 (MQ=255)
                                           cTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTa                                                                                                            <  2:2266710/150‑1 (MQ=255)
                                            ttCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTag                                                                                                           <  2:3796510/150‑1 (MQ=255)
                                             tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTaga                                                                                                          >  1:3556952/1‑150 (MQ=255)
                                                gCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaaga                                                                                                       <  2:1806384/150‑1 (MQ=255)
                                                  aCTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGaagaag                                                                                                     >  1:4053331/1‑150 (MQ=255)
                                                     gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                                                                  >  2:4344728/1‑150 (MQ=255)
                                                      aCGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGAc                                                                                                 <  1:1603777/150‑1 (MQ=255)
                                                       cGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACa                                                                                                >  2:1820805/1‑150 (MQ=255)
                                                        gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                                                               <  2:134356/150‑1 (MQ=255)
                                                         tGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                                                              >  1:845584/1‑150 (MQ=255)
                                                            tGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTAGTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                                                                           <  2:4232754/150‑1 (MQ=255)
                                                                 gTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                                                                                      >  1:3090283/1‑150 (MQ=255)
                                                                   ccGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                                                                                    >  1:4057806/1‑150 (MQ=255)
                                                                      gCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                                                                 <  2:3999704/150‑1 (MQ=255)
                                                                      gCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                                                                 <  2:1160821/150‑1 (MQ=255)
                                                                        tCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGAt                                                                               <  2:2560049/150‑1 (MQ=255)
                                                                          cTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                                                             <  2:3697027/150‑1 (MQ=255)
                                                                          cTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                                                             <  1:2448550/150‑1 (MQ=255)
                                                                           tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            <  2:120230/150‑1 (MQ=255)
                                                                           tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAAt                                                                            >  1:491261/1‑150 (MQ=255)
                                                                              agaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCg                                                                         >  1:3845904/1‑150 (MQ=255)
                                                                               gaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGt                                                                        <  2:1577805/150‑1 (MQ=255)
                                                                               gaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGt                                                                        >  2:66234/1‑150 (MQ=255)
                                                                                 aCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTAt                                                                      >  1:683448/1‑150 (MQ=255)
                                                                                   ttCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCa                                                                    >  1:3169159/1‑150 (MQ=255)
                                                                                    tCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAt                                                                   <  2:2308476/150‑1 (MQ=255)
                                                                                     ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAtt                                                                  <  2:3713077/150‑1 (MQ=255)
                                                                                     ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAtt                                                                  >  1:442703/1‑150 (MQ=255)
                                                                                       aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                                                                >  1:373481/1‑150 (MQ=255)
                                                                                       aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                                                                >  1:3945934/1‑150 (MQ=255)
                                                                                       aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                                                                <  1:3444902/150‑1 (MQ=255)
                                                                                        ggATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata                                                               <  1:2687760/150‑1 (MQ=255)
                                                                                        ggATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata                                                               <  2:418953/150‑1 (MQ=255)
                                                                                        ggATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata                                                               <  2:786952/150‑1 (MQ=255)
                                                                                              ttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg                                                         <  1:2746642/150‑1 (MQ=255)
                                                                                                aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGcc                                                       >  2:3412269/1‑150 (MQ=255)
                                                                                                aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGcc                                                       <  1:576287/150‑1 (MQ=255)
                                                                                                    gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGc                                                   >  1:1236994/1‑150 (MQ=255)
                                                                                                       gAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTa                                                <  1:2399522/150‑1 (MQ=255)
                                                                                                        aaGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa                                               <  1:2604530/150‑1 (MQ=255)
                                                                                                           cGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATg                                            >  2:2367607/1‑150 (MQ=255)
                                                                                                            gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTACACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                           <  1:1895895/150‑1 (MQ=255)
                                                                                                                  ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTTTCATTATAGAATTGCCACGCCTAAATGAAAACca                                     <  2:4623797/150‑1 (MQ=255)
                                                                                                                  ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACca                                     >  2:1297609/1‑150 (MQ=255)
                                                                                                                   tCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcac                                    <  2:3258826/150‑1 (MQ=255)
                                                                                                                    cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcaca                                   <  2:3547748/150‑1 (MQ=255)
                                                                                                                      ggCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAg                                 <  1:3835691/150‑1 (MQ=255)
                                                                                                                       gCCAGCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGt                                >  1:1078215/1‑150 (MQ=255)
                                                                                                                           gCTTTCTTTAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATTCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                            <  2:4304647/150‑1 (MQ=255)
                                                                                                                            ctttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTg                           <  2:148059/150‑1 (MQ=255)
                                                                                                                             tttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                          <  1:936271/150‑1 (MQ=255)
                                                                                                                              ttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                         >  1:2380745/1‑150 (MQ=255)
                                                                                                                              ttctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                         <  2:3344556/150‑1 (MQ=255)
                                                                                                                               tctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCa                        <  1:1993407/150‑1 (MQ=255)
                                                                                                                               tctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCa                        <  2:50423/150‑1 (MQ=255)
                                                                                                                                ctttAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCaa                       <  1:960178/150‑1 (MQ=255)
                                                                                                                                 tttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCaaa                      >  2:3839960/1‑150 (MQ=255)
                                                                                                                                   tAGGGCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAata                    >  2:169867/1‑150 (MQ=255)
                                                                                                                                     gggCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                  <  2:275555/150‑1 (MQ=255)
                                                                                                                                     gggCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                  >  2:2331932/1‑150 (MQ=255)
                                                                                                                                      ggCGGTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTg                 >  1:2194488/1‑150 (MQ=255)
                                                                                                                                         ggTATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta              >  2:693620/1‑150 (MQ=255)
                                                                                                                                           tATTGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaat            <  1:1366464/150‑1 (MQ=255)
                                                                                                                                              tGTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaa         <  2:1204580/150‑1 (MQ=255)
                                                                                                                                               gTTATCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaaa        >  1:1714356/1‑150 (MQ=255)
                                                                                                                                                  aTCCGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGACACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAAct     >  2:3670756/1‑149 (MQ=255)
                                                                                                                                                     aGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTgg  <  2:3733231/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGG  >  CP009273/705948‑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: