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% 228.2 / NA 71*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 (35/36);  total (35/36)

TCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCT  >  NC_000913/4604070‑4604322
                                                                                                                                      |                                                                                                                        
tCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                    <  2:743790/139‑1 (MQ=255)
     gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                               <  2:732730/139‑1 (MQ=255)
     gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                               >  1:24075/1‑139 (MQ=255)
         cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                           >  2:685003/1‑139 (MQ=255)
          gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                          <  2:24075/139‑1 (MQ=255)
           cTGAAATGTATGCATGGATAAAGAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑ATAGAAATGTCGCACTCATAAAATAATGCGTCATGagt                                                                                                         <  2:458015/139‑1 (MQ=255)
               aaTGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                     <  2:229948/139‑1 (MQ=255)
                 tGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                                   <  1:647963/139‑1 (MQ=255)
                  gTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                  <  1:570906/139‑1 (MQ=255)
                    aTGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGt                                                                                                >  2:157106/1‑139 (MQ=255)
                     tGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTa                                                                                               >  2:282220/1‑139 (MQ=255)
                      gCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                               <  1:231341/122‑1 (MQ=255)
                      gCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                               >  2:231341/1‑122 (MQ=255)
                        aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                            <  1:678793/139‑1 (MQ=255)
                        aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                            >  2:678793/1‑139 (MQ=255)
                        aTGGATAAATAGCGCGCAGGGTCCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                            >  1:634838/1‑139 (MQ=255)
                               aaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCa                                                                                                              <  2:43282/114‑1 (MQ=255)
                               aaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCa                                                                                                              >  1:43282/1‑114 (MQ=255)
                                    cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                                >  2:422233/1‑139 (MQ=255)
                                    cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                                >  1:212815/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  2:387545/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  2:38788/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  2:766026/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  2:575321/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  2:612738/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                              >  1:144554/1‑139 (MQ=255)
                                      cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATCGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑ACAGAAATGTCGCACTCATAAAAAACTGCGTCATGAGTAGTATCCGTATTGAGAGGTTATTCCGa                                                                              >  1:447175/1‑139 (MQ=255)
                                                 agagAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                   <  1:343656/139‑1 (MQ=255)
                                                 agagAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                   <  1:344212/139‑1 (MQ=255)
                                                  gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                  <  1:387545/139‑1 (MQ=255)
                                                    gaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                                <  2:554098/139‑1 (MQ=255)
                                                     aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAc                                                               >  1:411538/1‑139 (MQ=255)
                                                     aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAc                                                               >  1:193522/1‑139 (MQ=255)
                                                       cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAc                                                               >  1:637283/1‑139 (MQ=255)
                                                       cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                     >  1:511109/1‑99 (MQ=255)
                                                       cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                     <  2:511109/99‑1 (MQ=255)
                                                           acttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtat                                                         <  2:168797/139‑1 (MQ=255)
                                                            cttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtata                                                        <  2:438670/139‑1 (MQ=255)
                                                              tGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAg                                                      <  1:282220/139‑1 (MQ=255)
                                                               gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGc                                                     <  1:16381/139‑1 (MQ=255)
                                                                 cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                   <  1:422233/139‑1 (MQ=255)
                                                                 cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACGCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                   <  1:575321/139‑1 (MQ=255)
                                                                 cTTCTGTTTATGGACATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACACAGAAAAAAATGCGTCATGAGTAGCATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                   <  2:490705/139‑1 (MQ=255)
                                                                  ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                  >  1:440745/1‑123 (MQ=255)
                                                                  ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                  <  2:440745/123‑1 (MQ=255)
                                                                  ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                  <  2:634838/139‑1 (MQ=255)
                                                                  ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                  <  2:193522/139‑1 (MQ=255)
                                                                  ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                  <  2:53229/139‑1 (MQ=255)
                                                                    cTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCAt                                                >  2:341551/1‑139 (MQ=255)
                                                                      gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                              <  2:546064/139‑1 (MQ=255)
                                                                               gaTGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTGATGag                                                                                                          <  1:601874/70‑1 (MQ=255)
                                                                               gaTGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                          >  2:601874/1‑70 (MQ=255)
                                                                               gaTGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAg                                     >  1:740617/1‑139 (MQ=255)
                                                                                aTGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGt                                    <  1:494334/139‑1 (MQ=255)
                                                                                 tGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                                   >  2:369021/1‑139 (MQ=255)
                                                                                  gCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCa                                  >  1:237698/1‑139 (MQ=255)
                                                                                              aaCACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCa                      <  1:157106/139‑1 (MQ=255)
                                                                                                cacaGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAgg                    >  2:671898/1‑139 (MQ=255)
                                                                                                cacaGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAgg                    <  1:685003/139‑1 (MQ=255)
                                                                                                   aGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                 >  2:561376/1‑139 (MQ=255)
                                                                                                   aGCCGAATTAAGA‑‑AGAGAAATGTCGCACTAATAAAAAAATGCGTCATGAGTAGTATAGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGCAGTCATTTCAGGATg                 >  1:365170/1‑139 (MQ=255)
                                                                                                     ccGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc               <  1:571695/139‑1 (MQ=255)
                                                                                                     ccGAATTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGATTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc               >  2:391616/1‑139 (MQ=255)
                                                                                                        aaTTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                  <  2:45365/53‑1 (MQ=255)
                                                                                                        aaTTAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                  >  1:45365/1‑53 (MQ=255)
                                                                                                           tAAGA‑‑AGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcg         <  2:135913/139‑1 (MQ=255)
                                                                                                                    agaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAtct  >  1:359146/1‑139 (MQ=255)
                                                                                                                     gaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATaa                                                       <  2:615923/85‑1 (MQ=255)
                                                                                                                     gaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATaa                                                       >  1:615923/1‑85 (MQ=255)
                                                                                                                     gaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                  <  1:326581/90‑1 (MQ=255)
                                                                                                                     gaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                  >  2:326581/1‑90 (MQ=255)
                                                                                                                                      |                                                                                                                        
TCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGA‑‑AGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCT  >  NC_000913/4604070‑4604322

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