BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,089,241	G→A	R30H (CGC→CAC)	b2946 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,089,241	0	G	A	98.9%	1260.2 / ‑2.5	356	R30H (CGC→CAC)	b2946	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (3/1); new base A (187/165); total (190/166)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 368 total aligned reads displayed.

CTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAA  >  NC_000913/3089151‑3089330

 tATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGt                                                                                        <  2:3611107/93‑1 (MQ=255)

 tATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGt                                                                                        >  2:3858089/1‑93 (MQ=255)

   tGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCCACCATATCGGGCACGTAc                                                                                      >  2:4197341/1‑93 (MQ=255)

   tGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTAc                                                                                      <  2:1728558/93‑1 (MQ=255)

   tGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTAc                                                                                      <  1:2426692/93‑1 (MQ=255)

   tGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTAc                                                                                      <  1:1903291/93‑1 (MQ=255)

   tGCGTAGCCCCCGCATTTATCATCCTGAACCACTGTCCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTAc                                                                                      <  1:306208/93‑1 (MQ=255)

      gTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTgc                                                                                   >  1:4057747/1‑93 (MQ=255)

      gTATCCCCCGCATGTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTgc                                                                                   <  1:1134428/93‑1 (MQ=255)

         tCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCa                                                                                >  2:4975699/1‑93 (MQ=255)

         tCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCa                                                                                <  2:2250249/93‑1 (MQ=255)

           ccccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATg                                                                              <  1:1038711/93‑1 (MQ=255)

           ccccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATg                                                                              <  1:3566279/93‑1 (MQ=255)

            cccGCATTTATCATCCTGACCCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATgg                                                                             <  1:1440335/93‑1 (MQ=255)

            cccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATgg                                                                             <  2:5111474/93‑1 (MQ=255)

            cccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATgg                                                                             <  1:2455614/93‑1 (MQ=255)

             ccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATggg                                                                            <  1:1869610/93‑1 (MQ=255)

             ccGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATggg                                                                            >  1:2681270/1‑93 (MQ=255)

             ccGCATTTATCATCCTGAACCACTGACCAGCCATTCCCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATggg                                                                            <  2:4785991/93‑1 (MQ=255)

              cGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATgggg                                                                           <  1:7549095/93‑1 (MQ=255)

              cGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATgggg                                                                           >  1:6510288/1‑93 (MQ=255)

               ccATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGc                                                                          <  1:2048292/92‑1 (MQ=255)

                cATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGcc                                                                         >  1:7928501/1‑93 (MQ=255)

                 aTTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCg                                                                        <  1:2108915/93‑1 (MQ=255)

                  tttATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgg                                                                       <  1:3491051/93‑1 (MQ=255)

                  tttATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgg                                                                       >  1:7316208/1‑93 (MQ=255)

                  tttATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgg                                                                       >  2:2468870/1‑93 (MQ=255)

                  tttATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgg                                                                       >  2:2101262/1‑93 (MQ=255)

                   ttATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCggc                                                                      >  1:5200555/1‑92 (MQ=255)

                    tATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgggg                                                                     <  1:4975829/93‑1 (MQ=255)

                    tATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCgggg                                                                     <  1:5572019/93‑1 (MQ=255)

                      tcatcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCa                                                                   >  1:2542145/1‑93 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGTGCATGGGGCCGGGGCAGGc                                                                >  2:7241911/1‑93 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGc                                                                >  2:4880372/1‑93 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGc                                                                >  2:7599932/1‑93 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGc                                                                >  2:3026051/1‑93 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGc                                                                <  2:5533634/93‑1 (MQ=255)

                         tcCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGc                                                                <  1:494308/93‑1 (MQ=255)

                          ccTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGTGCCTGGGCAGGct                                                               >  2:4156518/1‑92 (MQ=255)

                           cTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGt                                                              <  1:1332469/93‑1 (MQ=255)

                           cTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGt                                                              >  1:1182960/1‑93 (MQ=255)

                           cTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGt                                                              >  2:6014544/1‑93 (MQ=255)

                            tGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGtt                                                             >  2:4547483/1‑93 (MQ=255)

                            tGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGtt                                                             >  1:7101084/1‑93 (MQ=255)

                             gAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTg                                                            >  2:1909970/1‑93 (MQ=255)

                             gAACCACGGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTg                                                            <  1:2740551/93‑1 (MQ=255)

                              aaCCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGc                                                           <  2:6619501/93‑1 (MQ=255)

                               aCCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCa                                                          >  2:3282676/1‑93 (MQ=255)

                               aCCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCa                                                          >  1:5740527/1‑93 (MQ=255)

                               aCCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGGTGCa                                                          >  2:2860193/1‑93 (MQ=255)

                                ccACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGGGTTGCat                                                         >  2:5803772/1‑92 (MQ=255)

                                ccACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCaa                                                         >  2:1409763/1‑93 (MQ=255)

                                 cACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAAt                                                        <  1:4400437/93‑1 (MQ=255)

                                  aCTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAAtt                                                       >  2:1810660/1‑93 (MQ=255)

                                   cTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTg                                                      <  1:4283125/93‑1 (MQ=255)

                                    tGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGt                                                     >  2:6736127/1‑93 (MQ=255)

                                    tGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGt                                                     <  2:7858208/93‑1 (MQ=255)

                                     gACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGtt                                                    >  2:6973004/1‑93 (MQ=255)

                                       ccAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGTGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTg                                                  >  1:3196207/1‑93 (MQ=255)

                                       ccAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTg                                                  >  2:1394705/1‑93 (MQ=255)

                                         aGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGAc                                                <  1:1142049/93‑1 (MQ=255)

                                          gCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACg                                               <  2:3509125/93‑1 (MQ=255)

                                          gCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACg                                               >  1:5869513/1‑93 (MQ=255)

                                           ccATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACgg                                              >  2:2877649/1‑93 (MQ=255)

                                            cATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGt                                             <  2:7184591/93‑1 (MQ=255)

                                              ttCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAg                                           <  1:176552/93‑1 (MQ=255)

                                              ttCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAg                                           <  2:378676/93‑1 (MQ=255)

                                                  cacaTCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAAc                                       <  1:6882861/93‑1 (MQ=255)

                                                   acaTCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAAcc                                      <  1:7157722/93‑1 (MQ=255)

                                                     aTCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAg                                    <  1:7686220/93‑1 (MQ=255)

                                                     aTCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAg                                    <  2:6437408/93‑1 (MQ=255)

                                                      tCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAgg                                   >  1:456643/1‑93 (MQ=255)

                                                      tCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGAGGCCGGGGCAGGCGTTGCAATTGTTTGACGGCAGCAACCAgg                                   >  2:3702754/1‑93 (MQ=255)

                                                        gcgcTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTc                                 >  2:3949847/1‑93 (MQ=255)

                                                        gcgcTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTc                                 >  1:4184373/1‑93 (MQ=255)

                                                         cgcTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCt                                >  1:6234567/1‑93 (MQ=255)

                                                         cgcTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCt                                <  2:1066723/93‑1 (MQ=255)

                                                          gcTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCtt                               <  2:136248/93‑1 (MQ=255)

                                                          gcTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCtt                               <  1:5469017/93‑1 (MQ=255)

                                                           cTTTGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCttt                              <  1:4668742/93‑1 (MQ=255)

                                                            tttGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTg                             >  2:5306704/1‑93 (MQ=255)

                                                            tttGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTg                             >  1:3279666/1‑93 (MQ=255)

                                                             ttGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGTTAGCAACCAGGTCTTTGa                            <  2:7744606/93‑1 (MQ=255)

                                                             ttGCGAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGa                            >  2:1245126/1‑93 (MQ=255)

                                                                 gAAGATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGcc                        >  1:3689380/1‑93 (MQ=255)

                                                                    gATGCCGCCAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGaa                     >  1:2220321/1‑93 (MQ=255)

                                                                        ccgccAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTa                 <  2:341780/93‑1 (MQ=255)

                                                                         cgccAACCATATCGGGCACGTATTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTAc                >  2:4460409/1‑93 (MQ=255)

                                                                          gccAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTAcc               <  1:5830383/93‑1 (MQ=255)

                                                                          gccAACCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTAcc               <  1:6642261/93‑1 (MQ=255)

                                                                             aaCCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAgc            >  2:5369362/1‑93 (MQ=255)

                                                                              aCCATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAgcg           <  2:7505649/93‑1 (MQ=255)

                                                                               ccATATCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAgcgc          >  2:5625088/1‑93 (MQ=255)

                                                                                 atatCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCa        <  2:439657/93‑1 (MQ=255)

                                                                                 atatCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCa        <  2:7679148/93‑1 (MQ=255)

                                                                                 atatCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCa        >  2:7020091/1‑93 (MQ=255)

                                                                                   atCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGc      <  1:3128655/93‑1 (MQ=255)

                                                                                   atCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGc      >  1:6555354/1‑93 (MQ=255)

                                                                                    tCGGGCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAATCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCa     >  1:4832915/1‑93 (MQ=255)

                                                                                      gggCACGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCACCaaa   >  2:1502048/1‑93 (MQ=255)

CTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACCATATCGGGCGCGTACTGCGCATGGGGCCGGGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAA  >  NC_000913/3089151‑3089330

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