BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	394,595	T→C	D306G (GAT→GGT)	yaiH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	394,595	0	T	C	99.0%	1047.1 / ‑2.4	312	D306G (GAT→GGT)	yaiH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (154/155); total (156/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 324 total aligned reads displayed.

GTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTG  >  NC_000913/394505‑394684

gTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCAccc                                                                                         <  1:7497487/93‑1 (MQ=255)

gTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCAccc                                                                                         >  2:4541223/1‑93 (MQ=255)

gTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCAccc                                                                                         <  2:6538709/93‑1 (MQ=255)

gTAATGAAACCGCCACCACCGCCAGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCAccc                                                                                         <  2:4763125/93‑1 (MQ=255)

  aaTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGc                                                                                       <  2:1111621/93‑1 (MQ=255)

      aaaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTg                                                                                   <  1:2757504/93‑1 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGc                                                                                  >  1:372070/1‑93 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGc                                                                                  <  2:7764288/93‑1 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGc                                                                                  >  2:3660832/1‑93 (MQ=255)

        aCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGCGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGcc                                                                                 >  1:7348328/1‑93 (MQ=255)

         ccGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCg                                                                                <  1:3027638/93‑1 (MQ=255)

         ccGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCg                                                                                <  2:4634360/93‑1 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgg                                                                               >  1:4474438/1‑93 (MQ=255)

           gCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCggg                                                                              >  2:1514049/1‑93 (MQ=255)

           gCCACCACCGCCTGTCTTCTGAATAAGCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCggg                                                                              <  1:4578298/93‑1 (MQ=255)

            ccaccaCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGTCATATATACCCAACCTAAACAAAGCGCACCCGCTTTGCCggag                                                                             >  2:7086963/1‑91 (MQ=255)

            ccaccaCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgggg                                                                             >  1:7945071/1‑93 (MQ=255)

            ccaccaCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgggg                                                                             <  1:7175570/93‑1 (MQ=255)

              accaccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGAc                                                                           >  1:466742/1‑93 (MQ=255)

               ccaccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACa                                                                          >  2:5015799/1‑93 (MQ=255)

                caccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACAt                                                                         >  2:7128317/1‑93 (MQ=255)

                caccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACAt                                                                         >  2:76689/1‑93 (MQ=255)

                  ccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATcc                                                                       <  1:4827408/93‑1 (MQ=255)

                  ccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATcc                                                                       >  2:5888887/1‑93 (MQ=255)

                   cGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCa                                                                      <  2:401996/93‑1 (MQ=255)

                   cGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCa                                                                      <  2:4399514/93‑1 (MQ=255)

                   cGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCa                                                                      <  2:7870934/93‑1 (MQ=255)

                    gCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCAt                                                                     >  2:2534604/1‑93 (MQ=255)

                      cTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGc                                                                   >  2:390987/1‑93 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGATGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  >  2:6363268/1‑93 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCGAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  <  1:1152334/93‑1 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  <  1:7797489/93‑1 (MQ=255)

                        gTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCa                                                                 >  2:8178205/1‑93 (MQ=255)

                          cttctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAAt                                                               <  2:1841684/93‑1 (MQ=255)

                          cttctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAAt                                                               >  1:634116/1‑93 (MQ=255)

                             ctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCAc                                                            <  1:5322630/93‑1 (MQ=255)

                               gAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACtt                                                          >  1:6626343/1‑93 (MQ=255)

                                aataatCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACttt                                                         <  1:206203/93‑1 (MQ=255)

                                  taatCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCg                                                       <  1:1057769/93‑1 (MQ=255)

                                     tCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTaa                                                    <  2:8426217/93‑1 (MQ=255)

                                      ccccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTaaa                                                   >  1:3859931/1‑93 (MQ=255)

                                       cccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTacac                                                  >  1:8594145/1‑93 (MQ=255)

                                       cccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAAc                                                  <  2:6697566/93‑1 (MQ=255)

                                         cGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTg                                                >  2:6944690/1‑93 (MQ=255)

                                          ggACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgc                                               <  2:6624899/93‑1 (MQ=255)

                                          ggACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgc                                               <  1:227278/93‑1 (MQ=255)

                                           gacgccCCTCTTTCGGGGCCATATATCCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcg                                              <  1:6708128/93‑1 (MQ=255)

                                            acgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgc                                             >  1:4166592/1‑93 (MQ=255)

                                              gacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAACCTgcgcgc                                           >  1:580229/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  2:2052911/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  1:6392908/1‑93 (MQ=255)

                                                cccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATAACTTTCGTAAACTgcgcgcgc                                         <  1:2127761/93‑1 (MQ=255)

                                                 ccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCt                                        <  1:1467915/93‑1 (MQ=255)

                                                 catctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCt                                        <  1:8681959/91‑1 (MQ=255)

                                                   tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGa                                      >  2:2015166/1‑93 (MQ=255)

                                                   tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGa                                      >  1:896554/1‑93 (MQ=255)

                                                   tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGa                                      >  2:4191306/1‑93 (MQ=255)

                                                   tctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGa                                      <  1:8410601/93‑1 (MQ=255)

                                                    ctTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGAt                                     >  1:3663775/1‑93 (MQ=255)

                                                     tttCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATa                                    <  1:6919105/93‑1 (MQ=255)

                                                      ttCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATaa                                   <  1:6768940/93‑1 (MQ=255)

                                                      ttCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATaa                                   <  2:8541525/93‑1 (MQ=255)

                                                      ttCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATaa                                   <  2:7254403/93‑1 (MQ=255)

                                                       tCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATaaa                                  >  2:6468742/1‑93 (MQ=255)

                                                           ggCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCaa                              <  2:2579146/93‑1 (MQ=255)

                                                           ggCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCaa                              <  1:2549004/93‑1 (MQ=255)

                                                            gCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAAc                             <  2:6900052/93‑1 (MQ=255)

                                                             ccATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACg                            >  1:4923461/1‑93 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACCCGTTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          >  1:564300/1‑93 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          <  2:8597481/93‑1 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          <  1:7028012/93‑1 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACACGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTACGCGCGCTGATAAATCAACGTc                          >  2:1218245/1‑93 (MQ=255)

                                                               atatatACCCAACCAAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          <  2:5525972/93‑1 (MQ=255)

                                                                 atataCCCAACCTAAACCAAGCGCCCCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTg                        <  2:4798665/93‑1 (MQ=255)

                                                                 atataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTg                        >  1:8154108/1‑93 (MQ=255)

                                                                 atataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTg                        >  2:4474045/1‑93 (MQ=255)

                                                                  tataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGc                       <  1:3318933/93‑1 (MQ=255)

                                                                    taCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCAt                     <  2:2291751/93‑1 (MQ=255)

                                                                      cccAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGc                   >  1:1412145/1‑93 (MQ=255)

                                                                      cccAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGc                   >  2:1450156/1‑93 (MQ=255)

                                                                        cAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGa                 >  2:5274973/1‑93 (MQ=255)

                                                                         aaCCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGAt                <  1:1945227/93‑1 (MQ=255)

                                                                          aCCTAAACCAAGCGCCCCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATc               <  1:1361038/93‑1 (MQ=255)

                                                                          aCCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATc               <  1:5048107/93‑1 (MQ=255)

                                                                            cTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTg             >  1:4795484/1‑93 (MQ=255)

                                                                            cTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTg             <  2:2703202/93‑1 (MQ=255)

                                                                             caaaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGc            <  2:6311222/92‑1 (MQ=255)

                                                                              aaaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGcc           <  2:2843609/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGGTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCTTGCGATCTGCCt          <  2:1971742/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          <  2:8732932/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          <  1:6496378/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          <  1:3048274/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          <  1:5549451/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          >  1:674628/1‑93 (MQ=255)

                                                                                 ccAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTgg        >  1:6085305/1‑93 (MQ=255)

                                                                                   aaGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGtt      >  2:3486956/1‑93 (MQ=255)

                                                                                   aaGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGtt      >  2:3386439/1‑93 (MQ=255)

                                                                                    aGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTa     <  2:3159629/93‑1 (MQ=255)

                                                                                     gcgcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTAc    >  2:5721485/1‑93 (MQ=255)

                                                                                     gcgcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTAc    >  1:801436/1‑93 (MQ=255)

GTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTG  >  NC_000913/394505‑394684

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