Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,855,960 C→T S4N (AGC→AAC) 
*116* (TAG→TAA) 		 yidG ←
yidH ←		 inner membrane protein
DUF202 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,855,9600CT100.0% 257.5 / NA 79S4N (AGC→AAC) 
*116* (TAG→TAA) 		yidG
yidH		inner membrane protein
DUF202 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/41);  total (38/41)

ATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855828‑3856093
                                                                                                                                    |                                                                                                                                     
aTGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  1:1165892/139‑1 (MQ=255)
 tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCTTCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  1:200970/139‑1 (MQ=255)
  gCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGc                                                                                                                               >  2:1620420/1‑139 (MQ=255)
    ccGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCAt                                                                                                                             >  1:1499768/1‑139 (MQ=255)
      gCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           >  2:923940/1‑139 (MQ=255)
        cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATacaa                                                                                                                         >  2:1650531/1‑139 (MQ=255)
            tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATa                                                                                                                            >  2:1411103/1‑132 (MQ=255)
            tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATa                                                                                                                            <  1:1411103/132‑1 (MQ=255)
            tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:923940/139‑1 (MQ=255)
                aGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTc                                                                                                                 <  2:1016811/139‑1 (MQ=255)
                  ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTccc                                                                                                               >  1:41200/1‑139 (MQ=255)
                  ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTccc                                                                                                               >  2:1448531/1‑139 (MQ=255)
                          gATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGAt                                                                                                       <  2:41200/139‑1 (MQ=255)
                           aTAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATc                                                                                                      <  1:1650531/139‑1 (MQ=255)
                            tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  2:733516/1‑111 (MQ=255)
                            tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  1:733516/111‑1 (MQ=255)
                                 aaCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGcc                                                                                                >  2:667934/1‑139 (MQ=255)
                                    gCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:1448531/139‑1 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:1058434/139‑1 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:1620420/139‑1 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:1678393/139‑1 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  2:770992/1‑139 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:25134/139‑1 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGC‑‑CTGGCTTTGCGGTTATCCg                                                                                                                                 >  2:1129540/1‑96 (MQ=255)
                                         cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGC‑‑CTGGCTTTGCGGTTATCCg                                                                                                                                 <  1:1129540/96‑1 (MQ=255)
                                              cgcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAata                                                                                   >  1:954707/1‑139 (MQ=255)
                                                 cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:852712/131‑1 (MQ=255)
                                                 cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  1:852712/1‑131 (MQ=255)
                                                    aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGc                                                                             <  1:770992/139‑1 (MQ=255)
                                                                tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         >  1:1723321/1‑131 (MQ=255)
                                                                tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         <  2:1723321/131‑1 (MQ=255)
                                                                  cGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAg                                                               >  1:1376778/1‑139 (MQ=255)
                                                                     aaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCt                                                                            >  1:204017/1‑123 (MQ=255)
                                                                     aaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCt                                                                            <  2:204017/123‑1 (MQ=255)
                                                                     aaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg                                                            >  1:388984/1‑139 (MQ=255)
                                                                      aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  2:1376778/139‑1 (MQ=255)
                                                                        cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTg                                                         <  2:1499768/139‑1 (MQ=255)
                                                                          caccAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgt                                                       >  1:140591/3‑139 (MQ=255)
                                                                          cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgt                                                       >  2:1373088/1‑139 (MQ=255)
                                                                            ccAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:96352/104‑1 (MQ=255)
                                                                            ccAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  2:96352/1‑104 (MQ=255)
                                                                               gTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           <  1:564719/66‑1 (MQ=255)
                                                                               gTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           >  2:564719/1‑66 (MQ=255)
                                                                                tGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:1266207/95‑1 (MQ=255)
                                                                                tGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  2:1266207/1‑95 (MQ=255)
                                                                                tGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAggg                                                 <  1:1373088/139‑1 (MQ=255)
                                                                                       cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTgg                                                        <  2:373739/125‑1 (MQ=255)
                                                                                       cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTgg                                                        >  1:373739/1‑125 (MQ=255)
                                                                                         cTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:756828/61‑1 (MQ=255)
                                                                                         cTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  2:756828/1‑61 (MQ=255)
                                                                                         cTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:1747284/1‑62 (MQ=255)
                                                                                         cTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:1747284/62‑1 (MQ=255)
                                                                                          tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAc                                                                                           <  1:1418758/87‑1 (MQ=255)
                                                                                          tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAc                                                                                           >  2:1418758/1‑87 (MQ=255)
                                                                                          tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCct                                         >  1:1603482/1‑137 (MQ=255)
                                                                                          tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCct                                         <  2:1603482/137‑1 (MQ=255)
                                                                                            cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  2:673990/1‑83 (MQ=255)
                                                                                            cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:673990/83‑1 (MQ=255)
                                                                                            cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:1701654/83‑1 (MQ=255)
                                                                                            cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  2:1701654/1‑83 (MQ=255)
                                                                                             ggCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:1741082/116‑1 (MQ=255)
                                                                                             ggCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           >  2:1741082/1‑116 (MQ=255)
                                                                                                  cAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACg                               <  2:140591/139‑1 (MQ=255)
                                                                                                        tGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAtt                                                                                 <  1:602065/83‑1 (MQ=255)
                                                                                                        tGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAtt                                                                                 >  2:602065/1‑83 (MQ=255)
                                                                                                         gggTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCt                        <  1:667934/139‑1 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAAc                                                                                                          <  2:248892/51‑1 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAAc                                                                                                          >  1:248892/1‑51 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  1:175013/62‑1 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               >  2:175013/1‑62 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTc                   >  2:1625251/1‑138 (MQ=255)
                                                                                                               gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTc                   <  1:1625251/138‑1 (MQ=255)
                                                                                                                         tGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 <  2:129048/82‑1 (MQ=255)
                                                                                                                         tGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  1:129048/1‑82 (MQ=255)
                                                                                                                          ggCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  2:1165597/51‑1 (MQ=255)
                                                                                                                          ggCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               >  1:1165597/1‑51 (MQ=255)
                                                                                                                               tGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCa  >  1:645151/1‑139 (MQ=255)
                                                                                                                                gCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGTCCATCAATATTAAGc                                                                             <  1:1200598/63‑1 (MQ=255)
                                                                                                                                gCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGc                                                                             >  2:1200598/1‑63 (MQ=255)
                                                                                                                                    |                                                                                                                                     
ATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855828‑3856093

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