Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,819,488 C→T H205H (CAC→CAT) 
*561* (TAG→TAA) 		 yicG →
ligB ←		 UPF0126 family inner membrane protein
DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,819,4880CT100.0% 236.4 / NA 73H205H (CAC→CAT) 
*561* (TAG→TAA) 		yicG
ligB		UPF0126 family inner membrane protein
DNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/38);  total (35/38)

TGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAG  >  NC_000913/3819361‑3819614
                                                                                                                               |                                                                                                                              
tGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaaa                                                                                                                     >  1:335238/1‑139 (MQ=255)
  cAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                   <  1:1626823/139‑1 (MQ=255)
  cAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                   <  2:1500501/139‑1 (MQ=255)
  cAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                   <  1:217169/139‑1 (MQ=255)
                  aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGAt                                                                                                            >  1:760148/1‑130 (MQ=255)
                  aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGAt                                                                                                            <  2:760148/130‑1 (MQ=255)
                        atgTGGTAATTATCTCCACCCTGGTA‑TCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagccagc                                                                                            >  1:1121912/1‑139 (MQ=255)
                          gtgGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCagccagcc                                                                                           <  1:755611/139‑1 (MQ=255)
                            ggTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                         >  2:30828/1‑139 (MQ=255)
                            ggTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGAGCAGCCAGCCaa                                                                                         >  2:260833/1‑139 (MQ=255)
                                 ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    <  1:745493/139‑1 (MQ=255)
                                 ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    <  2:335238/139‑1 (MQ=255)
                                   aTCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                  >  2:641094/1‑139 (MQ=255)
                                    tctcCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAg                                                                                 >  2:1285271/1‑139 (MQ=255)
                                      tcCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                              <  1:1332865/108‑1 (MQ=255)
                                      tcCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                              >  2:1332865/1‑108 (MQ=255)
                                        cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                   >  2:107317/1‑101 (MQ=255)
                                        cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                   <  1:107317/101‑1 (MQ=255)
                                        cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                               >  1:1323724/1‑121 (MQ=255)
                                        cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGcagcc                                                                                               <  2:1323724/121‑1 (MQ=255)
                                        cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttc                                                                             >  1:33431/1‑139 (MQ=255)
                                          cccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttctt                                                                           >  1:1333247/1‑139 (MQ=255)
                                             tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                   <  2:1210433/112‑1 (MQ=255)
                                             tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                   >  1:1210433/1‑112 (MQ=255)
                                             tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    <  1:1541356/127‑1 (MQ=255)
                                             tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    >  2:1541356/1‑127 (MQ=255)
                                             tGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                        >  1:229821/1‑139 (MQ=255)
                                              ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                   >  1:709452/1‑111 (MQ=255)
                                              ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                   <  2:709452/111‑1 (MQ=255)
                                               gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                        >  2:521946/1‑89 (MQ=255)
                                               gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                        <  1:521946/89‑1 (MQ=255)
                                               gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAttt                                                                      >  1:331443/1‑139 (MQ=255)
                                                tATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTg                                                                     <  2:33431/139‑1 (MQ=255)
                                                 aTTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGAGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgc                                                                    <  1:260833/139‑1 (MQ=255)
                                                   tCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcgc                                                                  <  2:1342692/139‑1 (MQ=255)
                                                    cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaa                                                                                                                      <  2:925315/86‑1 (MQ=255)
                                                    cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCaa                                                                                                                      >  1:925315/1‑86 (MQ=255)
                                                     ggCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCAt                                                                <  1:641094/139‑1 (MQ=255)
                                                          ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                        <  1:1469101/126‑1 (MQ=255)
                                                          ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                        >  2:1469101/1‑126 (MQ=255)
                                                          ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                   >  2:380124/1‑131 (MQ=255)
                                                          ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                   <  1:380124/131‑1 (MQ=255)
                                                          ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                           <  2:229821/139‑1 (MQ=255)
                                                                       cTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                              <  2:1121912/139‑1 (MQ=255)
                                                                       cTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                              <  1:30828/139‑1 (MQ=255)
                                                                         ggCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                              >  2:906008/1‑137 (MQ=255)
                                                                         ggCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                              <  1:906008/137‑1 (MQ=255)
                                                                              tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                       <  1:1646585/107‑1 (MQ=255)
                                                                              tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                       >  2:1646585/1‑107 (MQ=255)
                                                                                cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    >  1:944032/1‑92 (MQ=255)
                                                                                cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    <  2:944032/92‑1 (MQ=255)
                                                                                       aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTcc                              <  1:1285271/139‑1 (MQ=255)
                                                                                            ggATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttc                                                                             >  2:99405/1‑87 (MQ=255)
                                                                                            ggATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttc                                                                             <  1:99405/87‑1 (MQ=255)
                                                                                                  ccGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                       >  2:1386803/1‑87 (MQ=255)
                                                                                                  ccGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                       <  1:1386803/87‑1 (MQ=255)
                                                                                                      tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                         >  2:241181/1‑65 (MQ=255)
                                                                                                      tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                         <  1:241181/65‑1 (MQ=255)
                                                                                                       cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcgc                                                                  >  1:882073/1‑87 (MQ=255)
                                                                                                       cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcgc                                                                  <  2:882073/87‑1 (MQ=255)
                                                                                                       cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                         <  1:1600171/128‑1 (MQ=255)
                                                                                                       cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                         >  2:1600171/1‑128 (MQ=255)
                                                                                                           tactacAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgct          <  2:1333247/139‑1 (MQ=255)
                                                                                                                cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCCGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAc                                             >  1:1414458/1‑99 (MQ=255)
                                                                                                                cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAc                                             <  2:1414458/99‑1 (MQ=255)
                                                                                                                  gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                         >  1:1003983/1‑53 (MQ=255)
                                                                                                                  gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCg                            >  2:1100250/1‑114 (MQ=255)
                                                                                                                  gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCg                            <  1:1100250/114‑1 (MQ=255)
                                                                                                                  accACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGTCAGCCAGCGaa                                                                                         <  2:1003983/52‑1 (MQ=255)
                                                                                                                   ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    >  2:321662/1‑57 (MQ=255)
                                                                                                                   ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                    <  1:321662/57‑1 (MQ=255)
                                                                                                                   ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAg  >  1:1187037/1‑139 (MQ=255)
                                                                                                                   ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAg  >  1:1090362/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                              
TGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAG  >  NC_000913/3819361‑3819614

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: ‑