BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	138,226	A→G	S382G (AGC→GGC)	yacK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	138,226	0	A	G	99.6%	887.6 / ‑4.7	260	S382G (AGC→GGC)	yacK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (132/127); total (132/128)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 272 total aligned reads displayed.

CAACTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAC  >  NC_000913/138136‑138315

 aaCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGcc                                                                                        <  1:1765767/93‑1 (MQ=255)

  aCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCa                                                                                       >  2:4701705/1‑93 (MQ=255)

  aCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCa                                                                                       <  2:3681309/93‑1 (MQ=255)

    tctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAga                                                                                     <  2:8733916/93‑1 (MQ=255)

    tctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAga                                                                                     >  2:2496764/1‑93 (MQ=255)

    tctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAga                                                                                     >  1:5446706/1‑93 (MQ=255)

     ctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgat                                                                                    <  1:1500674/93‑1 (MQ=255)

     ctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgat                                                                                    >  2:6590703/1‑93 (MQ=255)

      tctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatg                                                                                   >  2:2382874/1‑93 (MQ=255)

      tctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatg                                                                                   <  1:4793116/93‑1 (MQ=255)

      tctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatg                                                                                   >  1:2990181/1‑93 (MQ=255)

       ctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatga                                                                                  <  1:3923639/93‑1 (MQ=255)

        tATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatgat                                                                                 <  1:1201855/93‑1 (MQ=255)

        tATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatgat                                                                                 >  2:7976564/1‑93 (MQ=255)

         aTGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGatgatg                                                                                <  2:1138465/93‑1 (MQ=255)

         aTGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGatgatg                                                                                <  2:4221029/93‑1 (MQ=255)

               ccGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCat                                                                          >  1:7748004/1‑93 (MQ=255)

               ccGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCat                                                                          <  2:6641344/93‑1 (MQ=255)

                    gCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATggg                                                                     <  2:1432442/93‑1 (MQ=255)

                      tCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGc                                                                   <  2:7719814/93‑1 (MQ=255)

                      tCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGc                                                                   >  1:5921463/1‑93 (MQ=255)

                       cGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCa                                                                  <  1:3738803/93‑1 (MQ=255)

                       cGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCa                                                                  <  2:8095226/93‑1 (MQ=255)

                         atatGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACg                                                                >  2:361635/1‑93 (MQ=255)

                         atatGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACg                                                                >  2:3229926/1‑93 (MQ=255)

                           atgatgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGc                                                              <  2:4388302/93‑1 (MQ=255)

                            tgatgTGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCa                                                             >  1:7579587/1‑93 (MQ=255)

                            tgatgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCa                                                             >  1:7263407/1‑93 (MQ=255)

                               tgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAata                                                          <  2:4909946/93‑1 (MQ=255)

                               tgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAata                                                          <  1:848334/93‑1 (MQ=255)

                               tgGCGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAata                                                          <  1:6702642/93‑1 (MQ=255)

                                ggggATGCAGATTCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGTCCATATGGGGCACGGCAatat                                                         <  2:7322233/93‑1 (MQ=255)

                                 gggATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGCCCAGATGATGGGCAATATGGGGCACGGCAatata                                                        >  1:2696219/1‑92 (MQ=255)

                                  ggATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGa                                                       >  1:8799025/1‑93 (MQ=255)

                                   gATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGaa                                                      <  2:5256422/93‑1 (MQ=255)

                                   gATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGaa                                                      <  2:1960097/93‑1 (MQ=255)

                                     tGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATc                                                    <  1:1659005/93‑1 (MQ=255)

                                     tGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATc                                                    >  2:2106766/1‑93 (MQ=255)

                                      gCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCa                                                   <  2:7612209/93‑1 (MQ=255)

                                      gCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATCGGCCATATGGGGCACGGCAATATGCATCa                                                   >  1:8421199/1‑93 (MQ=255)

                                       gaGATGCTAATGGAGAAATAGCGCGATCAGGCGATGGCCGGAATGCATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCat                                                  <  2:334723/92‑1 (MQ=255)

                                        aGATGCTAATTGAAAAATATGGCGATCCGGCGATGGGCGGGATGGATCGCGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCata                                                 <  1:7340548/93‑1 (MQ=255)

                                        aGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCAAGATGATGGGCCATATGGGGCACGGCAATATGAATCata                                                 <  2:6336354/93‑1 (MQ=255)

                                          aTGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATg                                               <  2:5542877/93‑1 (MQ=255)

                                          aTGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATg                                               >  2:2170437/1‑93 (MQ=255)

                                          aTGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATg                                               >  2:3937922/1‑93 (MQ=255)

                                            gCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGaa                                             >  1:7391880/1‑93 (MQ=255)

                                            gCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGaa                                             <  1:6045900/93‑1 (MQ=255)

                                            gCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGaa                                             >  2:1634703/1‑93 (MQ=255)

                                             cTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAc                                            >  2:3522599/1‑93 (MQ=255)

                                              tAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAcc                                           >  2:5781007/1‑93 (MQ=255)

                                              tAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAca                                           >  1:2385115/1‑92 (MQ=255)

                                              tAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAca                                           >  1:7310553/1‑92 (MQ=255)

                                               aaTGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCa                                          >  2:2484812/1‑93 (MQ=255)

                                               aaTGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCa                                          <  1:8140985/93‑1 (MQ=255)

                                                aTGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAc                                         >  1:8285253/1‑93 (MQ=255)

                                                  ggAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcgg                                       <  1:3362219/93‑1 (MQ=255)

                                                  ggAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcgg                                       >  1:2540456/1‑93 (MQ=255)

                                                     gaAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcggcgg                                    <  1:6499263/93‑1 (MQ=255)

                                                      aaaTATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCggg                                   <  1:5771125/93‑1 (MQ=255)

                                                         tatGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAg                                >  1:8578826/1‑93 (MQ=255)

                                                          atGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGt                               >  2:1715175/1‑93 (MQ=255)

                                                           tGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGtt                              <  1:2033639/93‑1 (MQ=255)

                                                           tGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGtt                              <  1:4833104/93‑1 (MQ=255)

                                                              cGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGa                           <  1:7502561/93‑1 (MQ=255)

                                                               gATCAGGCGATGGCCGGGATGGATCACGGGCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGGTCGAt                          >  1:6154131/1‑93 (MQ=255)

                                                               gATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAt                          <  2:2769096/93‑1 (MQ=255)

                                                               gATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAt                          >  1:7051623/1‑93 (MQ=255)

                                                                aTCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAtt                         <  2:6216295/93‑1 (MQ=255)

                                                                aTCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGCCGGGAAGTTCGttt                         >  1:1101327/1‑93 (MQ=255)

                                                                 tCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCCATATGAATCATATGAACCACGGCGGGAAGTTCGAttt                        >  2:3027886/1‑93 (MQ=255)

                                                                  cAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTc                       <  2:4322240/93‑1 (MQ=255)

                                                                  cAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTc                       <  2:5832134/93‑1 (MQ=255)

                                                                  cAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGCTCGATTTc                       >  2:1192620/1‑93 (MQ=255)

                                                                   aGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTcc                      >  2:5429268/1‑93 (MQ=255)

                                                                   aGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTcc                      <  2:4130809/93‑1 (MQ=255)

                                                                    ggCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATGATAAGAACCACGGCGGGAAGTTCGATTTccc                     >  1:1317845/1‑92 (MQ=255)

                                                                    ggCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTcca                     >  2:12515/1‑93 (MQ=255)

                                                                       gATGGCCGGGATGGATCACGGTCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                  <  2:8673514/93‑1 (MQ=255)

                                                                       gATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                  >  1:6230078/1‑93 (MQ=255)

                                                                          ggCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCAACAATGc               >  1:5230551/1‑93 (MQ=255)

                                                                             cGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaa            >  2:8393651/1‑93 (MQ=255)

                                                                              gggATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaac           >  2:4793901/1‑93 (MQ=255)

                                                                              gggATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaaa           >  1:3460193/1‑92 (MQ=255)

                                                                               ggatggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCCTGCcaaca          >  2:7648525/1‑93 (MQ=255)

                                                                               ggatggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaaca          >  1:4041016/1‑93 (MQ=255)

                                                                                gatggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaacaa         >  2:6130338/1‑93 (MQ=255)

                                                                                gatggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCAATTTCCACCATGCcaacaa         >  1:2368496/1‑93 (MQ=255)

                                                                                 atggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACaaa        <  1:6295549/93‑1 (MQ=255)

                                                                                  tggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACaaaa       >  1:2839845/1‑93 (MQ=255)

                                                                                  tggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACaaaa       <  2:8497147/93‑1 (MQ=255)

                                                                                  tggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTCACCCCAAGCCAACaaaa       >  2:6879819/1‑93 (MQ=255)

                                                                                    gatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATc     >  2:1105346/1‑93 (MQ=255)

                                                                                     atCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCa    >  2:785827/1‑93 (MQ=255)

                                                                                      tCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCaa   >  1:8705486/1‑93 (MQ=255)

                                                                                      tCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCaa   >  1:458247/1‑93 (MQ=255)

                                                                                       cacgGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc  <  1:6647714/93‑1 (MQ=255)

                                                                                       cacgGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc  <  1:3861196/93‑1 (MQ=255)

                                                                                       cacgGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc  <  2:8425673/93‑1 (MQ=255)

                                                                                       cacgGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc  >  1:4872433/1‑93 (MQ=255)

CAACTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAC  >  NC_000913/138136‑138315

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`