Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 172,368 A→G M303V (ATG→GTG)  fhuB → iron(3+)‑hydroxamate import ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913172,3680AG99.1% 422.7 / ‑4.3 112M303V (ATG→GTG) fhuBiron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (50/61);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 118 total aligned reads displayed.

ATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCC  >  NC_000913/172221‑172515
                                                                                                                                                   |                                                                                                                                                   
aTCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTgtgggtgg                                                                                                                                                  <  2:306786/151‑1 (MQ=255)
 tCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGtgggtg                                                                                                                                                   >  2:1095908/1‑149 (MQ=255)
  cGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGtgggtgga                                                                                                                                                 >  1:1189889/1‑150 (MQ=255)
   ggCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGtgggtgga                                                                                                                                                 >  2:1493353/1‑149 (MQ=255)
    gCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAgt                                                                                                                                              <  1:92064/151‑1 (MQ=255)
      tttATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAgtgt                                                                                                                                            >  2:827696/1‑151 (MQ=255)
      tttATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAgtgt                                                                                                                                            >  2:263898/1‑151 (MQ=255)
       ttATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAgtgt                                                                                                                                            <  2:1158466/150‑1 (MQ=255)
       ttATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTc                                                                                                                                           >  2:1051928/1‑151 (MQ=255)
        tATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTc                                                                                                                                           <  1:1225165/150‑1 (MQ=255)
          tCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCAc                                                                                                                                        >  2:817888/1‑151 (MQ=255)
           cGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCAc                                                                                                                                        <  1:1181025/150‑1 (MQ=255)
            gggTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCg                                                                                                                                      >  1:83904/1‑151 (MQ=255)
                 gttCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTc                                                                                                                                  <  2:1603598/150‑1 (MQ=255)
                 gttCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCg                                                                                                                                 <  1:1442521/151‑1 (MQ=255)
                  ttCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCgg                                                                                                                                >  1:959190/1‑151 (MQ=255)
                    cgcgCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTc                                                                                                                              >  2:1070029/1‑151 (MQ=255)
                     gcgcCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCa                                                                                                                             <  2:1328719/151‑1 (MQ=255)
                       gccgcTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCAc                                                                                                                            >  1:1405062/1‑150 (MQ=255)
                        ccgcTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTg                                                                                                                          >  2:1524036/1‑151 (MQ=255)
                         cgcTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGc                                                                                                                         >  2:904858/1‑151 (MQ=255)
                             gctgGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGtt                                                                                                                      <  2:976537/150‑1 (MQ=255)
                               tgGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGAt                                                                                                                   >  2:736777/1‑151 (MQ=255)
                                 gCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATc                                                                                                                  >  2:1395017/1‑150 (MQ=255)
                                  cAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCgg                                                                                                                >  2:971232/1‑151 (MQ=255)
                                    aaaaTGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTg                                                                                                              <  2:859752/151‑1 (MQ=255)
                                      aaTGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTgcg                                                                                                            <  2:1073443/151‑1 (MQ=255)
                                       aTGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTgcgc                                                                                                           >  2:34499/1‑151 (MQ=255)
                                       aTGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTgcgc                                                                                                           <  2:1027000/151‑1 (MQ=255)
                                              gggCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCgctgct                                                                                                    >  1:185106/1‑151 (MQ=255)
                                                   cggcgTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACtgt                                                                                               <  2:1189979/151‑1 (MQ=255)
                                                   cggcgTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACtgt                                                                                               >  2:206335/1‑151 (MQ=255)
                                                    ggcgTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACtgtg                                                                                              >  2:1160331/1‑151 (MQ=255)
                                                     gcgTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTgg                                                                                             <  1:936483/151‑1 (MQ=255)
                                                      cgTCTGCTACCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGc                                                                                            <  2:759252/151‑1 (MQ=255)
                                                       gTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCt                                                                                           <  1:533017/151‑1 (MQ=255)
                                                           gctgcCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTg                                                                                       >  2:89286/1‑151 (MQ=255)
                                                           gctgcCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTg                                                                                       <  2:260320/151‑1 (MQ=255)
                                                            ctgcCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTg                                                                                       >  1:1577312/1‑150 (MQ=255)
                                                              gcCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCcg                                                                                    <  1:179961/151‑1 (MQ=255)
                                                               ccACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCcgc                                                                                   <  1:1551251/151‑1 (MQ=255)
                                                               ccACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCcgc                                                                                   <  2:794556/151‑1 (MQ=255)
                                                                cACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCcgcg                                                                                  <  2:957994/151‑1 (MQ=255)
                                                                   gACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGttt                                                                               <  1:856375/151‑1 (MQ=255)
                                                                    aCTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTa                                                                              <  1:263875/151‑1 (MQ=255)
                                                                    aCTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTa                                                                              <  2:1440575/151‑1 (MQ=255)
                                                                    aCTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTa                                                                              <  2:744522/151‑1 (MQ=255)
                                                                       gATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACg                                                                            <  1:1129939/150‑1 (MQ=255)
                                                                        aTGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACgca                                                                          <  2:308921/151‑1 (MQ=255)
                                                                             ggCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCAtt                                                                     <  2:1542358/151‑1 (MQ=255)
                                                                              gCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTa                                                                    <  1:1245772/151‑1 (MQ=255)
                                                                                gtcgtTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAgc                                                                  <  1:964399/151‑1 (MQ=255)
                                                                                gtcgtTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAgc                                                                  <  1:1394925/151‑1 (MQ=255)
                                                                                  cgtTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAgcgc                                                                >  1:777944/1‑151 (MQ=255)
                                                                                     tGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAgcgcgc                                                              <  1:1006317/150‑1 (MQ=255)
                                                                                      gATTGGCGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTCCCGCGTTTACGCAGCATTAgcgc                                                                >  1:1214373/1‑147 (MQ=255)
                                                                                          ggTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGata                                                        <  2:1714344/151‑1 (MQ=255)
                                                                                          ggTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGat                                                         <  2:471149/150‑1 (MQ=255)
                                                                                           gTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGatat                                                       <  2:1504310/151‑1 (MQ=255)
                                                                                             gcgcTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGa                                                     <  2:1176405/151‑1 (MQ=255)
                                                                                               gcTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAg                                                   <  1:1493252/151‑1 (MQ=255)
                                                                                               gcTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGa                                                     >  2:1603019/1‑149 (MQ=255)
                                                                                                cTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAgg                                                  >  1:849426/1‑151 (MQ=255)
                                                                                                 tGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGt                                                 >  2:407947/1‑151 (MQ=255)
                                                                                                   aTCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAgg                                                  <  1:693536/148‑1 (MQ=255)
                                                                                                    tCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGCTTACGCAGCATTAGCGCGCCGGATATGAAGGTCa                                               >  1:1091542/1‑150 (MQ=255)
                                                                                                        ctGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGAt                                          >  2:496940/1‑151 (MQ=255)
                                                                                                         tGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATc                                         >  2:295498/1‑151 (MQ=255)
                                                                                                          ggCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGa                                           <  2:1459817/148‑1 (MQ=255)
                                                                                                          ggCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGAt                                          <  2:186591/149‑1 (MQ=255)
                                                                                                          ggCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCg                                        <  1:937542/151‑1 (MQ=255)
                                                                                                            cTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCgtg                                      >  1:567985/1‑151 (MQ=255)
                                                                                                              ttCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCgtgt                                     <  1:817823/150‑1 (MQ=255)
                                                                                                                ccGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTcgc                                  >  1:53327/1‑151 (MQ=255)
                                                                                                                ccGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTcg                                   <  1:27615/150‑1 (MQ=255)
                                                                                                                 cGATCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTcgcg                                 <  1:1217000/151‑1 (MQ=255)
                                                                                                                   aTCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGc                               >  1:800926/1‑151 (MQ=255)
                                                                                                                   aTCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGc                               >  1:498577/1‑151 (MQ=255)
                                                                                                                    tCAAATCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCgg                                >  1:1705036/1‑149 (MQ=255)
                                                                                                                       aaTCATCCTCTGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGaa                           <  1:1323308/151‑1 (MQ=255)
                                                                                                                           ctcCTCTGGCTGACTCGAGTGTGGGTGGAAGTGTCGACCGGTTCGCTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCTTTAGCGCGCCGGCGATGCAGGGCAACGATCGTGTCTCGGCCGAACGcc                       >  1:411367/2‑151 (MQ=255)
                                                                                                                              ctctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCaa                     >  1:1104371/1‑150 (MQ=255)
                                                                                                                              ctctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCaa                     >  1:1465243/1‑150 (MQ=255)
                                                                                                                               tctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAAc                    >  1:1690689/1‑150 (MQ=255)
                                                                                                                                ctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAAc                    <  2:1301505/149‑1 (MQ=255)
                                                                                                                                ctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAAc                    <  2:602339/149‑1 (MQ=255)
                                                                                                                                ctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAAc                    <  1:150025/149‑1 (MQ=255)
                                                                                                                                ctGGCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAt                  >  2:1235384/1‑151 (MQ=255)
                                                                                                                                  ggCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtg                 >  1:97528/1‑150 (MQ=255)
                                                                                                                                   gCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtgtg               <  2:83908/151‑1 (MQ=255)
                                                                                                                                   gCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtgtg               <  2:1148161/151‑1 (MQ=255)
                                                                                                                                   gCTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtgtg               <  1:1160244/151‑1 (MQ=255)
                                                                                                                                    cTGACTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc              >  2:348543/1‑151 (MQ=255)
                                                                                                                                        cTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTgg           >  2:61209/1‑150 (MQ=255)
                                                                                                                                        cTCGCGTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTgg           >  1:1047220/1‑150 (MQ=255)
                                                                                                                                          cgcgTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCg         <  2:651733/150‑1 (MQ=255)
                                                                                                                                           gcgTGTGGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGtt       <  2:1315466/151‑1 (MQ=255)
                                                                                                                                              tgtgGGTGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGc    >  2:817686/1‑151 (MQ=255)
                                                                                                                                               gtgggtggAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGcc   >  2:542753/1‑151 (MQ=255)
                                                                                                                                                tgggtggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGccc  >  2:1425413/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCC  >  NC_000913/172221‑172515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: