Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,514,335 A→T 5.7% intergenic (‑27/+18) fecB ← / ← fecA ferric citrate ABC transporter periplasmic binding protein/TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,3350AT 5.7% 333.2 / 5.1 113intergenic (‑27/+18)fecB/fecAferric citrate ABC transporter periplasmic binding protein/TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Reads supporting (aligned to +/- strand):  ref base A (47/59);  new base T (2/5);  total (49/64)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

CACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  NC_000913/4514186‑4514484
                                                                                                                                                      |                                                                                                                                                     
cacaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCa                                                                                                                                                       >  2:1468005/1‑150 (MQ=255)
  caATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAAACACATACCAGCTAAAAGCCCGGCAAg                                                                                                                                                     >  2:152623/1‑150 (MQ=255)
   aaTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc                                                                                                                                                    >  1:1224772/1‑150 (MQ=255)
    aTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACTGATATATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGcc                                                                                                                                                   <  1:1682821/150‑1 (MQ=255)
       cGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCggg                                                                                                                                                <  2:221297/150‑1 (MQ=255)
          tgtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGctt                                                                                                                                             <  1:413738/150‑3 (MQ=255)
           gtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt                                                                                                                                            <  2:358440/150‑1 (MQ=255)
           gtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt                                                                                                                                            >  2:1022377/1‑150 (MQ=255)
            tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGATCCGTGGCGGCAAAGGCGGGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAACACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCtttt                                                                                                                                           <  1:348053/150‑2 (MQ=255)
              gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttag                                                                                                                                         <  2:804716/150‑4 (MQ=255)
              gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAac                                                                                                                                         <  1:340850/150‑1 (MQ=255)
                gTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                       <  1:1643502/150‑1 (MQ=255)
                   ttttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATc                                                                                                                                    >  2:848455/1‑150 (MQ=255)
                    tttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCa                                                                                                                                   <  1:1290472/150‑1 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  1:728997/1‑150 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGAGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  2:142617/1‑150 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGATGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               <  2:841679/150‑1 (MQ=255)
                          gtgtAAACGTGCCGTGTTCGTCCTGCACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGATAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACAGCTCAGAACtt                                                                                                                             >  1:847939/1‑150 (MQ=255)
                          gtgtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACtt                                                                                                                             <  1:701164/150‑1 (MQ=255)
                                  gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACTTCAGAACTTCAACGAcc                                                                                                                     <  1:608358/150‑1 (MQ=255)
                                  gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcc                                                                                                                     <  2:1393343/150‑1 (MQ=255)
                                     ccGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCt                                                                                                                  <  2:77909/150‑1 (MQ=255)
                                     ccGTGTTCGTCATGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCt                                                                                                                  >  2:407915/1‑150 (MQ=255)
                                        tgtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                               <  1:803547/150‑1 (MQ=255)
                                          ttCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCata                                                                                                             <  2:107404/150‑1 (MQ=255)
                                            cGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCCGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGCCGTTAACACATCAGAACTTCAACGACCCCTCCAgata                                                                                                           >  2:619682/1‑150 (MQ=255)
                                            cGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCAtata                                                                                                           >  2:423728/1‑150 (MQ=255)
                                               ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAACAGCCCGGCAAGCCGGGCGTTACCACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                        >  1:1364987/1‑150 (MQ=255)
                                               ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                        >  2:1525418/1‑150 (MQ=255)
                                               ccTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGAAGCAGTCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCAGGCATGCCGGGCGTTAACACATCAGAACTTCAACGAGCACTGCATATACAg                                                                                                        <  2:301821/150‑1 (MQ=255)
                                                       gTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGcgcg                                                                                                <  1:1194544/150‑1 (MQ=255)
                                                         ggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGc                                                                                              <  1:1137428/150‑1 (MQ=255)
                                                             gcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGc                                                                                          >  2:1008730/1‑150 (MQ=255)
                                                              cAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGcc                                                                                         <  1:784212/150‑1 (MQ=255)
                                                                aaGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTg                                                                                       <  2:375286/150‑1 (MQ=255)
                                                                 aGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCCGCACTGCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTg                                                                                       >  2:732502/1‑150 (MQ=255)
                                                                 aGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGc                                                                                      <  1:1067899/150‑1 (MQ=255)
                                                                  ggCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCa                                                                                     <  2:664536/150‑1 (MQ=255)
                                                                    ggTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATa                                                                                   <  2:1728165/149‑1 (MQ=255)
                                                                      tGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGTGTGCGCGGCTGGCCTGCATAGa                                                                                 >  2:1682660/1‑150 (MQ=255)
                                                                        gCTGATCACCAGCAGCAGGACTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATg                                                                               <  2:919290/150‑1 (MQ=255)
                                                                        gCTGATCACCAGCAGCAGGACTGCAAAAAGAAAACGGATAAATGCCAACATAATCACAATCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATg                                                                               <  1:19622/150‑1 (MQ=255)
                                                                            aTCACCAGCAGCAGTCCTGCAACAAGAAATCGGATCAATGCCAACATCATCACATTCCAGCTAAAAGCCCTGCACGCCGTGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCtt                                                                           <  1:292788/150‑1 (MQ=255)
                                                                              cACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttg                                                                         >  2:939647/1‑150 (MQ=255)
                                                                                 cagcagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttg                                                                      >  1:1600753/1‑150 (MQ=255)
                                                                                 cagcagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttg                                                                      <  1:1341314/150‑1 (MQ=255)
                                                                                  agcagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgt                                                                     >  2:1488629/1‑150 (MQ=255)
                                                                                    cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  1:198147/150‑1 (MQ=255)
                                                                                    cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  2:365256/150‑1 (MQ=255)
                                                                                     agcagGCCGGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTtgttgttgt                                                                  <  2:293038/150‑1 (MQ=255)
                                                                                       cagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcg                                                                <  2:124837/150‑1 (MQ=255)
                                                                                         ggCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcgtc                                                              >  2:154589/1‑150 (MQ=255)
                                                                                          gCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCa                                                             >  1:1704365/1‑150 (MQ=255)
                                                                                            cTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATa                                                           <  1:939647/150‑1 (MQ=255)
                                                                                             tGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATaa                                                          >  1:790712/1‑150 (MQ=255)
                                                                                              gCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGAGGCCTTTGTTGTTGTCGTCATAag                                                         >  2:253476/1‑150 (MQ=255)
                                                                                               cAAAAAGAAAACGGATAAATGCCAACATAATCAGATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAagc                                                        >  2:729516/1‑149 (MQ=255)
                                                                                               cAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGAAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAaga                                                        >  1:1300833/1‑150 (MQ=255)
                                                                                                   aaaGAAAATGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCg                                                     <  1:1583366/150‑1 (MQ=255)
                                                                                                         aaCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGaa                                               >  2:1594778/1‑150 (MQ=255)
                                                                                                          aCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAag                                              <  2:743302/150‑1 (MQ=255)
                                                                                                           cGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagt                                             <  2:1262297/150‑1 (MQ=255)
                                                                                                           cGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCGGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagt                                             <  2:85081/150‑1 (MQ=255)
                                                                                                             gATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAagtag                                           <  1:1004818/150‑1 (MQ=255)
                                                                                                                 aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAATAGCGGATGAAGTAGTCCt                                       >  1:296962/1‑150 (MQ=255)
                                                                                                                 aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCt                                       >  1:187065/1‑150 (MQ=255)
                                                                                                                   tGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTgg                                     <  2:1470290/150‑1 (MQ=255)
                                                                                                                    gCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCg                                                     >  2:1014931/1‑133 (MQ=255)
                                                                                                                    gCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCg                                                     <  1:1014931/133‑1 (MQ=255)
                                                                                                                     ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTc                                   <  1:197455/150‑1 (MQ=255)
                                                                                                                     ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTc                                   <  1:1345781/150‑1 (MQ=255)
                                                                                                                     ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATACATGCCTTTGTTGTTGTCGACATAAGAGCGGATGAAGTAGTCCTGGTc                                   >  1:868253/1‑150 (MQ=255)
                                                                                                                       aaCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGa                                 <  2:349426/150‑1 (MQ=255)
                                                                                                                        aCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGaa                                >  2:1256359/1‑150 (MQ=255)
                                                                                                                         cATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAg                               <  2:143501/150‑1 (MQ=255)
                                                                                                                         cATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAg                               <  2:1600753/150‑1 (MQ=255)
                                                                                                                          aTAATCACATTCCAGCTAAATGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCGTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGa                              <  1:152623/150‑1 (MQ=255)
                                                                                                                          aTAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGa                              >  1:579875/1‑150 (MQ=255)
                                                                                                                               cacaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttt                         <  2:233114/150‑1 (MQ=255)
                                                                                                                               cacaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttt                         <  1:569259/150‑1 (MQ=255)
                                                                                                                                 caTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttttt                       <  1:1517514/150‑1 (MQ=255)
                                                                                                                                  aTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                      <  1:1525418/150‑1 (MQ=255)
                                                                                                                                    tCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGGGGCTGGCCTGCATAGATGCCTTTGTTTTTGTCGTCATAAGAGCGGATGAATTAGTCCTGGTCGAAGATGTTTTTcac                    <  1:417415/150‑1 (MQ=255)
                                                                                                                                       aGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTGTGTTGTTGTCGTCATAAGCGCGGATGATGTAGTCCTGGTCGAAGATGTATTTCAcccc                 >  2:303615/1‑150 (MQ=255)
                                                                                                                                        gCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCg                <  1:582799/150‑1 (MQ=255)
                                                                                                                                         gttaacgCCCGGCTTGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGATTTTCACACCGa               >  2:1362341/7‑150 (MQ=255)
                                                                                                                                           aaaaGCCCGGCAAGCCGGGCGTTACCACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAt             >  2:1200692/1‑150 (MQ=255)
                                                                                                                                           aaaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAt             >  2:404492/1‑150 (MQ=255)
                                                                                                                                           aaaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAt             >  2:1218692/1‑150 (MQ=255)
                                                                                                                                             acgCCCGGCTTGCCGGGCGTTAACACATCAGATCTTCAACGACCCCTGCAGATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTGGTTGTTGGCGTCATAATAGCGGATGAAGTAGACCTGGTCGAAGATGTTTGTCACACCGAATGc           >  1:752945/3‑150 (MQ=255)
                                                                                                                                             aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                      <  2:1095080/139‑1 (MQ=255)
                                                                                                                                             aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                      >  1:1095080/1‑139 (MQ=255)
                                                                                                                                             aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGc           >  2:770512/1‑150 (MQ=255)
                                                                                                                                                cccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAg        <  2:1224772/150‑1 (MQ=255)
                                                                                                                                                 ccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAgg       >  1:72120/1‑150 (MQ=255)
                                                                                                                                                 ccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAgg       >  1:1368717/1‑150 (MQ=255)
                                                                                                                                                  cGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGCGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGt      <  1:80104/150‑1 (MQ=255)
                                                                                                                                                   ggCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGtt     >  1:1482677/1‑150 (MQ=255)
                                                                                                                                                     cAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa   <  1:407915/150‑1 (MQ=255)
                                                                                                                                                      aaGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAg  <  1:142617/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  NC_000913/4514186‑4514484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: