Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,608 T→C I182M (ATA→ATG yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,764,6080TC100.0% 304.5 / NA 85I182M (ATA→ATGyfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/41);  total (44/41)

GTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCAT  >  NC_000913/2764476‑2764738
                                                                                                                                    |                                                                                                                                  
gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTc                                                                                                                              >  2:1435988/1‑139 (MQ=255)
  catcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGAGCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                            >  2:1179629/1‑139 (MQ=255)
  catcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                            <  2:967422/139‑1 (MQ=255)
       gCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATcc                                                                                                                       <  2:699580/139‑1 (MQ=255)
          cTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTa                                                                                                                    >  2:1408708/1‑139 (MQ=255)
           tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                   >  2:749233/1‑139 (MQ=255)
           tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                   >  2:951668/1‑139 (MQ=255)
           tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                   <  1:1408708/139‑1 (MQ=255)
           tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                   >  1:891269/1‑139 (MQ=255)
           tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                   >  1:1454151/1‑139 (MQ=255)
                  gCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTcc                                                                                                            >  1:1647687/1‑139 (MQ=255)
                       gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGAt                                                                                                                         >  1:893461/1‑121 (MQ=255)
                       gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGAt                                                                                                                         <  2:893461/121‑1 (MQ=255)
                            tGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:572762/139‑1 (MQ=255)
                            tGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:1464560/139‑1 (MQ=255)
                            tGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:891269/139‑1 (MQ=255)
                              cTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAgg                                                                                                >  2:75486/1‑139 (MQ=255)
                                  cccTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                                                                            >  1:448737/1‑139 (MQ=255)
                                          cgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACt                                                                                    <  1:749233/139‑1 (MQ=255)
                                           gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACtt                                                                                   >  1:1729893/1‑139 (MQ=255)
                                               caGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                            <  2:1839120/94‑1 (MQ=255)
                                               caGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                            >  1:1839120/1‑94 (MQ=255)
                                                     cAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGc                                                                         <  1:1435988/139‑1 (MQ=255)
                                                       aGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:72898/112‑1 (MQ=255)
                                                       aGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  >  1:72898/1‑112 (MQ=255)
                                                          aTGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCAt                                                                    <  2:1729893/139‑1 (MQ=255)
                                                            gCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCg                                                                  >  1:42846/1‑139 (MQ=255)
                                                            gCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCg                                                                  >  2:1828744/1‑139 (MQ=255)
                                                              ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGTTGa                                                                                           <  2:267056/112‑1 (MQ=255)
                                                              ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                                           >  1:267056/1‑112 (MQ=255)
                                                               gCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:332586/104‑1 (MQ=255)
                                                               gCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  >  1:332586/1‑104 (MQ=255)
                                                                  ccAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCaa                                                            >  1:1590659/1‑139 (MQ=255)
                                                                        gTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                      <  2:1647687/139‑1 (MQ=255)
                                                                        gTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                      <  2:448737/139‑1 (MQ=255)
                                                                            tCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGc                                                  >  1:139464/1‑139 (MQ=255)
                                                                             cGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    >  1:1608970/1‑88 (MQ=255)
                                                                             cGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    <  2:1608970/88‑1 (MQ=255)
                                                                             cGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                     <  2:1321520/119‑1 (MQ=255)
                                                                             cGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTCCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                     >  1:1321520/1‑119 (MQ=255)
                                                                               gCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                                           >  1:370824/1‑95 (MQ=255)
                                                                               gCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGa                                                                                           <  2:370824/95‑1 (MQ=255)
                                                                                cGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTg                                                                                                                           >  2:126852/1‑62 (MQ=255)
                                                                                cGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTg                                                                                                                           <  1:126852/62‑1 (MQ=255)
                                                                                cGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    <  2:249436/85‑1 (MQ=255)
                                                                                cGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    >  1:249436/1‑85 (MQ=255)
                                                                                       ttCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAg                                       <  2:1454151/139‑1 (MQ=255)
                                                                                                gagaATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa                              >  2:1382857/1‑139 (MQ=255)
                                                                                                gagaATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa                              >  2:585039/1‑139 (MQ=255)
                                                                                                 agaATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                     <  2:1266156/99‑1 (MQ=255)
                                                                                                 agaATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                     >  1:1266156/1‑99 (MQ=255)
                                                                                                    aTGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGt                                      <  2:437501/127‑1 (MQ=255)
                                                                                                    aTGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGt                                      >  1:437501/1‑127 (MQ=255)
                                                                                                      gACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  <  2:1149007/65‑1 (MQ=255)
                                                                                                      gACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                  >  1:1149007/1‑65 (MQ=255)
                                                                                                        ccccATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGGCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAATATTGTc                      >  1:1254845/1‑139 (MQ=255)
                                                                                                        ccccATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTc                      <  1:1179629/139‑1 (MQ=255)
                                                                                                          ccATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                      <  2:528262/105‑1 (MQ=255)
                                                                                                          ccATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                      >  1:528262/1‑105 (MQ=255)
                                                                                                          ccATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCtg                    <  2:1737904/139‑1 (MQ=255)
                                                                                                                  aTCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACt                                                                                    <  1:1410599/67‑1 (MQ=255)
                                                                                                                  aTCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACt                                                                                    >  2:1410599/1‑67 (MQ=255)
                                                                                                                    cTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCtg                    >  2:158031/1‑129 (MQ=255)
                                                                                                                    cTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCtg                    <  1:158031/129‑1 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCa                                                                               >  2:1713409/1‑67 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCa                                                                               <  1:1713409/67‑1 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt                                                                >  2:1518930/1‑82 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt                                                                <  1:1518930/82‑1 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATAc                                            <  2:213748/102‑1 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATAc                                            >  1:213748/1‑102 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCg                <  1:626109/130‑1 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCg                >  2:626109/1‑130 (MQ=255)
                                                                                                                       gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTGTTGAgg      >  1:1462629/1‑139 (MQ=255)
                                                                                                                        ggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    >  2:991125/1‑45 (MQ=255)
                                                                                                                        ggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    <  1:991125/45‑1 (MQ=255)
                                                                                                                        ggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAgg      <  2:1590659/139‑1 (MQ=255)
                                                                                                                        ggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAgg      <  1:1382857/139‑1 (MQ=255)
                                                                                                                            gCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGtt                                 >  2:1621726/1‑108 (MQ=255)
                                                                                                                            gCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGtt                                 <  1:1621726/108‑1 (MQ=255)
                                                                                                                            gCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCAt  <  1:585039/139‑1 (MQ=255)
                                                                                                                            gCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCAt  >  1:1020649/1‑139 (MQ=255)
                                                                                                                             ccAGCCGCATGTTCTTGATCGCTAGAAGTTCCAGTTcgcg                                                                                                    >  1:1103813/1‑40 (MQ=255)
                                                                                                                             ccAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGt                                                                                                         >  2:482782/1‑35 (MQ=255)
                                                                                                                             ccAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGt                                                                                                         <  1:482782/35‑1 (MQ=255)
                                                                                                                             ccAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                                                    <  2:1103813/40‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCAT  >  NC_000913/2764476‑2764738

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