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% 247.8 / NA 76H205H (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 (40/36);  total (40/36)

CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621
                                                                                                                                    |                                                                                                                                      
cgcgCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGTCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGt                                                                                                                                  >  2:771499/1‑139 (MQ=255)
  cgcTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTc                                                                                                                                >  2:537472/1‑139 (MQ=255)
   gcTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCa                                                                                                                               <  1:663706/139‑1 (MQ=255)
   gcTGCAACACTATGTTTCTAACCATGAGGTGGAAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGTCATTAAGGTTCa                                                                                                                               >  1:875345/1‑139 (MQ=255)
                      aaCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTggg                                                                                                            >  1:1127270/1‑139 (MQ=255)
                       aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                                      >  2:513503/1‑128 (MQ=255)
                       aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                                      <  1:513503/128‑1 (MQ=255)
                       aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                           >  2:1130932/1‑139 (MQ=255)
                       aCCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGc                                                                                                           >  1:231381/1‑139 (MQ=255)
                                     aTTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGc                                                                                             >  2:296278/1‑139 (MQ=255)
                                      ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            <  2:231381/139‑1 (MQ=255)
                                      ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            <  1:537472/139‑1 (MQ=255)
                                      ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            <  1:296278/139‑1 (MQ=255)
                                      ttATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGACTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            <  1:1055471/139‑1 (MQ=255)
                                          ctcCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                           >  2:806236/1‑104 (MQ=255)
                                          ctcCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                           <  1:806236/104‑1 (MQ=255)
                                             cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTa                                                                                                                                      <  2:488909/90‑1 (MQ=255)
                                             cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTa                                                                                                                                      >  1:488909/1‑90 (MQ=255)
                                             cACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttc                                                                                     >  2:509790/1‑139 (MQ=255)
                                               cccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGcttctt                                                                                   >  2:295377/1‑139 (MQ=255)
                                                ccTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTg                                                                                  >  1:1170891/1‑139 (MQ=255)
                                                 cTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGa                                                                                 <  2:1170891/139‑1 (MQ=255)
                                                   ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                           >  2:1012768/1‑95 (MQ=255)
                                                   ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAAc                                                                                                                           <  1:1012768/95‑1 (MQ=255)
                                                   ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                               >  1:349749/1‑139 (MQ=255)
                                                   ggTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAtt                                                                               <  1:509790/139‑1 (MQ=255)
                                                    gTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAttt                                                                              >  1:966087/1‑139 (MQ=255)
                                                       ttCGTCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTgcg                                                                           <  1:771499/139‑1 (MQ=255)
                                                         cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                                 <  2:458201/115‑1 (MQ=255)
                                                         cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCaa                                                                                                 >  1:458201/1‑115 (MQ=255)
                                                         cGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCa                                                                         >  1:1304294/1‑139 (MQ=255)
                                                          ggCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCAt                                                                        >  1:598561/1‑139 (MQ=255)
                                                               ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                                   <  2:1437952/139‑1 (MQ=255)
                                                               ttttGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                                   <  2:1304294/139‑1 (MQ=255)
                                                                   gCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                               <  2:1393068/139‑1 (MQ=255)
                                                                              ggCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCagccagcc                                                                                                   >  1:889549/1‑92 (MQ=255)
                                                                              ggCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCagccagcc                                                                                                   <  2:889549/92‑1 (MQ=255)
                                                                                  cTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                          >  2:1333345/1‑97 (MQ=255)
                                                                                  cTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCa                                                                                          <  1:1333345/97‑1 (MQ=255)
                                                                                   tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            >  2:1382070/1‑94 (MQ=255)
                                                                                   tACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGcc                                                                                            <  1:1382070/94‑1 (MQ=255)
                                                                                       ccTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                           <  1:295377/139‑1 (MQ=255)
                                                                                          tAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGt                                        >  2:1168849/1‑139 (MQ=255)
                                                                                                 ggATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                                 >  2:466983/1‑139 (MQ=255)
                                                                                                 ggATTACCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                                 <  1:1130932/139‑1 (MQ=255)
                                                                                                      aCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                                   >  1:779465/1‑100 (MQ=255)
                                                                                                      aCCGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttcc                                                                   <  2:779465/100‑1 (MQ=255)
                                                                                                        cGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                           <  2:952736/106‑1 (MQ=255)
                                                                                                        cGGTCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                           >  1:952736/1‑106 (MQ=255)
                                                                                                           tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCgg                      <  2:503263/139‑1 (MQ=255)
                                                                                                           tCTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCgg                      >  1:503263/1‑139 (MQ=255)
                                                                                                            cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                      <  1:697163/107‑1 (MQ=255)
                                                                                                            cTTTTACTACAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                      >  2:697163/1‑107 (MQ=255)
                                                                                                                tactacAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAgct                 <  1:466983/139‑1 (MQ=255)
                                                                                                                 actacAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAgctg                <  1:1168849/139‑1 (MQ=255)
                                                                                                                    acAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                                      >  2:1386153/1‑35 (MQ=255)
                                                                                                                    acAGCCACGAAGGACATTAAGGTTCAAAACCtgtg                                                                                                                      <  1:1386153/35‑1 (MQ=255)
                                                                                                                     cAGCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCCGGCAGctgctg             >  2:941256/1‑139 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                                <  1:1140443/70‑1 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGAt                                                                                >  2:1140443/1‑70 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                               >  1:680730/1‑87 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                               <  2:680730/87‑1 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCCTTCAATAACCTGTCTGGCGCGTCCCGCTCTAGT‑CCCCGGCAGCTGCTGCCa          >  1:1357546/1‑139 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                                 >  2:314906/1‑101 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                                 <  1:314906/101‑1 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCa          >  2:1196213/1‑139 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCa          >  2:1078745/1‑139 (MQ=255)
                                                                                                                       gCCACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCa          >  2:1044426/1‑139 (MQ=255)
                                                                                                                        ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCtgctgc            >  2:57128/1‑136 (MQ=255)
                                                                                                                        ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCtgctgc            <  1:57128/136‑1 (MQ=255)
                                                                                                                        ccACGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAg         >  1:1020260/1‑139 (MQ=255)
                                                                                                                          aCGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAGaa       <  1:279294/139‑1 (MQ=255)
                                                                                                                           cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCtgctgc            <  2:1263103/133‑1 (MQ=255)
                                                                                                                           cGAAGGACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCtgctgc            >  1:1263103/1‑133 (MQ=255)
                                                                                                                               ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAGAACTGCt  <  2:875345/139‑1 (MQ=255)
                                                                                                                               ggACATTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAGAACTGCt  <  2:966087/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CGCGCTGCAACACTATGTTTCTAACCATGATGTGGTAATTATCTCCACCCTGGTATTCGGCTTTTTTGCCCGCTTACTGGCGCTACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGT‑CCCCGGCAGCTGCTGCCAGAACTGCT  >  NC_000913/3819356‑3819621

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