BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,837,057	G→A	L188L (CTC→CTT)	nlpA ←	lipoprotein‑28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,837,057	0	G	A	98.6%	737.3 / ‑2.2	214	L188L (CTC→CTT)	nlpA	lipoprotein‑28
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (96/115); total (97/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 219 total aligned reads displayed.

GCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTT  >  NC_000913/3836967‑3837146

gCTGTCGTGCACCGGAGAAAGCCCGTTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCACCCAGTACCCGCGGGAg                                                                                         <  2:2732372/93‑1 (MQ=255)

gCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAg                                                                                         <  2:1901426/93‑1 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGc                                                                                        <  1:2962030/93‑1 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGc                                                                                        >  2:1357234/1‑93 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGc                                                                                        <  2:4545340/93‑1 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGc                                                                                        <  1:2210596/93‑1 (MQ=255)

  tGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCt                                                                                       >  1:3859003/1‑93 (MQ=255)

   gTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTg                                                                                      <  2:3274335/93‑1 (MQ=255)

         cACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGcccc                                                                                >  1:1948446/1‑93 (MQ=255)

          aCCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGccccc                                                                               >  1:3668741/1‑93 (MQ=255)

          aCCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGccccc                                                                               <  1:2775777/93‑1 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCt                                                                              <  1:740139/93‑1 (MQ=255)

            cGGAGAAAGCTCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTc                                                                             <  1:1908248/93‑1 (MQ=255)

            cGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTc                                                                             <  1:2879232/93‑1 (MQ=255)

             ggAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCg                                                                            >  1:310123/1‑93 (MQ=255)

              gagaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGa                                                                           >  1:3780914/1‑93 (MQ=255)

                gaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGt                                                                         >  2:668675/1‑93 (MQ=255)

                 aaaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGtt                                                                        >  2:1760572/1‑93 (MQ=255)

                  aaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTc                                                                       >  1:3496117/1‑93 (MQ=255)

                    gCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCa                                                                     >  1:266109/1‑93 (MQ=255)

                    gCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCa                                                                     <  2:868265/93‑1 (MQ=255)

                     ccgggTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCAt                                                                    <  2:4637687/90‑1 (MQ=255)

                     cccGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCAt                                                                    <  1:1583346/93‑1 (MQ=255)

                      ccGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATa                                                                   >  2:2339158/1‑93 (MQ=255)

                        ggTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAAt                                                                 >  1:4665745/1‑93 (MQ=255)

                          tCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCt                                                               <  1:4468343/93‑1 (MQ=255)

                          tCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCt                                                               <  1:1241960/93‑1 (MQ=255)

                          tCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCt                                                               >  1:3460118/1‑93 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGc                                                             <  2:418130/93‑1 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGc                                                             >  1:177388/1‑93 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGc                                                             >  2:2564003/1‑93 (MQ=255)

                             gctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCa                                                            <  1:2581355/93‑1 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCaaa                                                          >  2:3078736/1‑93 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGACCCCTCGAGTTCCATAATCTGCaaa                                                          >  2:2957803/1‑93 (MQ=255)

                                gAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAGAt                                                         >  2:4021833/1‑93 (MQ=255)

                                gAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAAt                                                         <  2:2661714/93‑1 (MQ=255)

                                 aaTGTAACTGGTGCTGATAATCGCTACATCATCCTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATg                                                        <  1:163441/93‑1 (MQ=255)

                                  aTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGa                                                       <  2:558112/93‑1 (MQ=255)

                                   tGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGAc                                                      <  2:4308664/93‑1 (MQ=255)

                                    gTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACg                                                     <  2:2623780/93‑1 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    <  1:2270018/93‑1 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    <  1:506693/93‑1 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    <  1:2361503/93‑1 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   >  1:4196210/1‑93 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   >  2:1914996/1‑93 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   <  1:3178352/93‑1 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAACGACGAg                                                   <  2:4669021/93‑1 (MQ=255)

                                       aGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAgg                                                  <  1:934860/93‑1 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 <  2:3605228/93‑1 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 <  2:3458339/93‑1 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 >  1:4298231/1‑93 (MQ=255)

                                           gtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTa                                              <  2:1646732/93‑1 (MQ=255)

                                            tgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTAt                                             <  2:401998/93‑1 (MQ=255)

                                               tGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCgg                                          >  2:1373689/1‑93 (MQ=255)

                                                gATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGt                                         >  2:803888/1‑93 (MQ=255)

                                                 aTAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGtt                                        >  2:3335466/1‑93 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  1:1755224/93‑1 (MQ=255)

                                                  tAATCGCTACAGCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  2:2513344/93‑1 (MQ=255)

                                                   aaTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtat                                      <  2:351698/93‑1 (MQ=255)

                                                   aaTCGCTACATCAACTTTAGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtat                                      <  2:3763883/93‑1 (MQ=255)

                                                    aTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtata                                     >  1:3419620/1‑93 (MQ=255)

                                                     tCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTatat                                    >  2:2487765/1‑93 (MQ=255)

                                                       gCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATct                                  >  2:2425746/1‑92 (MQ=255)

                                                       gCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATcc                                  <  1:1840551/93‑1 (MQ=255)

                                                       ccTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATcc                                  <  2:3430903/92‑1 (MQ=255)

                                                         tACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAg                                >  1:4129824/1‑93 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               <  1:3250791/93‑1 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               >  1:2615263/1‑93 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               >  2:1066208/1‑93 (MQ=255)

                                                            aTCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgcgc                             <  2:2010307/93‑1 (MQ=255)

                                                            aTCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgcgc                             >  2:1469266/1‑93 (MQ=255)

                                                            aTCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgcgc                             >  2:1279003/1‑93 (MQ=255)

                                                             tCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGcgcg                            <  1:1416096/93‑1 (MQ=255)

                                                              cAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCgg                           <  1:579649/93‑1 (MQ=255)

                                                               atcTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          <  1:4404572/91‑1 (MQ=255)

                                                               aaCTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          >  2:478456/1‑93 (MQ=255)

                                                                    tGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtaggta                     <  1:3990046/93‑1 (MQ=255)

                                                                    tGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtaggta                     >  1:422837/1‑93 (MQ=255)

                                                                     gggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaa                    >  1:216249/1‑93 (MQ=255)

                                                                      ggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaat                   >  1:1123746/1‑93 (MQ=255)

                                                                      ggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaat                   <  1:3983007/93‑1 (MQ=255)

                                                                        gtcatcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaataa                 <  1:3294538/92‑1 (MQ=255)

                                                                        atcatcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaataa                 <  2:3529319/93‑1 (MQ=255)

                                                                        atcatcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaataa                 >  2:3877941/1‑93 (MQ=255)

                                                                         tcatcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAc                >  2:2067874/1‑93 (MQ=255)

                                                                           atcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAccc              <  1:2125019/93‑1 (MQ=255)

                                                                           atcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAccc              <  2:411777/93‑1 (MQ=255)

                                                                             ccAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTTCAAATGTCGATGATTATCCGCTATATCCAGCGCGGATGGTAATAACCCtt            >  2:1134319/1‑93 (MQ=255)

                                                                             ccAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtt            >  2:1329443/1‑93 (MQ=255)

                                                                             ccAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtt            <  2:4681138/93‑1 (MQ=255)

                                                                              cAGTACGCGCGGAAGCTGCGACCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGGAGGTAATAACCCttt           >  2:3242355/1‑93 (MQ=255)

                                                                                gTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTc         <  1:4251697/93‑1 (MQ=255)

                                                                                  aCGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTccc       <  2:2256716/93‑1 (MQ=255)

                                                                                     cgcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCctct    >  1:4028440/1‑93 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   <  2:1083743/93‑1 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   <  1:1526248/93‑1 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   >  2:2917217/1‑93 (MQ=255)

                                                                                       cgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATTCAGCGCGGTAGGTAATAACCCTTTTCCCTCttt  >  1:4564867/1‑93 (MQ=255)

                                                                                       cgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCttt  <  2:366159/93‑1 (MQ=255)

                                                                                       cgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCttt  >  1:664869/1‑93 (MQ=255)

GCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTT  >  NC_000913/3836967‑3837146

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: `‑`