Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,998 (A)7→6 pseudogene (359/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,9920A.100.0% 233.8 / NA 56pseudogene (353/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (25/31);  total (25/31)

TATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTG  >  NC_000913/528854‑529122
                                                                                                                                          |                                                                                                                                  
tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCTAGATACa                                                                                                                                    >  1:868485/1‑139 (MQ=255)
tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                                                                                                                    >  1:1332824/1‑139 (MQ=255)
tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                                                                                                                    >  1:1405946/1‑139 (MQ=255)
tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                                                                                                                    >  1:600826/1‑139 (MQ=255)
tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                                                                                                                    >  1:812439/1‑139 (MQ=255)
  taCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTTACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaa                                                                                                                                  >  1:684966/1‑139 (MQ=255)
     acaAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCTGATACaacaaa                                                                                                                               >  1:979547/1‑139 (MQ=255)
     acaAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaa                                                                                                                               >  2:1208014/1‑139 (MQ=255)
     acaAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGAAGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGCGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaa                                                                                                                               >  2:389182/1‑139 (MQ=255)
       aaaTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagt                                                                                                                             >  1:1483435/1‑137 (MQ=255)
        aaTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            <  2:1405946/139‑4 (MQ=255)
         aTGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtgt                                                                                                                           >  2:971555/1‑135 (MQ=255)
          tGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGtgtg                                                                                                                         >  2:1058100/1‑139 (MQ=255)
            tAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTgtggt                                                                                                                       <  2:684966/139‑1 (MQ=255)
               aagaTGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTa                                                                                                                    >  2:770425/1‑139 (MQ=255)
                   tgatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            <  2:783989/128‑4 (MQ=255)
                   tgatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            >  1:783989/1‑125 (MQ=255)
                   tgatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTgtggt                                                                                                                       <  2:1179273/132‑1 (MQ=255)
                   tgatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAata                                                                                                                <  1:389182/139‑1 (MQ=255)
                   tgatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATAAAAGAATGATCCAGATACAAAAAAGTgtggt                                                                                                                       >  1:1179273/1‑132 (MQ=255)
                    gatTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAatat                                                                                                               >  2:1428699/1‑139 (MQ=255)
                      ttattTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATgg                                                                                                             >  1:374944/1‑139 (MQ=255)
                                      tttCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCaa                                                                                             <  1:1399954/139‑1 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAata                                                                                                                <  1:150394/115‑1 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAata                                                                                                                >  2:150394/1‑115 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGa                                                                                        <  1:1058100/139‑1 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGa                                                                                        <  2:374944/139‑1 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGa                                                                                        <  1:1428699/139‑1 (MQ=255)
                                           cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGAGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGa                                                                                        <  2:979547/139‑1 (MQ=255)
                                                  atgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTgtggt                                                                                                                       <  2:555256/101‑1 (MQ=255)
                                                  atgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTgtggt                                                                                                                       >  1:555256/1‑101 (MQ=255)
                                                    gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTaaa                                                                               >  2:652467/1‑139 (MQ=255)
                                                        atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAg                                                                                                  >  2:138989/1‑116 (MQ=255)
                                                        atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAg                                                                                                  <  1:138989/116‑1 (MQ=255)
                                                        atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCtt                                                                           >  2:741514/1‑139 (MQ=255)
                                                        atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCtt                                                                           >  2:84540/1‑139 (MQ=255)
                                                          gggggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCtttt                                                                         >  2:225094/1‑139 (MQ=255)
                                                                     tAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTg                                                              >  1:623414/1‑139 (MQ=255)
                                                                      aaGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGa                                                             >  2:161529/1‑139 (MQ=255)
                                                                           aaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTa                                                        <  1:652467/139‑1 (MQ=255)
                                                                             cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGt                                                               <  2:11170/130‑1 (MQ=255)
                                                                             cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGt                                                               >  1:11170/1‑130 (MQ=255)
                                                                             cATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAAt                                                      >  1:45781/1‑139 (MQ=255)
                                                                              aTTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATg                                                     <  2:812439/139‑1 (MQ=255)
                                                                                   cAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACGAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACg                                                <  1:801772/139‑1 (MQ=255)
                                                                                     gTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                              <  2:1332824/139‑1 (MQ=255)
                                                                                     gTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                              <  1:84540/139‑1 (MQ=255)
                                                                                      tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTttg                                                                    <  2:74517/116‑1 (MQ=255)
                                                                                      tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTttg                                                                    >  1:74517/1‑116 (MQ=255)
                                                                                      tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCat                                             <  2:1483435/139‑1 (MQ=255)
                                                                                        aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGt                                                               >  2:1273733/1‑119 (MQ=255)
                                                                                        aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGt                                                               <  1:1273733/119‑1 (MQ=255)
                                                                                        aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTCTGTTGGTGAAATTAATGTTACGCatta                                           >  2:869079/1‑139 (MQ=255)
                                                                                          ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCattatt                                         <  1:741514/139‑1 (MQ=255)
                                                                                                ggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            <  2:12950/51‑4 (MQ=255)
                                                                                                ggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            >  1:12950/1‑48 (MQ=255)
                                                                                                  aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            >  1:794285/1‑46 (MQ=255)
                                                                                                  aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaaaaaagtg                                                                                                                            <  2:794285/49‑4 (MQ=255)
                                                                                                  aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCat                                             >  2:581299/1‑127 (MQ=255)
                                                                                                  aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                 <  1:971555/139‑1 (MQ=255)
                                                                                                  aaCAAGTGGCGAAAATCGTATCAAAGAATGAGCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATAACAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCat                                             <  1:581299/127‑1 (MQ=255)
                                                                                                       gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                            <  2:868485/139‑1 (MQ=255)
                                                                                                       gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                            <  1:773302/139‑1 (MQ=255)
                                                                                                       gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                            >  2:773302/1‑139 (MQ=255)
                                                                                                          gCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                         >  2:1458816/1‑139 (MQ=255)
                                                                                                          gCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                         >  1:229016/1‑139 (MQ=255)
                                                                                                           tgAAAATCGTATCAAAGAATGATCCAGATACAAAAAAGTGTGGGGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                        >  2:1328513/2‑139 (MQ=255)
                                                                                                                       cAAAGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt            <  1:770425/139‑1 (MQ=255)
                                                                                                                         aaGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa          <  2:600826/139‑1 (MQ=255)
                                                                                                                         aaGAATGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa          <  2:45781/139‑1 (MQ=255)
                                                                                                                            aaTGATCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc       <  1:302436/139‑1 (MQ=255)
                                                                                                                                 tCCAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg  <  2:781352/139‑1 (MQ=255)
                                                                                                                                   cAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                              >  2:825140/1‑93 (MQ=255)
                                                                                                                                   cAGATACAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                              <  1:825140/93‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
TATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTG  >  NC_000913/528854‑529122

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