BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,584,676	A→G	S35G (AGC→GGC)	yhhA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,584,676	0	A	G	85.4%	637.4 / 92.3	247	S35G (AGC→GGC)	yhhA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (21/15); new base G (105/105); total (127/120)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 269 total aligned reads displayed.

TTCTGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCA  >  NC_000913/3584584‑3584767

   tGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACccggcc                                                                                          <  1:4289643/93‑1 (MQ=255)

   tGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACACTAACCAGCCGGGGTATCAGATACccggcc                                                                                          >  2:4350515/1‑93 (MQ=255)

    gATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACccggcca                                                                                         >  1:2516295/1‑93 (MQ=255)

     aTTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACccagccag                                                                                        <  1:5558181/93‑1 (MQ=255)

      ttCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCcggccagc                                                                                       <  1:4455970/93‑1 (MQ=255)

      ttCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCcggccagc                                                                                       >  1:2052304/1‑93 (MQ=255)

         ttACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCagcagc                                                                                    >  1:3185201/1‑93 (MQ=255)

          tACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCg                                                                                   <  2:4185427/93‑1 (MQ=255)

            cGGCCCTCCGGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGa                                                                                 <  1:3077270/93‑1 (MQ=255)

            cGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGa                                                                                 <  2:1623378/93‑1 (MQ=255)

             ggCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGAt                                                                                >  2:4657577/1‑93 (MQ=255)

              gcccTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATg                                                                               <  2:1992968/90‑1 (MQ=255)

              gCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATGCCCGGCCAGCAGCGGATg                                                                               <  2:81814/93‑1 (MQ=255)

              gCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATg                                                                               >  2:4571187/1‑93 (MQ=255)

                ccTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATCCTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCa                                                                             <  2:477954/92‑1 (MQ=255)

                 cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCaa                                                                            >  2:3250066/1‑93 (MQ=255)

                 cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCaa                                                                            >  1:4953624/1‑93 (MQ=255)

                 cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCaa                                                                            >  1:4825385/1‑93 (MQ=255)

                 cTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCaa                                                                            >  1:3531358/1‑93 (MQ=255)

                  tCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCaaa                                                                           >  2:4397007/1‑93 (MQ=255)

                  tCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCaaa                                                                           >  2:2422602/1‑93 (MQ=255)

                   ccTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAAAAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAAc                                                                          <  1:3397680/93‑1 (MQ=255)

                      gCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCa                                                                       >  1:124200/1‑93 (MQ=255)

                       ccGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAg                                                                      <  1:175928/93‑1 (MQ=255)

                       ccGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGCATCAGATACCCGGCCAGCAGCGGATGCAAACCCAg                                                                      >  2:477301/1‑93 (MQ=255)

                        cGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGa                                                                     <  1:3737925/93‑1 (MQ=255)

                        cGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGa                                                                     <  1:4095430/93‑1 (MQ=255)

                        cGTTTGTCGGCTTTGCACAGCCAATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGa                                                                     <  1:1427263/93‑1 (MQ=255)

                         gTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGAt                                                                    <  1:1069918/93‑1 (MQ=255)

                          tttGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATg                                                                   >  1:4680359/1‑93 (MQ=255)

                          tttGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATg                                                                   <  2:5030584/93‑1 (MQ=255)

                          tttGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATg                                                                   >  1:1918318/1‑93 (MQ=255)

                            tGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCa                                                                 >  1:6180402/1‑93 (MQ=255)

                             gTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAg                                                                <  2:2344163/93‑1 (MQ=255)

                              tCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGCATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGa                                                               >  1:3677815/1‑93 (MQ=255)

                                ggCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACt                                                             >  2:5621332/1‑93 (MQ=255)

                                ggCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACt                                                             >  2:770529/1‑93 (MQ=255)

                                ggCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACACAGATGCAGAcc                                                             >  2:3944017/1‑92 (MQ=255)

                                 gCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTc                                                            <  2:5485416/93‑1 (MQ=255)

                                  cTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCa                                                           <  1:1011683/93‑1 (MQ=255)

                                  cTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCa                                                           <  1:2301722/93‑1 (MQ=255)

                                  cTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCa                                                           <  2:5817153/93‑1 (MQ=255)

                                   tttGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAg                                                          <  2:1555767/93‑1 (MQ=255)

                                    ttGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGc                                                         <  1:638428/93‑1 (MQ=255)

                                    ttGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGc                                                         >  2:3431279/1‑93 (MQ=255)

                                     tGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCa                                                        <  1:4854494/93‑1 (MQ=255)

                                      gCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCaa                                                       >  2:3462578/1‑93 (MQ=255)

                                        acaGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAAt                                                     >  1:2139554/1‑93 (MQ=255)

                                         caGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATc                                                    <  2:3709851/93‑1 (MQ=255)

                                         caGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATc                                                    >  1:2434839/1‑93 (MQ=255)

                                         caGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATc                                                    <  1:5017921/93‑1 (MQ=255)

                                          aGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATcc                                                   >  1:1826934/1‑93 (MQ=255)

                                            cccATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAg                                                 >  2:3083420/1‑93 (MQ=255)

                                             ccATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGc                                                >  2:4492678/1‑93 (MQ=255)

                                              cATTAATCCTCTGAACAACCCTAACCAGCCGGGGTATCAGAAACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCa                                               <  2:6548836/93‑1 (MQ=255)

                                              cATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCa                                               >  1:3260797/1‑93 (MQ=255)

                                                 tAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaa                                            >  1:1985230/1‑93 (MQ=255)

                                                  aaTACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCaaaaa                                           >  2:4470715/1‑93 (MQ=255)

                                                   aTACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAg                                          <  2:138024/93‑1 (MQ=255)

                                                    tACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAgg                                         >  2:3810591/1‑93 (MQ=255)

                                                      ctctGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAGCCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                       <  1:2631560/93‑1 (MQ=255)

                                                      ctctGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                       >  1:2939694/1‑93 (MQ=255)

                                                      ctctGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                       <  1:523950/93‑1 (MQ=255)

                                                      ctctGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                       <  2:6506947/93‑1 (MQ=255)

                                                      ctctGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGa                                       >  1:3303038/1‑93 (MQ=255)

                                                         tGAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGc                                    <  1:6241594/93‑1 (MQ=255)

                                                         tGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGc                                    >  1:3286168/1‑93 (MQ=255)

                                                          gAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCt                                   <  1:4526012/93‑1 (MQ=255)

                                                          gAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCt                                   <  1:5423695/93‑1 (MQ=255)

                                                          gAACAACCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCt                                   >  2:1666006/1‑93 (MQ=255)

                                                           aacaacCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTg                                  >  2:275569/1‑93 (MQ=255)

                                                           aacaacCCTAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTg                                  <  1:157390/93‑1 (MQ=255)

                                                                   tAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCa                          <  2:5147401/93‑1 (MQ=255)

                                                                   tAACCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCa                          <  1:4785075/93‑1 (MQ=255)

                                                                     aCCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                        <  1:3831623/93‑1 (MQ=255)

                                                                     aCCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                        >  2:5386481/1‑93 (MQ=255)

                                                                     aCCAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAAc                        >  1:882115/1‑93 (MQ=255)

                                                                      ccAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACt                       <  2:4070701/93‑1 (MQ=255)

                                                                       cAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTg                      >  2:1870884/1‑93 (MQ=255)

                                                                       cAGCCGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTg                      >  2:2508392/1‑93 (MQ=255)

                                                                          ccGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaa                   >  2:1194722/1‑93 (MQ=255)

                                                                           cGGGGTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGaaaa                  <  2:5447217/93‑1 (MQ=255)

                                                                             gggTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACg                >  1:6297746/1‑93 (MQ=255)

                                                                              ggTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGc               >  2:3227109/1‑93 (MQ=255)

                                                                               gTATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCa              >  1:580021/1‑93 (MQ=255)

                                                                                tATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaa             >  2:607757/1‑93 (MQ=255)

                                                                                tATCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaa             >  1:5962293/1‑93 (MQ=255)

                                                                                 aTCAGATACCCGGCGAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACTCaaa            <  1:642799/93‑1 (MQ=255)

                                                                                 aTCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCaaa            <  1:6371103/93‑1 (MQ=255)

                                                                                  tCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAAc           >  1:5475440/1‑93 (MQ=255)

                                                                                  tCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAAc           <  2:6023067/93‑1 (MQ=255)

                                                                                  tCAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAAc           >  1:96999/1‑93 (MQ=255)

                                                                                   cAGATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACt          <  2:2019919/93‑1 (MQ=255)

                                                                                     gATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCa        >  1:5819280/1‑93 (MQ=255)

                                                                                     gATACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCa        >  2:6259159/1‑93 (MQ=255)

                                                                                      aTACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAg       <  1:3777548/93‑1 (MQ=255)

                                                                                       tACCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGt      <  2:3877183/93‑1 (MQ=255)

                                                                                        aCCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGtt     <  2:5781041/93‑1 (MQ=255)

                                                                                        aCCCGGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGtt     <  2:6378733/93‑1 (MQ=255)

                                                                                        aCCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGtt     <  2:2835097/93‑1 (MQ=255)

TTCTGATTCTTACGGCACTCCTGCCGTTTGTCGGCTTTGCACAGCCCATTAATACTCTGAACAACCCTAACCAGCCGGGGTATCAGATACCCAGCCAGCAGCGGATGCAAACCCAGATGCAGACTCAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCA  >  NC_000913/3584584‑3584767

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