Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,202 G→A *242* (TAG→TAA) 
E15K (GAA→AAA) 		 yjjQ →
bglJ →		 putative transcriptional regulator
bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,2020GA100.0% 254.4 / NA 77*242* (TAG→TAA) 
E15K (GAA→AAA) 		yjjQ
bglJ		putative transcriptional regulator
bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/44);  total (33/44)

GCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCG  >  NC_000913/4604075‑4604315
                                                                                                                               |                                                                                                                 
gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCTGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        <  2:1622727/139‑1 (MQ=255)
   aCGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATg                                                                                                     <  1:1273731/139‑1 (MQ=255)
   aCGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATg                                                                                                     <  2:1201654/139‑1 (MQ=255)
     gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                   <  2:1031339/139‑1 (MQ=255)
     gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                   <  1:1935433/139‑1 (MQ=255)
          aaTGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                              <  1:1671682/139‑1 (MQ=255)
            tGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                            <  1:549230/139‑1 (MQ=255)
                   aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCGGAATAAATAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCAGTATTg                                                                                     <  1:379771/139‑1 (MQ=255)
                   aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                     <  2:1994235/139‑1 (MQ=255)
                     ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                             >  1:180760/1‑129 (MQ=255)
                     ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                             <  2:180760/129‑1 (MQ=255)
                      gATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAaaaaaa                                                                                                              <  2:613333/111‑1 (MQ=255)
                      gATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAaaaaaa                                                                                                              >  1:613333/1‑111 (MQ=255)
                         aaaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGtt                                                                               <  2:1127527/139‑1 (MQ=255)
                           aTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                     <  2:1975480/131‑1 (MQ=255)
                           aTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                     >  1:1975480/1‑131 (MQ=255)
                               cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        <  2:1952953/108‑1 (MQ=255)
                               cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        >  1:1952953/1‑108 (MQ=255)
                               cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                            >  2:1551054/1‑120 (MQ=255)
                               cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                            <  1:1551054/120‑1 (MQ=255)
                               cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                         >  1:492119/1‑139 (MQ=255)
                                 cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                       >  1:329610/1‑139 (MQ=255)
                                    aGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGt                                                                                >  1:1688282/1‑127 (MQ=255)
                                    aGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGt                                                                                <  2:1688282/127‑1 (MQ=255)
                                       gTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGttt                                                                 <  1:487433/139‑1 (MQ=255)
                                         gCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAAAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        <  1:990499/98‑1 (MQ=255)
                                         gCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAAAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        >  2:990499/1‑98 (MQ=255)
                                                aacttaacttGCCTTCTGTTTATGGAGATGCCGTAGGATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                     <  2:1260380/110‑1 (MQ=255)
                                                aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATTAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                     >  1:1260380/1‑110 (MQ=255)
                                                  cttaacttGCCTTCTGTTTTTGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAaaaaaa                                                                                                              >  1:1105519/1‑83 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAaaaaaa                                                                                                              <  2:1105519/83‑1 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                             <  2:3057/84‑1 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                             >  1:3057/1‑84 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  <  1:1048673/95‑1 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  >  2:1048673/1‑95 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAccc                                                      >  1:337542/1‑139 (MQ=255)
                                                  cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAccc                                                      >  1:723349/1‑139 (MQ=255)
                                                       cttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtata                                                 <  2:329610/139‑1 (MQ=255)
                                                           ccTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCt                                             >  1:609403/1‑139 (MQ=255)
                                                            cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            <  1:1660885/139‑1 (MQ=255)
                                                            cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            <  2:337542/139‑1 (MQ=255)
                                                             ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                           <  1:1823384/139‑1 (MQ=255)
                                                               cTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCAt                                         >  2:338849/1‑139 (MQ=255)
                                                                           aTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGt                             <  2:492119/139‑1 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        >  1:768045/1‑63 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        <  2:768045/63‑1 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  <  1:1480294/69‑1 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  >  2:1480294/1‑69 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            >  1:1905050/1‑123 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            >  1:1995958/1‑123 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            <  2:1995958/123‑1 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                            <  2:1905050/123‑1 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                            >  2:923304/1‑139 (MQ=255)
                                                                            tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                            >  1:759859/1‑139 (MQ=255)
                                                                             gCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                       <  2:598253/127‑1 (MQ=255)
                                                                             gCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                       >  1:598253/1‑127 (MQ=255)
                                                                              ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAg                          <  2:723349/139‑1 (MQ=255)
                                                                                gCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                              <  2:750281/69‑1 (MQ=255)
                                                                                gCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                              >  1:750281/1‑69 (MQ=255)
                                                                                      tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATg                                                                                                     >  1:1908468/1‑56 (MQ=255)
                                                                                      tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATg                                                                                                     <  2:1908468/56‑1 (MQ=255)
                                                                                         aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                             <  2:1907491/45‑1 (MQ=255)
                                                                                         aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                             >  1:1907491/1‑45 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        <  2:1414764/45‑1 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                        >  1:1414764/1‑45 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTta                                                   >  2:924335/1‑98 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTta                                                   <  1:924335/98‑1 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCt                                             <  2:924037/104‑1 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCt                                             >  1:924037/1‑104 (MQ=255)
                                                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg          >  1:318129/1‑139 (MQ=255)
                                                                                              aGCCGAATTAAAAAGAGAAATGTCGCACTCATAAAAAAATGCGt                                                                                                         <  1:1706079/44‑1 (MQ=255)
                                                                                              aGCCGAATTAAAAAGAGAAATGTCGCACTCATAAAAAAATGCGt                                                                                                         >  2:1706079/1‑44 (MQ=255)
                                                                                                ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc        <  1:338849/139‑1 (MQ=255)
                                                                                                 cGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCa       <  2:318129/139‑1 (MQ=255)
                                                                                                   aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                           <  1:1882863/85‑1 (MQ=255)
                                                                                                   aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                           >  2:1882863/1‑85 (MQ=255)
                                                                                                      tAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcg  <  2:609403/139‑1 (MQ=255)
                                                                                                          aagaGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  >  2:1433674/1‑39 (MQ=255)
                                                                                                          aagaGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                  <  1:1433674/39‑1 (MQ=255)
                                                                                                          aagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                            >  1:1992526/1‑77 (MQ=255)
                                                                                                          aagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                            <  2:1992526/77‑1 (MQ=255)
                                                                                                                               |                                                                                                                 
GCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCG  >  NC_000913/4604075‑4604315

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