BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,238	G→T	L133I (CTA→ATA)	yjjI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,238	0	G	T	99.8%	1640.5 / ‑4.2	445	coding (397/1551 nt)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (243/201); minor base . (0/1); total (243/202)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 452 total aligned reads displayed.

CGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGG  >  NC_000913/4614148‑4614327

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         <  2:7097386/93‑1 (MQ=255)

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         >  1:2488544/1‑93 (MQ=255)

 gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAt                                                                                        >  2:887818/1‑93 (MQ=255)

  atatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                                       >  2:9114613/1‑93 (MQ=255)

  atatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                                       >  1:234256/1‑93 (MQ=255)

   tatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTc                                                                                      >  2:3002607/1‑93 (MQ=255)

     ttGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                                    >  2:8315752/1‑93 (MQ=255)

     ttGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGGAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                                    >  2:5299503/1‑93 (MQ=255)

       ggCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAc                                                                                  <  2:2243986/93‑1 (MQ=255)

        gCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACa                                                                                 <  2:1207878/93‑1 (MQ=255)

         cGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                                                <  2:3798062/93‑1 (MQ=255)

           tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                              <  2:8592163/93‑1 (MQ=255)

           tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                              <  2:3413088/93‑1 (MQ=255)

              aTAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGt                                                                           <  2:7814778/93‑1 (MQ=255)

               tagaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTCTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                                          <  1:2918693/90‑1 (MQ=255)

               tAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                                          <  2:474105/93‑1 (MQ=255)

                aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTg                                                                         <  1:9380267/93‑1 (MQ=255)

                   ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAGTTTTAACATACGGTTGcaa                                                                      <  2:9561839/93‑1 (MQ=255)

                   ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaa                                                                      >  2:6982818/1‑93 (MQ=255)

                    gCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                     <  1:8133527/93‑1 (MQ=255)

                      gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                                   <  1:4700736/93‑1 (MQ=255)

                      gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                                   <  1:2094435/93‑1 (MQ=255)

                      gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                                   <  2:8145108/93‑1 (MQ=255)

                      gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                                   >  1:2451043/1‑93 (MQ=255)

                       tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                                  <  2:2582091/93‑1 (MQ=255)

                       tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                                  <  2:6459335/93‑1 (MQ=255)

                         tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                                                <  2:9303550/93‑1 (MQ=255)

                          ggCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCa                                                               <  1:2992506/93‑1 (MQ=255)

                            cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATc                                                             <  2:5377179/93‑1 (MQ=255)

                            cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATc                                                             <  1:6307652/93‑1 (MQ=255)

                             aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                            <  1:7842723/93‑1 (MQ=255)

                               ggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAg                                                          >  1:7082764/1‑93 (MQ=255)

                                gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                         <  1:4072727/93‑1 (MQ=255)

                                gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                         >  2:283764/1‑93 (MQ=255)

                                gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                         >  1:6951929/1‑93 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        <  1:6629534/93‑1 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        >  1:528651/1‑93 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        <  2:7671716/93‑1 (MQ=255)

                                  tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                       <  2:1491385/93‑1 (MQ=255)

                                      tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcccc                                                   >  2:4657630/1‑93 (MQ=255)

                                        gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                 <  2:3498390/93‑1 (MQ=255)

                                        gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                 >  1:428262/1‑93 (MQ=255)

                                         aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                                                <  1:1495505/93‑1 (MQ=255)

                                           gTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                              <  2:5473242/93‑1 (MQ=255)

                                             aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGa                                            >  1:1185456/1‑93 (MQ=255)

                                              aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                           >  2:5769961/1‑93 (MQ=255)

                                              aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                           <  2:6720157/93‑1 (MQ=255)

                                               cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                                          <  2:1807669/93‑1 (MQ=255)

                                                 ccAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCgg                                        >  1:9857451/1‑93 (MQ=255)

                                                   aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGct                                      >  1:7317694/1‑92 (MQ=255)

                                                   aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                                      >  1:9004532/1‑93 (MQ=255)

                                                   aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                                      >  1:1933242/1‑93 (MQ=255)

                                                    gtaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAAATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                                     >  2:10017424/3‑93 (MQ=255)

                                                    gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                                     >  1:8627881/1‑93 (MQ=255)

                                                     aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                                    >  1:8914674/1‑93 (MQ=255)

                                                      aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                   >  1:9611475/1‑93 (MQ=255)

                                                      aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                   <  1:6396250/93‑1 (MQ=255)

                                                      aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                   >  2:4839216/1‑93 (MQ=255)

                                                       aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                                  >  1:1291303/1‑93 (MQ=255)

                                                       aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                                  <  2:2737621/93‑1 (MQ=255)

                                                        cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                 >  1:1066508/1‑93 (MQ=255)

                                                        cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                 <  1:6046364/93‑1 (MQ=255)

                                                        cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                 <  1:7497256/93‑1 (MQ=255)

                                                        cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                                 <  2:2004052/93‑1 (MQ=255)

                                                         gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg                                >  1:6171223/1‑93 (MQ=255)

                                                            ttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa                             <  2:4150081/93‑1 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  2:9430523/93‑1 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  2:7158938/93‑1 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGCTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  2:7595811/93‑1 (MQ=255)

                                                               ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg                          >  2:8682918/1‑93 (MQ=255)

                                                                tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGa                         <  2:4888087/93‑1 (MQ=255)

                                                                tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGGTGCCCCAGGTAGACCGGCATCGATGTGAACGa                         >  1:555998/1‑93 (MQ=255)

                                                                  gAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                       >  1:7963673/1‑93 (MQ=255)

                                                                   aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACgg                      <  1:3922626/93‑1 (MQ=255)

                                                                    aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt                     >  2:6063969/1‑93 (MQ=255)

                                                                    aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt                     <  1:8053563/93‑1 (MQ=255)

                                                                     cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTa                    >  1:3572699/1‑93 (MQ=255)

                                                                      aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                   >  1:3841904/1‑93 (MQ=255)

                                                                         gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtg                >  1:4539499/1‑93 (MQ=255)

                                                                          aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg               >  2:5054622/1‑93 (MQ=255)

                                                                          aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg               >  2:4203226/1‑93 (MQ=255)

                                                                          aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg               >  1:7697689/1‑93 (MQ=255)

                                                                           tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt              >  2:3665482/1‑93 (MQ=255)

                                                                              cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTggtggt           >  1:3797291/1‑93 (MQ=255)

                                                                               gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTa          >  2:6344120/1‑93 (MQ=255)

                                                                               gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTa          <  2:8024981/93‑1 (MQ=255)

                                                                                tCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaa         >  1:7204810/1‑93 (MQ=255)

                                                                                tCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaa         <  1:9220063/93‑1 (MQ=255)

                                                                                 cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        >  1:1623747/1‑93 (MQ=255)

                                                                                 cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        >  2:8415944/1‑93 (MQ=255)

                                                                                 cTTGAGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        <  1:8332931/93‑1 (MQ=255)

                                                                                  ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAg       <  1:3346039/93‑1 (MQ=255)

                                                                                   tgtgTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa      <  1:4451792/93‑1 (MQ=255)

                                                                                    ttgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa     <  1:2904879/92‑1 (MQ=255)

                                                                                    gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa     >  2:3364502/1‑93 (MQ=255)

                                                                                    gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa     <  1:10008444/93‑1 (MQ=255)

                                                                                     tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt    <  1:4238133/93‑1 (MQ=255)

                                                                                      gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg   <  1:669811/93‑1 (MQ=255)

                                                                                       ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg  <  2:7280708/93‑1 (MQ=255)

                                                                                       ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg  >  2:4304813/1‑93 (MQ=255)

CGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGG  >  NC_000913/4614148‑4614327

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