BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,840,766	G→A	C243C (TGC→TGT)	yicO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,840,766	0	G	A	81.9%	746.0 / 116.7	293	C243C (TGC→TGT)	yicO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (23/30); new base A (122/118); total (145/148)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 297 total aligned reads displayed.

CGCTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTAAC  >  NC_000913/3840675‑3840856

 gcTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCacaaca                                                                                          <  2:5043946/93‑1 (MQ=255)

 gcTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCa                                                                                          >  1:6159394/1‑93 (MQ=255)

 gcTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCa                                                                                          >  1:1549089/1‑93 (MQ=255)

   tCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAg                                                                                        >  2:4916961/1‑93 (MQ=255)

    cAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGa                                                                                       >  2:359796/1‑93 (MQ=255)

     aaaTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGAc                                                                                      <  1:5063945/93‑1 (MQ=255)

     aaaTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGAc                                                                                      >  2:2671977/1‑93 (MQ=255)

       aTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGt                                                                                    <  2:3541067/93‑1 (MQ=255)

       aTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGt                                                                                    >  1:6641170/1‑93 (MQ=255)

       aTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGt                                                                                    <  1:3514906/93‑1 (MQ=255)

        tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                   >  1:7254853/1‑93 (MQ=255)

        tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                   >  2:3820863/1‑93 (MQ=255)

         cTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTca                                                                                  <  2:2523231/93‑1 (MQ=255)

          tACTTGACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                 >  1:3223875/1‑93 (MQ=255)

          tACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                 <  2:3184951/93‑1 (MQ=255)

           aCTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcacc                                                                                >  2:3044247/1‑93 (MQ=255)

            cTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcacca                                                                               >  2:1934805/1‑93 (MQ=255)

             ttCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTcaccac                                                                              >  2:2277265/1‑93 (MQ=255)

             ttCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcaccac                                                                              >  2:3267324/1‑93 (MQ=255)

                 ccAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATa                                                                          <  2:854937/93‑1 (MQ=255)

                  cAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAg                                                                         >  1:690805/1‑93 (MQ=255)

                   aaTGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGa                                                                        <  2:5802544/93‑1 (MQ=255)

                     tGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGaaa                                                                      >  2:6799522/1‑93 (MQ=255)

                     tGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGaaa                                                                      <  1:2875902/93‑1 (MQ=255)

                           cgcTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCa                                                                <  2:1821166/93‑1 (MQ=255)

                            gcTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCAc                                                               <  1:4911924/93‑1 (MQ=255)

                            gcTAATATCAGGCGCAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATATTCCACAACAAGACGTCACCACAATAGAAACCAGCAc                                                               <  2:4461381/93‑1 (MQ=255)

                             cTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCAcc                                                              >  1:1840802/1‑93 (MQ=255)

                              tAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcg                                                             >  1:2812298/1‑93 (MQ=255)

                               aaTATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcgc                                                            <  1:2298464/93‑1 (MQ=255)

                                 tatCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCgg                                                          >  1:2041133/1‑93 (MQ=255)

                                 tatCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCgg                                                          >  2:5626828/1‑93 (MQ=255)

                                   tCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCa                                                        >  1:2374308/1‑93 (MQ=255)

                                    cAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCAt                                                       <  2:4309990/93‑1 (MQ=255)

                                     aGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATg                                                      >  1:145051/1‑93 (MQ=255)

                                      ggCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGa                                                     >  1:5652818/1‑93 (MQ=255)

                                       gCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGaa                                                    <  1:6266985/93‑1 (MQ=255)

                                          gaatggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAAt                                                 >  2:1723924/1‑93 (MQ=255)

                                          gaatggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAAt                                                 >  1:1202595/1‑93 (MQ=255)

                                          gaatggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAAt                                                 >  1:3984936/1‑93 (MQ=255)

                                           aatggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATg                                                <  1:4374953/93‑1 (MQ=255)

                                            atggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGa                                               >  1:5561173/1‑93 (MQ=255)

                                            atggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGa                                               >  1:6837902/1‑93 (MQ=255)

                                            atggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGa                                               <  2:5533162/93‑1 (MQ=255)

                                            atggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGa                                               >  2:7155495/1‑93 (MQ=255)

                                              ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                             <  1:472385/93‑1 (MQ=255)

                                              ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                             >  1:3033233/1‑93 (MQ=255)

                                              ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                             <  2:7229885/93‑1 (MQ=255)

                                               gaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGt                                            <  2:1370723/93‑1 (MQ=255)

                                               gaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGt                                            <  1:4203101/93‑1 (MQ=255)

                                               gaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGt                                            >  2:309647/1‑93 (MQ=255)

                                                 atAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtga                                          <  2:6358263/93‑1 (MQ=255)

                                                 atAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtga                                          <  2:6075167/93‑1 (MQ=255)

                                                  tAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtgat                                         <  2:6926286/93‑1 (MQ=255)

                                                    gACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtgatga                                       <  2:6527683/93‑1 (MQ=255)

                                                      cgccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGaca                                     <  2:1089715/93‑1 (MQ=255)

                                                       gccgcTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaa                                    >  1:4560897/1‑93 (MQ=255)

                                                       gccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaa                                    >  1:3737379/1‑93 (MQ=255)

                                                        ccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaac                                   >  2:5492451/1‑93 (MQ=255)

                                                         cgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaaca                                  >  1:3928418/1‑93 (MQ=255)

                                                          gcTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAacac                                 >  1:3349892/1‑93 (MQ=255)

                                                          gcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAacac                                 <  2:5773269/93‑1 (MQ=255)

                                                          gcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAacac                                 >  2:2172204/1‑93 (MQ=255)

                                                           cTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACg                                >  2:2871111/1‑93 (MQ=255)

                                                           cTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACg                                >  2:1678002/1‑93 (MQ=255)

                                                            tAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACgg                               >  2:1880304/1‑93 (MQ=255)

                                                            tAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACgg                               <  2:4086721/93‑1 (MQ=255)

                                                             aaaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGt                              <  1:1108601/93‑1 (MQ=255)

                                                             aaaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGt                              >  2:6219475/1‑93 (MQ=255)

                                                             aaaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGt                              >  1:7016816/1‑93 (MQ=255)

                                                              aaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGtt                             >  2:1275505/1‑93 (MQ=255)

                                                               aaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTta                            <  2:5932157/93‑1 (MQ=255)

                                                               aaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTta                            <  1:40792/93‑1 (MQ=255)

                                                                aTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTtat                           <  1:956052/93‑1 (MQ=255)

                                                                 tGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTtata                          >  1:2147313/1‑93 (MQ=255)

                                                                    aCATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTataata                       >  2:2106600/1‑93 (MQ=255)

                                                                     cATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTAtaataa                      >  2:4235354/1‑93 (MQ=255)

                                                                      aTCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaa                     <  2:7357928/93‑1 (MQ=255)

                                                                      aTCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaa                     <  1:151379/93‑1 (MQ=255)

                                                                      aTCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaa                     >  2:6090989/1‑93 (MQ=255)

                                                                       tCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaa                    >  1:6330660/1‑93 (MQ=255)

                                                                       tCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaa                    <  2:2973873/93‑1 (MQ=255)

                                                                        cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                   >  2:3418528/1‑93 (MQ=255)

                                                                        cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                   >  2:6046238/1‑93 (MQ=255)

                                                                        cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                   <  1:5317793/93‑1 (MQ=255)

                                                                          ccGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAc                 >  1:2528047/1‑93 (MQ=255)

                                                                           cGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAcc                <  1:6247154/93‑1 (MQ=255)

                                                                              aaaaTAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTa             >  2:4650472/1‑93 (MQ=255)

                                                                               aaaTAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaa            <  2:1448048/93‑1 (MQ=255)

                                                                                aataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaa           <  2:3999717/93‑1 (MQ=255)

                                                                                 ataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaaa          <  2:6310888/93‑1 (MQ=255)

                                                                                 ataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaaa          <  2:472743/93‑1 (MQ=255)

                                                                                 ataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaaa          >  2:7483842/1‑93 (MQ=255)

                                                                                     tCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATAcc      <  1:19539/93‑1 (MQ=255)

                                                                                     tCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATAcc      >  1:752723/1‑93 (MQ=255)

                                                                                     tCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATAcc      <  1:7365713/93‑1 (MQ=255)

                                                                                      ccACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCt     <  2:531450/93‑1 (MQ=255)

                                                                                      ccACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAAAAAAAAACCCTAAAATACCt     <  2:15619/93‑1 (MQ=255)

                                                                                        acaacaAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTaa   >  1:2793439/1‑93 (MQ=255)

                                                                                         caGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTaac  <  2:866023/93‑1 (MQ=255)

CGCTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTAAC  >  NC_000913/3840675‑3840856

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