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,2020GA98.9% 291.2 / ‑3.9 89coding (726/726 nt)
coding (43/678 nt)
yjjQ
bglJ
putative transcriptional regulator
bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (43/45);  minor base . (1/0);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

CAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCG  >  NC_000913/4604076‑4604333
                                                                                                                              |                                                                                                                                   
cAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  1:6010/1‑138 (MQ=255)
cAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  2:6010/138‑1 (MQ=255)
   cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                                      >  1:404336/1‑139 (MQ=255)
     cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  1:84796/133‑1 (MQ=255)
     cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  2:84796/1‑133 (MQ=255)
     cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  1:693917/133‑1 (MQ=255)
     cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  2:693917/1‑133 (MQ=255)
     cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    <  2:404336/139‑1 (MQ=255)
       gAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagtag                                                                                                                  <  1:457832/139‑1 (MQ=255)
         aaTGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                <  1:335653/139‑1 (MQ=255)
           tGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                                              <  1:1039083/139‑1 (MQ=255)
               tGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  1:471074/123‑1 (MQ=255)
               tGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  2:471074/1‑123 (MQ=255)
                    ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    <  1:354236/124‑1 (MQ=255)
                    ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    >  2:354236/1‑124 (MQ=255)
                    ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTgag                                                                                                     >  1:415263/1‑139 (MQ=255)
                     gATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                               <  1:51402/112‑1 (MQ=255)
                     gATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                               >  2:51402/1‑112 (MQ=255)
                              cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    <  2:1152476/114‑1 (MQ=255)
                              cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    >  1:1152476/1‑114 (MQ=255)
                              cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                                           >  2:432060/1‑139 (MQ=255)
                                cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                               >  2:460350/1‑117 (MQ=255)
                                cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                               <  1:460350/117‑1 (MQ=255)
                                cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                                         >  1:1219293/1‑139 (MQ=255)
                                   aGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAg                                                                                      >  2:485500/1‑139 (MQ=255)
                                          aaGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    <  1:1213042/102‑1 (MQ=255)
                                          aaGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    >  2:1213042/1‑102 (MQ=255)
                                           agagAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                              <  2:1219293/139‑1 (MQ=255)
                                            gagaACTTAACTTGCCTTCTGTTTCTGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                             <  1:840851/139‑1 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                     >  2:122905/1‑94 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                     <  1:122905/94‑1 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                               >  1:1450078/1‑100 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                               <  2:1450078/100‑1 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                       >  1:35295/1‑108 (MQ=255)
                                                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                       <  2:35295/108‑1 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    <  1:1378281/87‑1 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                    >  2:1378281/1‑87 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGtt                                                                                                 >  1:389441/1‑106 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGtt                                                                                                 <  2:389441/106‑1 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAg                                                                 >  1:743283/1‑138 (MQ=255)
                                                         gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAg                                                                 <  2:743283/138‑1 (MQ=255)
                                                              cTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCAt                                                           >  2:321491/1‑139 (MQ=255)
                                                               tGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCa                                                                                                                         <  1:1049079/76‑1 (MQ=255)
                                                               tGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCa                                                                                                                         >  2:1049079/1‑76 (MQ=255)
                                                                         gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  2:247694/65‑1 (MQ=255)
                                                                         gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  1:247694/1‑65 (MQ=255)
                                                                           tGCCGCAGAATGGAACGCAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                                           >  1:331700/1‑110 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGACTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAtt                                                                                                                       <  1:926850/66‑2 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATg                                                                                                                       >  2:926850/1‑66 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                              <  1:1489179/107‑1 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                              >  2:1489179/1‑107 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                                           <  2:331700/110‑1 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                                           <  2:875198/110‑1 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                                           >  1:875198/1‑110 (MQ=255)
                                                                           tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                                              >  2:1348139/1‑139 (MQ=255)
                                                                               gCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAcc                                                        >  1:1423941/1‑125 (MQ=255)
                                                                               gCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAcc                                                        <  2:1423941/125‑1 (MQ=255)
                                                                                     tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTa                                                                            <  1:1319863/99‑1 (MQ=255)
                                                                                     tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTa                                                                            >  2:1319863/1‑99 (MQ=255)
                                                                                        aagagaGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAg                                         >  2:875243/6‑131 (MQ=255)
                                                                                        aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAg                                         <  1:875243/131‑1 (MQ=255)
                                                                                        aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCa                                 <  2:650312/139‑1 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          <  2:550460/45‑1 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                          >  1:550460/1‑45 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                        <  2:1211576/47‑1 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                        >  1:1211576/1‑47 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                              >  2:111589/1‑105 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                              <  1:111589/105‑1 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                            >  1:242598/1‑139 (MQ=255)
                                                                                             aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAATGCGTCATGAGTAATATTGGTATTTAGAATTTAATCAGAAAGTTTGCGGGTAAACCTTAAAAGCTCCAAACCTATTCCAGTCAGGAGTCATTTCAGGATgg                           >  1:1227992/1‑138 (MQ=255)
                                                                                               ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagtag                                                                                                                  <  2:910552/51‑1 (MQ=255)
                                                                                               ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagtag                                                                                                                  >  1:910552/1‑51 (MQ=255)
                                                                                               ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                             >  1:1022468/1‑104 (MQ=255)
                                                                                               ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                             <  2:1022468/104‑1 (MQ=255)
                                                                                               ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc                          <  2:242598/139‑1 (MQ=255)
                                                                                                  aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                             <  1:1482214/101‑1 (MQ=255)
                                                                                                  aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                             >  2:1482214/1‑101 (MQ=255)
                                                                                                  aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                            <  1:744535/134‑1 (MQ=255)
                                                                                                  aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                            >  2:744535/1‑134 (MQ=255)
                                                                                                       agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcgcg                  <  1:485500/139‑1 (MQ=255)
                                                                                                              aaaTGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                             <  2:1149942/41‑1 (MQ=255)
                                                                                                              aaaTGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                             >  1:1149942/1‑41 (MQ=255)
                                                                                                              aaaTGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                             >  1:462072/1‑73 (MQ=255)
                                                                                                              aaaTGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                             <  2:462072/73‑1 (MQ=255)
                                                                                                                   tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                                         <  1:1113860/56‑1 (MQ=255)
                                                                                                                   tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                                         >  2:1113860/1‑56 (MQ=255)
                                                                                                                     gCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAt               <  1:96915/128‑1 (MQ=255)
                                                                                                                     gCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAt               >  2:96915/1‑128 (MQ=255)
                                                                                                                       aCTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTgcggcg  <  1:321491/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                   
CAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCG  >  NC_000913/4604076‑4604333

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