Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,514,336 A→T 5.4% intergenic (‑28/+17) 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,3360AT 5.4% 319.3 / 4.6 113intergenic (‑28/+17)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/58);  new base T (2/5);  total (49/64)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 114 total aligned reads displayed.

CAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  NC_000913/4514188‑4514484
                                                                                                                                                     |                                                                                                                                                    
caATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAAACACATACCAGCTAAAAGCCCGGCAAg                                                                                                                                                     >  2:152623/1‑150 (MQ=255)
 aaTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc                                                                                                                                                    >  1:1224772/1‑150 (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                                                                                                                                         <  1:1611409/150‑4 (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)
                      gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  1:728997/1‑150 (MQ=255)
                      gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGAGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  2:142617/1‑150 (MQ=255)
                      gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGATGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               <  2:841679/150‑1 (MQ=255)
                       aGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACt                                                                                                                              >  2:613212/1‑150 (MQ=255)
                        gtgtAAACGTGCCGTGTTCGTCCTGCACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGATAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACAGCTCAGAACtt                                                                                                                             >  1:847939/1‑150 (MQ=255)
                        gtgtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACtt                                                                                                                             <  1:701164/150‑1 (MQ=255)
                            aaaCGTGCCGTGTTCGTCCTGATCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAc                                                                                                                         >  2:197455/1‑150 (MQ=255)
                                gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACTTCAGAACTTCAACGAcc                                                                                                                     <  1:608358/150‑1 (MQ=255)
                                gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcc                                                                                                                     <  1:277537/150‑1 (MQ=255)
                                gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcc                                                                                                                     <  2:1393343/150‑1 (MQ=255)
                                  gCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcccc                                                                                                                   <  1:848455/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)
                                             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:521494/150‑1 (MQ=255)
                                                       ggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGc                                                                                              <  1:1137428/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)
                                                                      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)
                                                                              ccAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgtt                                                                       >  1:1229987/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                                                                   <  2:365256/150‑1 (MQ=255)
                                                                                  cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  2:794215/150‑1 (MQ=255)
                                                                                  cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  1:198147/150‑1 (MQ=255)
                                                                                   agcagGCCGGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTtgttgttgt                                                                  <  2:293038/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)
                                                                                            gCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCTGCTAGCCGGGCGTTAACACATGAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAag                                                         <  1:727819/150‑1 (MQ=255)
                                                                                            gCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAag                                                         <  1:1468005/150‑1 (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)
                                                                                                               aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCt                                       >  1:187065/1‑150 (MQ=255)
                                                                                                               aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCAGATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTTGTCCt                                       >  2:118998/1‑150 (MQ=255)
                                                                                                                  gCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCg                                                     >  2:1014931/1‑133 (MQ=255)
                                                                                                                  gCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCg                                                     <  1:1014931/133‑1 (MQ=255)
                                                                                                                   ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTc                                   <  1:1345781/150‑1 (MQ=255)
                                                                                                                   ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTc                                   <  1:197455/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)
                                                                                                                             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)
                                                                                                                                aTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                      >  2:1289465/1‑150 (MQ=255)
                                                                                                                                  tCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGGGGCTGGCCTGCATAGATGCCTTTGTTTTTGTCGTCATAAGAGCGGATGAATTAGTCCTGGTCGAAGATGTTTTTcac                    <  1:417415/150‑1 (MQ=255)
                                                                                                                                   ccAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTGTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTcccc                   >  2:1515020/1‑147 (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)
                                                                                                                                         aaaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCAGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAt             >  2:323989/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)
                                                                                                                                                   cAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa   <  1:407915/150‑1 (MQ=255)
                                                                                                                                                    aaGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAg  <  1:142617/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  NC_000913/4514188‑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: