BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,097	A→G	M124V (ATG→GTG)	b3614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,097	0	A	G	97.9%	607.3 / 1.2	194	M124V (ATG→GTG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (93/96); total (97/97)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 199 total aligned reads displayed.

CTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTA  >  NC_000913/3786007‑3786186

ctacgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTg                                                                                         >  2:1183052/1‑93 (MQ=255)

ctacgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTg                                                                                         >  2:3784003/1‑93 (MQ=255)

ctacgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTg                                                                                         >  2:3938378/1‑93 (MQ=255)

      ccGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGgcagc                                                                                   <  1:1849535/93‑1 (MQ=255)

        ggAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCaa                                                                                 >  1:2632511/1‑93 (MQ=255)

         gagaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAg                                                                                >  1:1336432/1‑93 (MQ=255)

          agaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGa                                                                               <  1:4276433/93‑1 (MQ=255)

           gaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGCGGGCAGCAAGAt                                                                              >  1:300064/1‑93 (MQ=255)

              gAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGAc                                                                           <  1:2393145/93‑1 (MQ=255)

                gcagcGACGAAATTGAGCGCATTATTCGTAGTGTGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACct                                                                         >  1:1049751/1‑93 (MQ=255)

                   gcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTa                                                                      <  1:3378343/93‑1 (MQ=255)

                   gcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTa                                                                      <  1:2129419/93‑1 (MQ=255)

                     gacgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAAc                                                                    <  2:23136/93‑1 (MQ=255)

                      acgaAATTGAGTGCATTATTCGTAGTGCGGTCAATAACGTGCTCTATGCCGTGGGGAGCTACTACCACTTGGGCAGTAAGATGTCCTCTATcc                                                                   >  2:4202897/1‑93 (MQ=255)

                      acgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAAcc                                                                   >  1:65277/1‑93 (MQ=255)

                      acgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAAcc                                                                   >  2:1743377/1‑93 (MQ=255)

                       cgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCt                                                                  >  1:3730456/1‑93 (MQ=255)

                       cgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCt                                                                  >  1:525057/1‑93 (MQ=255)

                       cgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCt                                                                  >  1:3156656/1‑93 (MQ=255)

                        gaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTg                                                                 >  1:3280754/1‑93 (MQ=255)

                         aaaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGt                                                                <  1:2931213/93‑1 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               <  1:4739052/93‑1 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               <  2:4246203/93‑1 (MQ=255)

                           aTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGttt                                                              <  1:1889706/93‑1 (MQ=255)

                            ttGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTg                                                             <  2:2054422/93‑1 (MQ=255)

                            ttGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTg                                                             <  2:4460367/93‑1 (MQ=255)

                            ttGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTg                                                             >  1:158645/1‑93 (MQ=255)

                            tggAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTg                                                             <  2:103637/91‑1 (MQ=255)

                             tGATCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTgg                                                            >  1:2251348/1‑93 (MQ=255)

                             tGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTgg                                                            <  1:2789759/93‑1 (MQ=255)

                             tGAGCGCATTATTCGTAGTGCGGTCAATAACGCGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTgg                                                            <  2:380772/93‑1 (MQ=255)

                              gAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTCACCTGTTTGGt                                                           >  2:2779340/1‑93 (MQ=255)

                              gAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGt                                                           >  1:3649908/1‑93 (MQ=255)

                               aGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTa                                                          <  2:863931/93‑1 (MQ=255)

                                  gcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGc                                                       >  2:2941340/1‑93 (MQ=255)

                                  gcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGc                                                       <  1:3699461/93‑1 (MQ=255)

                                  gcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGc                                                       <  2:1417416/93‑1 (MQ=255)

                                   cATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCa                                                      >  2:2277678/1‑93 (MQ=255)

                                   cATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCa                                                      <  2:1577480/93‑1 (MQ=255)

                                    attattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAg                                                     >  2:2704069/1‑93 (MQ=255)

                                       attCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaaa                                                  >  2:1625818/1‑93 (MQ=255)

                                       attCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaaa                                                  <  2:1981430/93‑1 (MQ=255)

                                       attCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaaa                                                  <  1:3583216/93‑1 (MQ=255)

                                       attCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaaa                                                  <  1:1155245/93‑1 (MQ=255)

                                        ttCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAg                                                 <  1:2260585/93‑1 (MQ=255)

                                        ttCGTAGTGCGGTCAAAAACGTGCCCTATGCCGTGAGGATCAACATCCACGTGGGCAGAAAGATGACCTCTAACCTGTTTGGTATGCAGAAAg                                                 <  1:4216826/93‑1 (MQ=255)

                                          cGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTg                                               >  1:3020824/1‑93 (MQ=255)

                                            tAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGAt                                             <  2:788355/93‑1 (MQ=255)

                                            tAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGAt                                             >  1:3110628/1‑93 (MQ=255)

                                              gTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGc                                           >  1:4042157/1‑93 (MQ=255)

                                               tGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGAGAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGCTgag                                          >  1:2584219/1‑91 (MQ=255)

                                                 gggTCAATATGGTGCCCTATGCCGTGGGTATCAACGACCAGATGGGCAGCAAGATGACCGCTAAGCTGTTTGGTATGCAGAAAGTGATGCAgg                                        <  2:1526118/92‑1 (MQ=255)

                                                  ggTCAATAGCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGgc                                       <  1:4518308/93‑1 (MQ=255)

                                                  ggTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGgc                                       >  1:4383013/1‑93 (MQ=255)

                                                  ggTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGgc                                       <  1:3458090/93‑1 (MQ=255)

                                                    tCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGgcgc                                     <  1:2149915/93‑1 (MQ=255)

                                                      aaTAACGTGCCCTATGCGGTAGGGATCAACAACCACGGGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   <  2:875715/93‑1 (MQ=255)

                                                      aaTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   >  2:1362661/1‑93 (MQ=255)

                                                      aaTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   >  2:1749753/1‑93 (MQ=255)

                                                        tAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGa                                 >  1:2018048/1‑93 (MQ=255)

                                                        tAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGa                                 >  2:333749/1‑93 (MQ=255)

                                                         aaCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAg                                <  2:1374096/93‑1 (MQ=255)

                                                         aaCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAg                                >  1:3880417/1‑93 (MQ=255)

                                                         aaCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAg                                <  1:56572/93‑1 (MQ=255)

                                                          aCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGc                               >  2:1929699/1‑93 (MQ=255)

                                                          aCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGc                               <  1:3075516/93‑1 (MQ=255)

                                                           cGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCg                              <  2:1593783/93‑1 (MQ=255)

                                                           cGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCg                              <  1:4403080/93‑1 (MQ=255)

                                                            gTGCCCTATGCCGTGGGGATCAACAACCACGTGTGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGt                             <  2:3766029/93‑1 (MQ=255)

                                                            gTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGt                             <  1:397201/93‑1 (MQ=255)

                                                              gCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTa                           <  1:464065/93‑1 (MQ=255)

                                                              gCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTa                           >  1:3334505/1‑93 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         >  2:4358074/1‑93 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         >  1:763822/1‑93 (MQ=255)

                                                                 cTATGCCGTGGGGATCAACAACCACGTGGACAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACaa                        >  2:363154/1‑93 (MQ=255)

                                                                  tATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAAt                       <  1:446224/93‑1 (MQ=255)

                                                                  tATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAAt                       <  2:2752889/93‑1 (MQ=255)

                                                                    tGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCt                     >  1:3221771/1‑93 (MQ=255)

                                                                        gTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTACCCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTAc                 <  2:1445177/93‑1 (MQ=255)

                                                                        gTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTAc                 >  2:3012540/1‑93 (MQ=255)

                                                                         tGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACt                <  2:3531447/93‑1 (MQ=255)

                                                                         tGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACt                >  1:1427159/1‑93 (MQ=255)

                                                                          ggggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACtt               <  1:4165629/93‑1 (MQ=255)

                                                                          ggggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACtt               >  1:1617329/1‑93 (MQ=255)

                                                                          ggggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACtt               <  1:3086464/93‑1 (MQ=255)

                                                                           gggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTc              >  1:923958/1‑93 (MQ=255)

                                                                           gggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTc              <  2:2851231/93‑1 (MQ=255)

                                                                             gATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCt            >  2:578483/1‑93 (MQ=255)

                                                                             gATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCt            <  1:3813263/93‑1 (MQ=255)

                                                                             gATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCt            <  2:3330525/93‑1 (MQ=255)

                                                                              aTCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTc           <  2:1839177/93‑1 (MQ=255)

                                                                               tCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCg          <  1:1300244/93‑1 (MQ=255)

                                                                                caacaaCCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGa         <  2:3371853/93‑1 (MQ=255)

                                                                                caacaaCCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGa         <  2:2128918/93‑1 (MQ=255)

                                                                                 aacaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGAc        <  1:79866/93‑1 (MQ=255)

                                                                                  acaacCACGTGTGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACa       <  1:728236/93‑1 (MQ=255)

                                                                                    aacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGc     <  1:1706546/93‑1 (MQ=255)

                                                                                     acCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCg    >  1:2969137/1‑93 (MQ=255)

                                                                                       cacgTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTa  >  2:2488996/1‑93 (MQ=255)

                                                                                       cacgTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTa  >  2:4709855/1‑93 (MQ=255)

CTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTA  >  NC_000913/3786007‑3786186

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