BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,431,925	A→G	I79V (ATC→GTC)	b4212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,431,925	0	A	G	98.0%	604.5 / 0.6	195	I79V (ATC→GTC)	b4212	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/4); new base G (87/104); total (87/108)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 202 total aligned reads displayed.

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

aGCGGCCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTc                                                                                         <  1:1225696/93‑1 (MQ=255)

  cGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTGCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgc                                                                                       <  2:270731/93‑1 (MQ=255)

  cGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgc                                                                                       <  2:3485801/93‑1 (MQ=255)

  cGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgc                                                                                       <  1:1416690/93‑1 (MQ=255)

  cGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgc                                                                                       >  2:2852601/1‑93 (MQ=255)

    aCCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGt                                                                                     >  2:226147/1‑93 (MQ=255)

        gcggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATcc                                                                                 <  1:1247948/93‑1 (MQ=255)

        gcggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATcc                                                                                 <  2:3032779/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              <  1:2972198/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              <  1:733659/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggg                                                                              >  2:2270373/1‑92 (MQ=255)

            cgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtg                                                                             <  1:347596/93‑1 (MQ=255)

             gcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  1:1743878/93‑1 (MQ=255)

             gcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  2:1453790/93‑1 (MQ=255)

             gcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  2:3108829/93‑1 (MQ=255)

             gcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  2:4352305/93‑1 (MQ=255)

              cAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtggt                                                                           <  1:697628/93‑1 (MQ=255)

              cAAAATCGGTGGGGTGAGTAAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtggt                                                                           <  2:1350394/93‑1 (MQ=255)

                aaaTCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgc                                                                         <  2:3004251/93‑1 (MQ=255)

                aaaTCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgc                                                                         <  1:3365335/93‑1 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      <  1:4352734/93‑1 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      >  1:3714434/1‑93 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      <  2:1824665/93‑1 (MQ=255)

                    cGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgcc                                                                     <  2:1158724/93‑1 (MQ=255)

                    cGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgcc                                                                     <  2:4560222/93‑1 (MQ=255)

                     ggTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGccgccg                                                                    >  2:2099454/1‑93 (MQ=255)

                           gtgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtgtg                                                              >  2:4416727/1‑93 (MQ=255)

                            tgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTgg                                                             <  2:4312185/93‑1 (MQ=255)

                              agtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           >  1:1291599/1‑93 (MQ=255)

                              agtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           <  2:1247027/93‑1 (MQ=255)

                              agtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           >  1:1600235/1‑93 (MQ=255)

                               gtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAat                                                          >  2:4490585/1‑93 (MQ=255)

                                tgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAata                                                         >  1:3907652/1‑93 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        <  1:2392/93‑1 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        <  2:46698/93‑1 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        <  1:610554/93‑1 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        <  1:3964438/93‑1 (MQ=255)

                                  aaaaGATGGTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAAtata                                                       <  1:1554333/93‑1 (MQ=255)

                                  aaaaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAAtata                                                       <  2:3550953/93‑1 (MQ=255)

                                    aaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGc                                                     >  2:4185663/1‑93 (MQ=255)

                                    aaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGc                                                     >  1:1787363/1‑93 (MQ=255)

                                     aGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGcc                                                    <  2:3064551/93‑1 (MQ=255)

                                      gATGCTTGCGCAGTCGTTACAGGCGTTAGAACATGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCt                                                   <  1:4593245/93‑1 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  <  2:1047961/93‑1 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  >  1:304929/1‑93 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  <  2:472479/93‑1 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 <  1:3174638/93‑1 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 <  1:23144/93‑1 (MQ=255)

                                         gCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAc                                                <  2:3247717/93‑1 (MQ=255)

                                         gCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAc                                                >  1:1683571/1‑93 (MQ=255)

                                           ttGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgc                                              <  1:2733385/93‑1 (MQ=255)

                                           ttGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgc                                              >  1:5761/1‑93 (MQ=255)

                                             gcgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGGATATAGCCTCAcgccg                                            >  2:3907913/1‑93 (MQ=255)

                                             gcgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccg                                            >  1:157905/1‑93 (MQ=255)

                                               gcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCt                                          <  1:2788033/93‑1 (MQ=255)

                                                cAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTg                                         <  2:2650837/93‑1 (MQ=255)

                                                cAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTg                                         <  1:2157038/93‑1 (MQ=255)

                                                 aGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTgg                                        <  1:1031362/93‑1 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       >  2:2941061/1‑93 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       >  1:441845/1‑93 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       >  1:2454075/1‑93 (MQ=255)

                                                   tCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGc                                      <  2:724817/93‑1 (MQ=255)

                                                   tCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGc                                      <  1:3125256/93‑1 (MQ=255)

                                                    cGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCg                                     >  2:1570722/1‑93 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    <  2:4595000/93‑1 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:2476157/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   <  2:1462676/93‑1 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGCGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:165546/1‑93 (MQ=255)

                                                       tACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGc                                  <  2:3595291/93‑1 (MQ=255)

                                                       tACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGc                                  <  2:57456/93‑1 (MQ=255)

                                                        aCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCa                                 <  1:616000/93‑1 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               <  2:1999894/93‑1 (MQ=255)

                                                           ggCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGt                              >  1:2042566/1‑93 (MQ=255)

                                                           ggCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGt                              >  1:1867695/1‑93 (MQ=255)

                                                           ggCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGt                              <  1:2332351/93‑1 (MQ=255)

                                                           ggCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGt                              >  2:616815/1‑93 (MQ=255)

                                                             cGTTAGAACAGGATCGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGa                            <  1:9809/93‑1 (MQ=255)

                                                              gTTAGTACAGGATGGTTTTCTTAACCGTGTCGCGTATTCGGTGGTGCCGCCGCATGTGGAATATAGCCTCATGCCGCTGGTCGAGCAGGTGAg                           >  2:4337240/1‑93 (MQ=255)

                                                              gTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAg                           >  1:1109341/1‑93 (MQ=255)

                                                               ttAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGc                          >  2:1535369/1‑93 (MQ=255)

                                                               ttAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGc                          >  2:1072595/1‑93 (MQ=255)

                                                                 aGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGa                        >  1:1915377/1‑93 (MQ=255)

                                                                 aGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGa                        >  1:2025289/1‑93 (MQ=255)

                                                                  gAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaa                       >  1:220449/1‑93 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     <  1:4381758/93‑1 (MQ=255)

                                                                    aCAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaaaa                     >  1:4690490/1‑93 (MQ=255)

                                                                      aggaggGTTTTCGTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAgg                   <  1:4733725/93‑1 (MQ=255)

                                                                       ggATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGt                  <  2:307031/93‑1 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGCATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 >  1:1731640/1‑93 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 <  2:243749/93‑1 (MQ=255)

                                                                        gATGGTTTTCTTAACCGGGTCGCGTATCCGGTGGTTCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 <  2:664327/93‑1 (MQ=255)

                                                                          tGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgc               <  1:2464873/93‑1 (MQ=255)

                                                                            gTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCTAGCAGGTGAGCGAAAAGGTtgctg             <  2:2909436/93‑1 (MQ=255)

                                                                               ttCTTAACGGTGTCGCGTATTCGGTGGTGCCGCCGCATGTGGAGTATAGCCTAACGCGGCTCGGCCAACAGGTGAACGAAAGGGGTGCCGCAc          >  2:4072991/1‑93 (MQ=255)

                                                                               ttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCAc          <  1:3275006/93‑1 (MQ=255)

                                                                                   tAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATGGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGc      <  2:1755072/93‑1 (MQ=255)

                                                                                   tAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGc      >  2:3317611/1‑93 (MQ=255)

                                                                                    aaCCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCAATGGCg     >  2:4686497/1‑93 (MQ=255)

                                                                                       cGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAc  >  1:835213/1‑93 (MQ=255)

                                                                                       cGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAc  <  1:2361694/93‑1 (MQ=255)

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

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