Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 610,711 G→A 6.3% P595L (CCG→CTG)  fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913610,7110GA 6.3% 198.5 / 3.2 63P595L (CCG→CTG) fepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

TCGACAAATCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGTGAC  >  NC_000913/610568‑610856
                                                                                                                                               |                                                                                                                                                 
tCGACAAATCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACg                                                                                                                                              >  1:412397/1‑149 (MQ=255)
 cGACAAATCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGt                                                                                                                                             <  2:272345/149‑1 (MQ=255)
  gACAAATCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGtt                                                                                                                                            <  2:285356/149‑1 (MQ=255)
      aaTCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAAc                                                                                                                                        <  2:104392/149‑1 (MQ=255)
        ccTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGa                                                                                                                                      <  4:364848/148‑1 (MQ=255)
           tCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATcc                                                                                                                                   >  4:414754/1‑149 (MQ=255)
                 gcctgccAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAg                                                                                                                             >  4:292918/1‑149 (MQ=255)
                  cctgccAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGa                                                                                                                            >  4:341350/1‑149 (MQ=255)
                         gCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCaa                                                                                                                     >  1:529837/1‑149 (MQ=255)
                          cTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAac                                                                                                                    >  2:513677/1‑149 (MQ=255)
                                 ttGAGTTCGACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCt                                                                                                             >  3:523109/1‑149 (MQ=255)
                                  tGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCtt                                                                                                            <  3:165035/149‑1 (MQ=255)
                                    aGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTc                                                                                                          <  1:472949/149‑1 (MQ=255)
                                     gTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCg                                                                                                         <  4:70745/149‑1 (MQ=255)
                                       tCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGc                                                                                                       <  2:430786/149‑1 (MQ=255)
                                          aCGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACg                                                                                                    >  1:276905/1‑149 (MQ=255)
                                          aCGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACg                                                                                                    >  1:178747/1‑149 (MQ=255)
                                           cGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGt                                                                                                   <  3:311451/149‑1 (MQ=255)
                                                  tCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCa                                                                                            >  3:303481/1‑149 (MQ=255)
                                                  tCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCa                                                                                            >  4:525533/1‑149 (MQ=255)
                                                   ccGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCAc                                                                                           <  3:312058/149‑1 (MQ=255)
                                                    cGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACt                                                                                          <  3:108536/149‑1 (MQ=255)
                                                       gatgatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGAt                                                                                       >  1:428548/1‑149 (MQ=255)
                                                       gatgatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGAt                                                                                       >  1:320630/1‑149 (MQ=255)
                                                          gatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCAGTACGTTTAACGATCCtt                                                                                                                                 <  3:252632/104‑1 (MQ=255)
                                                          gatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCAGTACGTTTAACGATCCtt                                                                                                                                 <  3:85223/104‑1 (MQ=255)
                                                          gatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCAGTACGTTTAACGATCCtt                                                                                                                                 >  4:252632/1‑104 (MQ=255)
                                                          gatCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCAGTACGTTTAACGATCCtt                                                                                                                                 >  4:85223/1‑104 (MQ=255)
                                                             cGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGAt                                                                                 >  2:238239/1‑149 (MQ=255)
                                                              gACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATc                                                                                >  1:129415/1‑149 (MQ=255)
                                                                 aaaCGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGt                                                                             <  3:414754/149‑1 (MQ=255)
                                                                 aaaCGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGt                                                                             <  4:308204/149‑1 (MQ=255)
                                                                             gTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGtt                                                                 <  2:87024/149‑1 (MQ=255)
                                                                             gTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGtt                                                                 <  2:201737/149‑1 (MQ=255)
                                                                               ggTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGtttt                                                               <  3:142890/149‑1 (MQ=255)
                                                                                         tCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGc                                                     >  4:7944/1‑149 (MQ=255)
                                                                                         tCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGc                                                     <  3:292918/149‑1 (MQ=255)
                                                                                            ctctTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATc                                                                                <  4:77374/119‑1 (MQ=255)
                                                                                            ctctTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATc                                                                                >  3:77374/1‑119 (MQ=255)
                                                                                               ttaaGCATATAAGTGATGTTATTGGTCCACAGCACCGTTTCGCTAACTGGTGCGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCATACAGCCACATa                                               >  4:134683/4‑149 (MQ=255)
                                                                                                 cAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGc                                             >  4:215784/1‑149 (MQ=255)
                                                                                                    cATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCg                                          <  1:143581/149‑1 (MQ=255)
                                                                                                       ataAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCtt                                       <  3:341350/149‑1 (MQ=255)
                                                                                                          aGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCaa                                    <  1:577222/149‑1 (MQ=255)
                                                                                                          aGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCaa                                    <  4:465341/149‑1 (MQ=255)
                                                                                                           gTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAAt                                   <  1:144434/149‑1 (MQ=255)
                                                                                                            tGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATc                                  <  1:121210/149‑1 (MQ=255)
                                                                                                             gATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCt                                 <  1:97471/149‑1 (MQ=255)
                                                                                                                 ttATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCATGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTAtt                             <  2:442599/149‑1 (MQ=255)
                                                                                                                 ttATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCaa                                    <  4:299944/142‑1 (MQ=255)
                                                                                                                 ttATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCaa                                    >  3:299944/1‑142 (MQ=255)
                                                                                                                             aTCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGGCCTTCAACCACCGCTTTCGGCACGTTATCCCCATGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGtt                 <  4:357473/149‑1 (MQ=255)
                                                                                                                             aTCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGtt                 >  1:22279/1‑149 (MQ=255)
                                                                                                                             aTCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGtt                 >  1:153438/1‑149 (MQ=255)
                                                                                                                             aTCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGtt                 <  2:235052/149‑1 (MQ=255)
                                                                                                                              tCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTa                >  1:417024/1‑149 (MQ=255)
                                                                                                                                aCCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACg              >  3:477156/1‑149 (MQ=255)
                                                                                                                                   gTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGaa           >  2:24431/1‑149 (MQ=255)
                                                                                                                                         cTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGt     <  1:513677/149‑1 (MQ=255)
                                                                                                                                          tAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGCTg    >  2:349328/1‑149 (MQ=255)
                                                                                                                                            aCCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACCGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGTGAc  >  3:453359/1‑149 (MQ=255)
                                                                                                                                            aCCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGTGAc  <  3:525533/149‑1 (MQ=255)
                                                                                                                                            aCCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGTGAc  >  3:29014/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TCGACAAATCTTCCCGTGCCTGCCAGCTCAGCGTTGAGTTCAACGTATACTCCGGGATGATCGACAAACGGTCGCCCGTGGTTTTGTTTTCACTCTTCAGCATATAAGTGATGTTATTGGTCCACATCACCGTTTCGCTAACCGGTACGTTTAACGATCCTTCCAGACCTTCAACCACCGCTTTCGGCACGTTATCCCACTGATAGAGATCGGTGCCGACTGCGTTTTGCCCTACAGCCACATAGCCCGCTTCAATCTTATTGCGATAATCGTTACGGAACCAGGTGAC  >  NC_000913/610568‑610856

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