BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,062,187	A→G	L387S (TTG→TCG)	b3876 ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,062,187	0	A	G	88.9%	823.2 / 86.6	315	L387S (TTG→TCG)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base A (18/17); new base G (131/149); total (149/166)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

Only 100 of 323 total aligned reads displayed.

ATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCATTCACCG  >  NC_000913/4062097‑4062279

aTAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCgag                                                                                            >  1:3684678/1‑93 (MQ=255)

   aaaaTGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCgagaag                                                                                         <  2:6353439/93‑1 (MQ=255)

      aTGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCa                                                                                      >  1:5123379/1‑93 (MQ=255)

      aTGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCa                                                                                      <  1:1395244/93‑1 (MQ=255)

          gCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCaa                                                                                  >  1:7509669/1‑93 (MQ=255)

          gCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCaa                                                                                  <  1:528635/93‑1 (MQ=255)

          gCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCaa                                                                                  <  1:2937057/93‑1 (MQ=255)

          gCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCaa                                                                                  <  2:6518041/93‑1 (MQ=255)

           cTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAg                                                                                 <  1:6502034/93‑1 (MQ=255)

            tGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGt                                                                                <  2:22734/93‑1 (MQ=255)

             gACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTg                                                                               <  2:890229/93‑1 (MQ=255)

             gACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTg                                                                               >  2:21694/1‑93 (MQ=255)

              aCGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGc                                                                              <  1:2877259/93‑1 (MQ=255)

                gCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCgg                                                                            >  2:4630353/1‑93 (MQ=255)

                gCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCgg                                                                            >  1:4882120/1‑93 (MQ=255)

                   aaCCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAg                                                                         >  2:2556928/1‑93 (MQ=255)

                      ccTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGaaa                                                                      <  2:1706984/93‑1 (MQ=255)

                      ccTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGaaa                                                                      <  1:5353461/93‑1 (MQ=255)

                       cTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAt                                                                     <  1:5703249/93‑1 (MQ=255)

                       cTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAt                                                                     >  2:1053654/1‑93 (MQ=255)

                       cTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAt                                                                     >  2:2974600/1‑93 (MQ=255)

                       cTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAt                                                                     >  2:1972709/1‑93 (MQ=255)

                        tGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAtt                                                                    <  1:84636/93‑1 (MQ=255)

                        tGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAtt                                                                    >  2:7203423/1‑93 (MQ=255)

                        tGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAtt                                                                    <  2:1441537/93‑1 (MQ=255)

                        tGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAAtt                                                                    >  2:558234/1‑93 (MQ=255)

                         gTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAttt                                                                   <  2:4175622/93‑1 (MQ=255)

                         gTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAttt                                                                   >  1:125217/1‑93 (MQ=255)

                         gTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAAttt                                                                   >  1:1898916/1‑93 (MQ=255)

                         gTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAAttt                                                                   >  1:6965907/1‑93 (MQ=255)

                           aaTGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTa                                                                 <  1:1707142/93‑1 (MQ=255)

                           aaTGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTa                                                                 >  1:7044036/1‑93 (MQ=255)

                             tGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACg                                                               >  1:733107/1‑93 (MQ=255)

                             tGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACg                                                               >  2:7346785/1‑93 (MQ=255)

                                 gCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaag                                                           >  2:3141747/1‑93 (MQ=255)

                                  cATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaaga                                                          <  2:4996303/93‑1 (MQ=255)

                                  cATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaaga                                                          >  1:6859585/1‑93 (MQ=255)

                                   aTCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaagaa                                                         >  1:3886594/1‑93 (MQ=255)

                                   aTCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaagaa                                                         <  1:2007100/93‑1 (MQ=255)

                                    tCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaagaag                                                        <  1:1049987/93‑1 (MQ=255)

                                    tCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAaagaag                                                        <  2:3365032/93‑1 (MQ=255)

                                    tCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCAAGCAAGTGCGCCAGAAATTTTACGAaaaaag                                                        >  1:7104251/1‑93 (MQ=255)

                                      gCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGt                                                      >  2:4429567/1‑93 (MQ=255)

                                        tttGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGtat                                                    >  2:5874811/1‑93 (MQ=255)

                                        tttGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGtat                                                    >  2:3803564/1‑93 (MQ=255)

                                        tttGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGtat                                                    <  1:5419316/93‑1 (MQ=255)

                                         ttGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGtata                                                   <  1:5560801/93‑1 (MQ=255)

                                           gcgcGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAAc                                                 <  2:5199560/93‑1 (MQ=255)

                                           gcgcGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAAc                                                 >  2:5769605/1‑93 (MQ=255)

                                             gcgATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCg                                               >  2:1516663/1‑93 (MQ=255)

                                               gATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGt                                             >  2:47763/1‑93 (MQ=255)

                                               gATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGt                                             >  1:4381641/1‑93 (MQ=255)

                                                aTGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGta                                            <  1:3380452/93‑1 (MQ=255)

                                                 tGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGtat                                           <  2:2796062/93‑1 (MQ=255)

                                                 tGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGtat                                           >  1:6241683/1‑93 (MQ=255)

                                                 tGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGtat                                           >  1:4519564/1‑93 (MQ=255)

                                                 tGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGtat                                           >  2:834887/1‑93 (MQ=255)

                                                   ttGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATaa                                         >  2:1319314/1‑93 (MQ=255)

                                                    tGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATaaa                                        <  1:4326224/93‑1 (MQ=255)

                                                      ggAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACa                                      <  2:61878/93‑1 (MQ=255)

                                                      ggAACATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACa                                      <  1:3784637/93‑1 (MQ=255)

                                                        aaCATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGa                                    <  1:3051958/93‑1 (MQ=255)

                                                         aCATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGaa                                   <  1:4470754/93‑1 (MQ=255)

                                                         aCATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGaa                                   <  1:3110905/93‑1 (MQ=255)

                                                         aCATAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGaa                                   >  1:6431093/1‑93 (MQ=255)

                                                         aCATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGaa                                   <  2:6937992/93‑1 (MQ=255)

                                                            tAACCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCt                                >  1:3450226/1‑93 (MQ=255)

                                                             aaCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCtt                               >  2:2714629/1‑93 (MQ=255)

                                                             aaCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCtt                               <  1:4249230/93‑1 (MQ=255)

                                                              aCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTc                              <  2:3695916/93‑1 (MQ=255)

                                                              aCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTc                              <  1:1545878/93‑1 (MQ=255)

                                                              aCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTc                              <  1:262914/93‑1 (MQ=255)

                                                              aCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTc                              <  1:6548898/93‑1 (MQ=255)

                                                              aCCAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTc                              <  1:6027960/93‑1 (MQ=255)

                                                               ccAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCg                             <  1:6066430/93‑1 (MQ=255)

                                                               ccAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCg                             <  1:7445736/93‑1 (MQ=255)

                                                                cAATTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCgg                            >  2:5467870/1‑93 (MQ=255)

                                                                  aTTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCa                          >  2:4832308/1‑93 (MQ=255)

                                                                  aTTTGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCa                          <  1:6005374/93‑1 (MQ=255)

                                                                   tttGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCAc                         >  1:3276121/1‑93 (MQ=255)

                                                                    ttGCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACg                        >  1:2140587/1‑93 (MQ=255)

                                                                      gCGTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGaa                      <  2:2454896/93‑1 (MQ=255)

                                                                        gTCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATa                    <  1:2555054/93‑1 (MQ=255)

                                                                         tCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATAc                   <  2:5128438/93‑1 (MQ=255)

                                                                         tCAGCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATAc                   >  2:7459994/1‑93 (MQ=255)

                                                                            gCATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCgg                <  2:2098004/93‑1 (MQ=255)

                                                                            gCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCgg                <  2:7130062/93‑1 (MQ=255)

                                                                             cATAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGt               <  1:2184941/93‑1 (MQ=255)

                                                                              aTAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGtt              >  2:5960267/1‑93 (MQ=255)

                                                                               tAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGttt             >  2:6143037/1‑93 (MQ=255)

                                                                               tAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGttt             >  1:6910338/1‑93 (MQ=255)

                                                                               tAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGttt             >  2:672189/1‑93 (MQ=255)

                                                                               tAATGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGttt             >  1:3988768/1‑93 (MQ=255)

                                                                                aaTGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGtttt            <  2:3901778/93‑1 (MQ=255)

                                                                                  tGCCTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTcc          <  2:3258578/93‑1 (MQ=255)

                                                                                    ccTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCAt        >  1:3591483/1‑93 (MQ=255)

                                                                                     cTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCAtt       <  2:3432029/93‑1 (MQ=255)

                                                                                     cTGGCGAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCAtt       >  1:4959714/1‑93 (MQ=255)

                                                                                      tGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCATTc      <  2:4256740/93‑1 (MQ=255)

                                                                                       ggCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCATTCa     <  2:7019158/93‑1 (MQ=255)

ATAAAAATGAGCTGACGCAAACCCTGTAATGTCGCATCGCTTTGCGCGATGTTGGGAACATAACCAATTTGCGTCAGCATAATGCCTGGCAAGAAGCCAGCAAGTGCGGCAGAAATTTTACGAAAGAAGGTATAACCGGTATAAACAGAACCTTCGGCACGAATACCGGTTTTCCATTCACCG  >  NC_000913/4062097‑4062279

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