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	99.7%	2229.7 / ‑5.0	608	L188L (CTC→CTT)	nlpA	lipoprotein‑28
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (270/335); total (272/336)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

Only 100 of 629 total aligned reads displayed.

GCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTT  >  NC_000913/3836967‑3837146

gCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAg                                                                                          >  1:2765359/1‑93 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑c                                                                                        >  2:1962593/1‑93 (MQ=255)

       tGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGcc                                                                                  >  2:447569/1‑93 (MQ=255)

       tGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGcc                                                                                  >  1:4627795/1‑93 (MQ=255)

        gCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGccc                                                                                 >  1:5576507/1‑93 (MQ=255)

          aCCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGccccc                                                                               <  1:8307475/93‑1 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCt                                                                              >  2:5796887/1‑93 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCt                                                                              <  1:8756052/93‑1 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGGGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCt                                                                              <  1:9688401/93‑1 (MQ=255)

              gagaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGa                                                                           >  1:10034548/1‑93 (MQ=255)

              gagaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGa                                                                           <  1:7021326/93‑1 (MQ=255)

               agaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAg                                                                          <  2:777665/93‑1 (MQ=255)

                gaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGt                                                                         >  2:4150856/1‑93 (MQ=255)

                gaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGt                                                                         <  1:4442949/93‑1 (MQ=255)

                 aaaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGtt                                                                        <  1:2896344/93‑1 (MQ=255)

                 aaaGCCCGGTCTGCTGAATGTAAGTGATGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGtt                                                                        <  2:1511736/93‑1 (MQ=255)

                  aaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTc                                                                       <  1:9404991/93‑1 (MQ=255)

                   aGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTcc                                                                      <  2:7985517/93‑1 (MQ=255)

                       cGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATaa                                                                  <  1:3518109/93‑1 (MQ=255)

                        ggTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAAt                                                                 <  1:2337277/93‑1 (MQ=255)

                        ggTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAAt                                                                 >  1:5835911/1‑93 (MQ=255)

                         gTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATc                                                                >  2:3276839/1‑93 (MQ=255)

                         gTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATc                                                                <  1:2650896/93‑1 (MQ=255)

                          tCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCt                                                               >  1:486536/1‑93 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGc                                                             >  2:10046462/1‑93 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCaaa                                                          >  2:46470/1‑93 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCaaa                                                          <  2:299724/93‑1 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCaaa                                                          <  1:3452465/93‑1 (MQ=255)

                                gAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAAt                                                         <  1:7271947/93‑1 (MQ=255)

                                 aaTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTGCGCGCGGAAG‑TTGCGCCCCCTCGAGTTCCATAATCTGCAAATg                                                        <  2:7483155/93‑1 (MQ=255)

                                 aaTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGACCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATg                                                        <  2:9457082/93‑1 (MQ=255)

                                  aTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGa                                                       <  2:645891/93‑1 (MQ=255)

                                  aTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGa                                                       <  1:6312506/93‑1 (MQ=255)

                                   tGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGAc                                                      >  1:5313101/1‑93 (MQ=255)

                                   tGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGAc                                                      >  1:7910611/1‑93 (MQ=255)

                                    gTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACg                                                     >  1:3302938/1‑93 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    >  2:3952317/1‑93 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    >  1:9951886/1‑93 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   <  1:2251800/93‑1 (MQ=255)

                                       aGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAgg                                                  >  1:164906/1‑93 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 >  1:7816860/1‑93 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 <  2:3363808/93‑1 (MQ=255)

                                         tggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGAt                                                >  1:5363526/1‑93 (MQ=255)

                                         tggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGAt                                                >  2:2338728/1‑93 (MQ=255)

                                           gtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTa                                              <  1:4331703/93‑1 (MQ=255)

                                            tgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTAt                                             <  2:1714961/93‑1 (MQ=255)

                                            tgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTAt                                             >  1:5193958/1‑93 (MQ=255)

                                            tgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTAt                                             <  2:6915154/93‑1 (MQ=255)

                                            tgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTAt                                             <  1:2108208/93‑1 (MQ=255)

                                             gCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATc                                            >  1:8998053/1‑93 (MQ=255)

                                              cTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCg                                           <  2:9588586/93‑1 (MQ=255)

                                              cTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCg                                           >  2:1339884/1‑93 (MQ=255)

                                               tGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑ATGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATAgg                                          >  1:6433726/1‑93 (MQ=255)

                                                gATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGt                                         >  2:7012299/1‑93 (MQ=255)

                                                gATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGt                                         >  2:5428573/1‑93 (MQ=255)

                                                 aTAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGtt                                        >  1:7615618/1‑93 (MQ=255)

                                                 aTAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGtt                                        <  2:8848321/93‑1 (MQ=255)

                                                  tactCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  1:6833821/90‑1 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  2:9071210/93‑1 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  2:2153212/93‑1 (MQ=255)

                                                   aaTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtat                                      <  2:2594904/93‑1 (MQ=255)

                                                    aTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtata                                     >  1:3917850/1‑93 (MQ=255)

                                                     tCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTatat                                    >  2:3838216/1‑93 (MQ=255)

                                                      cGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATc                                   <  2:6571084/93‑1 (MQ=255)

                                                      cGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATc                                   <  2:5064928/93‑1 (MQ=255)

                                                       gCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATcc                                  >  1:4094783/1‑93 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCa                                 >  2:8092121/1‑93 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCa                                 >  2:2471256/1‑93 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGACGATTATCGGTCATATTCa                                 >  2:3549721/1‑93 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               >  2:6696767/1‑93 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               <  2:8255904/93‑1 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               <  1:5639713/93‑1 (MQ=255)

                                                            aTCAACTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTACCGGTTATATCCAgcgc                             <  2:8472068/93‑1 (MQ=255)

                                                               aaCTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          >  1:1962241/1‑93 (MQ=255)

                                                               aaCTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          <  1:9568843/93‑1 (MQ=255)

                                                                aCTTTGGGATCATCCAGTACGCGCGGAAGCCTGC‑CCCCCGCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggta                         <  1:7984455/93‑1 (MQ=255)

                                                                aCTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggta                         >  2:7256704/1‑93 (MQ=255)

                                                                 cTTTGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtag                        <  2:2803418/93‑1 (MQ=255)

                                                                   ttGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtaggt                      >  2:874140/1‑93 (MQ=255)

                                                                    tGGGATCATCCAGTACGCGTGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtaggta                     <  1:3105930/93‑1 (MQ=255)

                                                                    tGGGATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtaggta                     <  1:60/93‑1 (MQ=255)

                                                                     gggATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaa                    <  1:1793412/93‑1 (MQ=255)

                                                                     gggATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaa                    >  2:4873705/1‑93 (MQ=255)

                                                                      ggATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaat                   <  1:6647927/93‑1 (MQ=255)

                                                                       gATCATCCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaata                  >  2:9278778/1‑93 (MQ=255)

                                                                        atcatcCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaataa                 <  1:3710932/93‑1 (MQ=255)

                                                                         tcatcCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAc                <  2:6269391/93‑1 (MQ=255)

                                                                         tcatcCAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAc                >  1:2532567/1‑93 (MQ=255)

                                                                             ccAGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtt            <  1:1227153/93‑1 (MQ=255)

                                                                               aGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtttt          <  1:8438112/93‑1 (MQ=255)

                                                                               aGTACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtttt          <  1:6562207/93‑1 (MQ=255)

                                                                                 tACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        >  2:6133590/1‑93 (MQ=255)

                                                                                 tACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        >  2:8586842/1‑93 (MQ=255)

                                                                                 tACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        >  1:8144454/1‑93 (MQ=255)

                                                                                 tACGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        >  1:6135831/1‑93 (MQ=255)

                                                                                  aCGCGCGGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTccc       <  2:4553570/93‑1 (MQ=255)

                                                                                   cgcgcgGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCct      <  1:9107869/93‑1 (MQ=255)

                                                                                    gcgcgGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCctt     >  1:3035315/1‑92 (MQ=255)

                                                                                      gcgGGAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATAGCCAGCGCGGTAGGTAATAACCCTTTTCcttctt   >  1:1429926/1‑93 (MQ=255)

                                                                                      gcgGAAG‑CTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   >  1:7448298/1‑93 (MQ=255)

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