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.

CACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  NC_000913/4514186‑4514484
                                                                                                                                                      |                                                                                                                                                     
cacaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCa                                                                                                                                                       >  2:1468005/1‑150 (MQ=255)
  caATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAAACACATACCAGCTAAAAGCCCGGCAAg                                                                                                                                                     >  2:152623/1‑150 (MQ=255)
   aaTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc                                                                                                                                                    >  1:1224772/1‑150 (MQ=255)
    aTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACTGATATATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGcc                                                                                                                                                   <  1:1682821/150‑1 (MQ=255)
       cGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCggg                                                                                                                                                <  2:221297/150‑1 (MQ=255)
          tgtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGctt                                                                                                                                             <  1:413738/150‑3 (MQ=255)
           gtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt                                                                                                                                            <  2:358440/150‑1 (MQ=255)
           gtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt                                                                                                                                            >  2:1022377/1‑150 (MQ=255)
            tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGATCCGTGGCGGCAAAGGCGGGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAACACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCtttt                                                                                                                                           <  1:348053/150‑2 (MQ=255)
              gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttag                                                                                                                                         <  2:804716/150‑4 (MQ=255)
              gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAac                                                                                                                                         <  1:340850/150‑1 (MQ=255)
                gTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                       <  1:1643502/150‑1 (MQ=255)
                   ttttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATc                                                                                                                                    >  2:848455/1‑150 (MQ=255)
                    tttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCa                                                                                                                                   <  1:1290472/150‑1 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  1:728997/1‑150 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGAGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               >  2:142617/1‑150 (MQ=255)
                        gAGTGTAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGATGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAAc                                                                                                                               <  2:841679/150‑1 (MQ=255)
                          gtgtAAACGTGCCGTGTTCGTCCTGCACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGATAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACAGCTCAGAACtt                                                                                                                             >  1:847939/1‑150 (MQ=255)
                          gtgtAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACtt                                                                                                                             <  1:701164/150‑1 (MQ=255)
                                  gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACTTCAGAACTTCAACGAcc                                                                                                                     <  1:608358/150‑1 (MQ=255)
                                  gTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAcc                                                                                                                     <  2:1393343/150‑1 (MQ=255)
                                     ccGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCt                                                                                                                  <  2:77909/150‑1 (MQ=255)
                                     ccGTGTTCGTCATGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCt                                                                                                                  >  2:407915/1‑150 (MQ=255)
                                        tgtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                               <  1:803547/150‑1 (MQ=255)
                                          ttCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCata                                                                                                             <  2:107404/150‑1 (MQ=255)
                                            cGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCCGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGCCGTTAACACATCAGAACTTCAACGACCCCTCCAgata                                                                                                           >  2:619682/1‑150 (MQ=255)
                                            cGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCAtata                                                                                                           >  2:423728/1‑150 (MQ=255)
                                               ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAACAGCCCGGCAAGCCGGGCGTTACCACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                        >  1:1364987/1‑150 (MQ=255)
                                               ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAg                                                                                                        >  2:1525418/1‑150 (MQ=255)
                                               ccTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGAAGCAGTCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCAGGCATGCCGGGCGTTAACACATCAGAACTTCAACGAGCACTGCATATACAg                                                                                                        <  2:301821/150‑1 (MQ=255)
                                                       gTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGcgcg                                                                                                <  1:1194544/150‑1 (MQ=255)
                                                         ggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGc                                                                                              <  1:1137428/150‑1 (MQ=255)
                                                             gcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGc                                                                                          >  2:1008730/1‑150 (MQ=255)
                                                              cAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGcc                                                                                         <  1:784212/150‑1 (MQ=255)
                                                                aaGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTg                                                                                       <  2:375286/150‑1 (MQ=255)
                                                                 aGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCCGCACTGCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTg                                                                                       >  2:732502/1‑150 (MQ=255)
                                                                 aGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGc                                                                                      <  1:1067899/150‑1 (MQ=255)
                                                                  ggCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCa                                                                                     <  2:664536/150‑1 (MQ=255)
                                                                    ggTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATa                                                                                   <  2:1728165/149‑1 (MQ=255)
                                                                      tGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGTGTGCGCGGCTGGCCTGCATAGa                                                                                 >  2:1682660/1‑150 (MQ=255)
                                                                        gCTGATCACCAGCAGCAGGACTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATg                                                                               <  2:919290/150‑1 (MQ=255)
                                                                        gCTGATCACCAGCAGCAGGACTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACAATCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATg                                                                               <  1:19622/150‑1 (MQ=255)
                                                                            aTCACCAGCAGCAGTCCTGC‑AACAAGAAATCGGATCAATGCCAACATCATCACATTCCAGCTAAAAGCCCTGCACGCCGTGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCtt                                                                           <  1:292788/150‑1 (MQ=255)
                                                                              cACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttg                                                                         >  2:939647/1‑150 (MQ=255)
                                                                                 cagcagcagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttg                                                                      >  1:1600753/1‑150 (MQ=255)
                                                                                 cagcagcagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttg                                                                      <  1:1341314/150‑1 (MQ=255)
                                                                                  agcagcagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgt                                                                     >  2:1488629/1‑150 (MQ=255)
                                                                                    cagcagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  1:198147/150‑1 (MQ=255)
                                                                                    cagcagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgttg                                                                   <  2:365256/150‑1 (MQ=255)
                                                                                     agcagGCCGGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTtgttgttgt                                                                  <  2:293038/150‑1 (MQ=255)
                                                                                       cagGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcg                                                                <  2:124837/150‑1 (MQ=255)
                                                                                         ggCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcgtc                                                              >  2:154589/1‑150 (MQ=255)
                                                                                          gCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCa                                                             >  1:1704365/1‑150 (MQ=255)
                                                                                            cTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATa                                                           <  1:939647/150‑1 (MQ=255)
                                                                                             tGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATaa                                                          >  1:790712/1‑150 (MQ=255)
                                                                                              gC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGAGGCCTTTGTTGTTGTCGTCATAag                                                         >  2:253476/1‑150 (MQ=255)
                                                                                               c‑AAAAAGAAAACGGATAAATGCCAACATAATCAGATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAagc                                                        >  2:729516/1‑149 (MQ=255)
                                                                                               c‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGAAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAaga                                                        >  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)
                                                                                                                                                      |                                                                                                                                                     
CACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGC‑AAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAG  >  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: ‑