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% 281.9 / NA 85S4N (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 (41/44);  total (41/44)

CATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855827‑3856093
                                                                                                                                     |                                                                                                                                     
cATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATcc                                                                                                                                  >  2:91449/1‑139 (MQ=255)
 aTGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:682991/139‑1 (MQ=255)
 aTGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  1:516709/139‑1 (MQ=255)
  tGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCgg                                                                                                                                <  1:91449/139‑1 (MQ=255)
     ccGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCAt                                                                                                                             >  2:119093/1‑139 (MQ=255)
        ccTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATaca                                                                                                                          >  2:408436/1‑139 (MQ=255)
             tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:626814/1‑139 (MQ=255)
               ccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAg                                                                                                                   <  1:950561/139‑1 (MQ=255)
                cAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGt                                                                                                                  <  1:66037/139‑1 (MQ=255)
                   ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  2:247630/1‑133 (MQ=255)
                   ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  1:247630/133‑1 (MQ=255)
                      tgtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     <  2:724322/130‑1 (MQ=255)
                      tgtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACAcc                                                                                                                     >  1:724322/1‑130 (MQ=255)
                          tGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGa                                                                                                        >  1:265402/1‑139 (MQ=255)
                                ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCg                                                                                                  >  2:1256135/1‑139 (MQ=255)
                                 cAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGc                                                                                                 <  1:626814/139‑1 (MQ=255)
                                   aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                                                                                               <  1:1256135/139‑1 (MQ=255)
                                    cGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc                                                                                              >  2:1169113/1‑139 (MQ=255)
                                     gCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  2:528244/139‑1 (MQ=255)
                                         tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCa                                                                                         >  1:225909/1‑139 (MQ=255)
                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  1:1062930/1‑139 (MQ=255)
                                          cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  1:119093/139‑1 (MQ=255)
                                                gcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat                                                                                  <  2:225909/139‑1 (MQ=255)
                                                  cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATa                                                                                                            >  1:1117705/1‑111 (MQ=255)
                                                  cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATa                                                                                                            <  2:1117705/111‑1 (MQ=255)
                                                    tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  1:1258644/1‑99 (MQ=255)
                                                    tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  2:1258644/99‑1 (MQ=255)
                                                       cccAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTa                                                                           <  1:1169113/139‑1 (MQ=255)
                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:471974/91‑1 (MQ=255)
                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:1258995/91‑1 (MQ=255)
                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  2:471974/1‑91 (MQ=255)
                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  2:1258995/1‑91 (MQ=255)
                                                            cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                                                      >  2:674513/1‑139 (MQ=255)
                                                               ggTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc                                                                                       <  1:469032/119‑1 (MQ=255)
                                                               ggTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc                                                                                       >  2:469032/1‑119 (MQ=255)
                                                                gTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaat                                                                  >  2:1161250/1‑139 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  1:692375/1‑75 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  2:692375/75‑1 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 <  1:1049189/75‑1 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCg                                                                                                                                 >  2:1049189/1‑75 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  2:1126823/1‑111 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:1126823/111‑1 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         >  2:1346952/1‑131 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                                                                         <  1:1346952/131‑1 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  2:306205/1‑139 (MQ=255)
                                                                 tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  1:380302/1‑139 (MQ=255)
                                                                       aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:1333681/139‑1 (MQ=255)
                                                                       aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  2:265402/139‑1 (MQ=255)
                                                                       aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  1:408436/139‑1 (MQ=255)
                                                                       aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt                                                           <  2:1062930/139‑1 (MQ=255)
                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  2:418155/101‑1 (MQ=255)
                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  1:418155/1‑101 (MQ=255)
                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  1:522721/1‑129 (MQ=255)
                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 <  2:522721/129‑1 (MQ=255)
                                                                           cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgt                                                       >  2:119697/1‑139 (MQ=255)
                                                                             ccAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtgt                                                     <  1:674513/139‑1 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           >  1:341965/1‑60 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATac                                                                                                                           <  2:341965/60‑1 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  2:70353/65‑1 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  1:70353/1‑65 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  1:317050/90‑1 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  2:317050/1‑90 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        >  1:264484/1‑95 (MQ=255)
                                                                                      tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                                                                                        <  2:264484/95‑1 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACCCCAGTCCCATAACGATCa                                                                                                     >  2:668307/1‑80 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCa                                                                                                     <  1:668307/80‑1 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             <  2:1255568/88‑1 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACg                                                                                             >  1:1255568/1‑88 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 >  2:311569/1‑116 (MQ=255)
                                                                                        cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata                                                                 <  1:311569/116‑1 (MQ=255)
                                                                                           tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg                                                            <  2:847077/118‑1 (MQ=255)
                                                                                           tCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg                                                            >  1:847077/1‑118 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      >  2:620891/1‑39 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAacac                                                                                                                      <  1:620891/39‑1 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCa                                                                                                              <  1:754205/47‑1 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCa                                                                                                              >  2:754205/1‑47 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAtt                                                                                 <  2:1090320/76‑1 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAtt                                                                                 >  1:1090320/1‑76 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtg                                                      <  2:1260824/103‑1 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGgtg                                                      >  1:1260824/1‑103 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGt                                            >  2:720313/1‑113 (MQ=255)
                                                                                                                gCAATGCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGt                                            <  1:720313/113‑1 (MQ=255)
                                                                                                                     gCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTa                                                                           >  1:907218/1‑77 (MQ=255)
                                                                                                                     gCGTCTGGCTTTGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTa                                                                           <  2:907218/77‑1 (MQ=255)
                                                                                                                                tGCGGTTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCa  <  1:119697/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CATGCCCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATCGCTTTTTCATTGCGTAACCAGCGCA  >  NC_000913/3855827‑3856093

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