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	88.9%	700.6 / 72.4	271	S382G (AGC→GGC)	yacK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (15/15); new base G (113/128); total (128/143)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 280 total aligned reads displayed.

TGCAACTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACG  >  NC_000913/138134‑138316

tGCAACTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATcaca                                                                                            >  1:5688045/1‑93 (MQ=255)

   aaCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGcc                                                                                         >  1:2561786/1‑93 (MQ=255)

    aCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCa                                                                                        >  1:2069991/1‑93 (MQ=255)

    aCTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGCAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCa                                                                                        >  1:3958004/1‑93 (MQ=255)

      tctctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCTGGATGGATCACGGCCAga                                                                                      >  1:363917/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

         ctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatga                                                                                   >  2:7148778/1‑93 (MQ=255)

         ctATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAgatga                                                                                   >  2:3421044/1‑93 (MQ=255)

           aTGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGatgatg                                                                                 >  1:1020788/1‑93 (MQ=255)

            tGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATgg                                                                                <  2:1384522/93‑1 (MQ=255)

              gACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGc                                                                              <  1:4350628/93‑1 (MQ=255)

               aCCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGcc                                                                             >  2:5950513/1‑93 (MQ=255)

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

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

                   gATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCatat                                                                         >  2:2888810/1‑93 (MQ=255)

                   gATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCatat                                                                         <  2:1457510/93‑1 (MQ=255)

                     tGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATgg                                                                       <  1:7379084/93‑1 (MQ=255)

                     tGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATgg                                                                       >  1:4829478/1‑93 (MQ=255)

                          gATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCAc                                                                  >  2:2357894/1‑93 (MQ=255)

                          gATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCAc                                                                  >  1:1179723/1‑93 (MQ=255)

                           atatGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACg                                                                 >  2:5207250/1‑93 (MQ=255)

                             atgatgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGc                                                               >  2:2178573/1‑93 (MQ=255)

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

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

                               gatgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCaa                                                             <  1:3235572/93‑1 (MQ=255)

                                atgGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAat                                                            <  1:1796229/93‑1 (MQ=255)

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

                                   gggATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATg                                                         >  2:180839/1‑93 (MQ=255)

                                   gggATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATg                                                         >  2:5522222/1‑93 (MQ=255)

                                    ggATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGa                                                        <  2:464762/93‑1 (MQ=255)

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

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

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

                                         cAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCat                                                   >  1:6753909/1‑93 (MQ=255)

                                          aGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCata                                                  >  2:389201/1‑93 (MQ=255)

                                          aGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCata                                                  >  2:5251276/1‑93 (MQ=255)

                                           gATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCatat                                                 <  2:1031023/93‑1 (MQ=255)

                                           gATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCatat                                                 >  2:3716493/1‑93 (MQ=255)

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

                                                tAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAcc                                            <  1:6234864/93‑1 (MQ=255)

                                                tAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAAcc                                            <  1:5395361/93‑1 (MQ=255)

                                                  aTGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAc                                          <  2:7072923/93‑1 (MQ=255)

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

                                                  aTGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAc                                          <  2:7328753/93‑1 (MQ=255)

                                                   tGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcg                                         <  1:895535/93‑1 (MQ=255)

                                                   tGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcg                                         >  2:684628/1‑93 (MQ=255)

                                                   tGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcg                                         <  1:3109688/93‑1 (MQ=255)

                                                   tGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcg                                         >  2:6521500/1‑93 (MQ=255)

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

                                                     gagaAATATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcggc                                       <  1:1776377/93‑1 (MQ=255)

                                                       gaAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCAcggcgg                                     <  1:6530234/93‑1 (MQ=255)

                                                         aaTATGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGa                                   >  2:4688435/1‑93 (MQ=255)

                                                         aaTATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGa                                   >  1:7211577/1‑93 (MQ=255)

                                                          atatGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGaa                                  <  1:7474097/93‑1 (MQ=255)

                                                          atatGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGaa                                  <  1:6131213/93‑1 (MQ=255)

                                                            atGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGt                                <  1:5088942/93‑1 (MQ=255)

                                                            atGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGt                                <  1:6459598/93‑1 (MQ=255)

                                                             tGGCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGGAGtt                               >  2:6606778/1‑93 (MQ=255)

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

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

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

                                                              ggCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTc                              <  1:7082087/93‑1 (MQ=255)

                                                              ggCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATACGAATCATATGAACCACGGCGGGAAGTTc                              >  2:3060672/1‑93 (MQ=255)

                                                               gCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCg                             >  2:5384634/1‑93 (MQ=255)

                                                               gCGATCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCg                             >  1:1584203/1‑93 (MQ=255)

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

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

                                                                  aTCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAtt                          >  2:4345408/1‑93 (MQ=255)

                                                                   tCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAttt                         >  1:104732/1‑93 (MQ=255)

                                                                   tCAGGCGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAttt                         <  1:2712826/93‑1 (MQ=255)

                                                                   tCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGAttt                         <  1:6356437/93‑1 (MQ=255)

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

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

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

                                                                      ggCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTcca                      >  2:570008/1‑93 (MQ=255)

                                                                       gNGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccan                     <  2:5030025/93‑2 (MQ=255)

                                                                        cGATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacc                    <  2:4530267/93‑1 (MQ=255)

                                                                         gATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                   <  2:7292474/93‑1 (MQ=255)

                                                                         gATGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                   <  1:2369483/93‑1 (MQ=255)

                                                                         gATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTccacca                   >  1:1944708/1‑93 (MQ=255)

                                                                          aTGGCCGGGATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCAt                  >  2:4842949/1‑93 (MQ=255)

                                                                                gggATGGATCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCcaac            <  2:1863732/93‑1 (MQ=255)

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

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

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

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

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

                                                                                     ggatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAAt       <  1:6431114/93‑1 (MQ=255)

                                                                                      gatCACGGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATc      <  1:5588838/93‑1 (MQ=255)

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

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

                                                                                         cacgGCCAGATGATGGGCCATATGGAGCACGGCAATATGAATCATATGAACCACGGCGGGAAGATCGATTTCCACCATGCCAACAAAATCAAc   >  2:2015693/1‑93 (MQ=255)

                                                                                         cacaGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAAc   <  1:2862425/93‑1 (MQ=255)

TGCAACTCTCTATGGACCCGATGCTCGATATGATGGGGATGCAGATGCTAATGGAGAAATATGGCGATCAGGCGATGGCCGGGATGGATCACAGCCAGATGATGGGCCATATGGGGCACGGCAATATGAATCATATGAACCACGGCGGGAAGTTCGATTTCCACCATGCCAACAAAATCAACG  >  NC_000913/138134‑138316

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