Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,517 C→T P129L (CCT→CTT) 
*201* (TAG→TAA) 
tfaP →
tfaE ←
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,5170CT100.0% 300.0 / NA 98P129L (CCT→CTT) 
*201* (TAG→TAA) 
tfaP
tfaE
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (57/41);  total (57/41)

TAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAT  >  NC_000913/1208383‑1208652
                                                                                                                                        |                                                                                                                                       
tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACg                                                                                                                                     >  1:864514/1‑139 (MQ=255)
    aGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTc                                                                                                                                 <  2:1075550/139‑1 (MQ=255)
      aCAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                               <  2:679752/139‑1 (MQ=255)
          tGGCAGGGTAAGGCTGCATTGGGAAGACTTGAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAg                                                                                                                           >  2:1307541/1‑139 (MQ=255)
            gCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtt                                                                                                                         >  1:558654/1‑139 (MQ=255)
                    aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                >  1:1276405/1‑139 (MQ=255)
                          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                          >  1:232388/1‑139 (MQ=255)
                          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                          >  2:273842/1‑139 (MQ=255)
                                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  1:1523008/1‑139 (MQ=255)
                                    aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGt                                                                                                <  2:1255207/139‑1 (MQ=255)
                                    aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGt                                                                                                >  2:735332/1‑139 (MQ=255)
                                     cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                               <  1:273842/139‑1 (MQ=255)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtt                                                                                                                         >  2:1778805/1‑107 (MQ=18)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtt                                                                                                                         <  1:1778805/107‑1 (MQ=18)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                      >  2:546847/1‑139 (MQ=255)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                      >  1:1766574/1‑139 (MQ=255)
                                                tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                                <  2:1120068/96‑1 (MQ=255)
                                                tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                                >  1:1120068/1‑96 (MQ=255)
                                                  aGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  1:1495242/1‑121 (MQ=255)
                                                  aGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  2:1495242/121‑1 (MQ=255)
                                                   gCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGc                                                                                 >  1:1236970/1‑139 (MQ=255)
                                                      gAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGt                                                                              >  1:844564/1‑139 (MQ=255)
                                                             tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                       >  1:1268869/1‑139 (MQ=255)
                                                             tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                       >  2:1747757/1‑139 (MQ=255)
                                                             tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                       >  1:529696/1‑139 (MQ=255)
                                                              ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                      <  2:1523008/139‑1 (MQ=255)
                                                                 aTTTGTGGCTGGATTATCTGGAAGCACCGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                                   >  1:1212665/1‑139 (MQ=255)
                                                                  tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                                  >  2:1664242/1‑139 (MQ=255)
                                                                  tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                                  >  1:1185756/1‑139 (MQ=255)
                                                                   ttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                 <  2:1236970/139‑1 (MQ=255)
                                                                       ggCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                             >  2:1074611/1‑139 (MQ=255)
                                                                        gCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  2:561906/1‑99 (MQ=255)
                                                                        gCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  1:561906/99‑1 (MQ=255)
                                                                        gCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                            >  1:688881/1‑139 (MQ=255)
                                                                         cTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGc                                                                                              <  2:703432/104‑1 (MQ=255)
                                                                         cTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGc                                                                                              >  1:703432/1‑104 (MQ=255)
                                                                            gATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                        >  2:1466719/1‑139 (MQ=255)
                                                                             aTTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                       <  2:1185756/139‑1 (MQ=255)
                                                                             aTTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                       <  1:735332/139‑1 (MQ=255)
                                                                              ttATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                      <  1:769916/139‑1 (MQ=255)
                                                                                 tCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  2:428996/1‑90 (MQ=18)
                                                                                 tCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  1:428996/90‑1 (MQ=18)
                                                                                      aaGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtta                                                                                                                        >  2:1670522/1‑65 (MQ=255)
                                                                                      aaGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtta                                                                                                                        <  1:1670522/66‑2 (MQ=255)
                                                                                      aaGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTaa                                              >  1:209964/1‑139 (MQ=255)
                                                                                          aCTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                 <  1:610471/116‑1 (MQ=255)
                                                                                          aCTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                 >  2:610471/1‑116 (MQ=255)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  1:1517529/80‑1 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  2:1517529/1‑80 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  1:133532/1‑80 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  2:133532/80‑1 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                  <  1:798477/82‑1 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                  >  2:798477/1‑82 (MQ=17)
                                                                                           cTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTg                                         >  1:1354209/1‑139 (MQ=255)
                                                                                            tGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                        <  1:1307541/139‑1 (MQ=255)
                                                                                            tGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                        >  2:1699940/1‑139 (MQ=255)
                                                                                                   ggTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTc                                 >  2:393815/1‑139 (MQ=255)
                                                                                                    gTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                               >  2:496656/1‑76 (MQ=14)
                                                                                                    gTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                               <  1:496656/76‑1 (MQ=14)
                                                                                                    gTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                        <  2:317117/131‑1 (MQ=255)
                                                                                                    gTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                        >  1:317117/1‑131 (MQ=255)
                                                                                                    gTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                                >  1:485799/1‑139 (MQ=255)
                                                                                                      tGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttagg                                                                                                                      <  2:1748033/52‑4 (MQ=14)
                                                                                                      tGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttagg                                                                                                                      >  1:1748033/1‑49 (MQ=17)
                                                                                                       gAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACGTCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                             <  2:864514/139‑1 (MQ=255)
                                                                                                       gAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                             <  2:1212665/139‑1 (MQ=255)
                                                                                                        aT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                 >  1:329487/1‑102 (MQ=255)
                                                                                                        aT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                 <  2:329487/102‑1 (MQ=255)
                                                                                                        aT‑‑ACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct                            >  2:507/1‑139 (MQ=255)
                                                                                                              aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                        <  1:507/139‑1 (MQ=255)
                                                                                                              aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                        <  2:529696/139‑1 (MQ=255)
                                                                                                              aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                        <  1:225517/139‑1 (MQ=255)
                                                                                                              aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                        <  1:1466719/139‑1 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTaa                                                              <  2:1414343/99‑1 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTaa                                                              >  1:1414343/1‑99 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                >  2:1087283/1‑113 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                <  1:1087283/113‑1 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                             >  1:71603/1‑132 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                             <  2:71603/132‑1 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                      >  2:501007/1‑139 (MQ=255)
                                                                                                                ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                      >  2:430424/1‑139 (MQ=255)
                                                                                                                 cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                      >  2:1316149/1‑106 (MQ=255)
                                                                                                                 cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                      <  1:1316149/106‑1 (MQ=255)
                                                                                                                    tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATg                                                                                                              >  1:1079207/1‑47 (MQ=14)
                                                                                                                    tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATg                                                                                                              <  2:1079207/47‑1 (MQ=14)
                                                                                                                    tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    >  2:1054347/1‑57 (MQ=11)
                                                                                                                    tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                    <  1:1054347/57‑1 (MQ=11)
                                                                                                                    tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                  >  2:215987/1‑139 (MQ=255)
                                                                                                                      cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                        <  1:97789/67‑1 (MQ=11)
                                                                                                                      cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                        >  2:97789/1‑67 (MQ=11)
                                                                                                                        ccAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGt              >  1:1441062/1‑139 (MQ=255)
                                                                                                                         cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc             >  2:1699207/1‑139 (MQ=255)
                                                                                                                         cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc             <  2:1441062/139‑1 (MQ=255)
                                                                                                                         cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc             >  2:135244/1‑139 (MQ=255)
                                                                                                                         cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc             >  2:135243/1‑139 (MQ=255)
                                                                                                                         cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc             >  1:97654/1‑139 (MQ=255)
                                                                                                                             tatTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGc         <  1:1664242/139‑1 (MQ=255)
                                                                                                                                    tGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAt  <  2:97654/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAT‑‑ACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAT  >  NC_000913/1208383‑1208652

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