BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,008,642	C→T	A198A (GCC→GCT)	b2873 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,008,642	0	C	T	98.7%	533.2 / ‑3.8	152	A198A (GCC→GCT)	b2873	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (72/78); total (73/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 158 total aligned reads displayed.

ATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCAC  >  NC_000913/3008552‑3008730

aTTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATcgctgc                                                                                        >  1:4309052/1‑93 (MQ=255)

aTTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATcgctgc                                                                                        >  2:1039338/1‑93 (MQ=255)

 ttACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCa                                                                                       <  1:1641696/93‑1 (MQ=255)

  tACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAg                                                                                      <  1:1286612/93‑1 (MQ=255)

  tACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAg                                                                                      <  2:2968069/93‑1 (MQ=255)

  tACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAg                                                                                      <  2:4476765/93‑1 (MQ=255)

  tACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAg                                                                                      <  2:4395798/93‑1 (MQ=255)

  tACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAg                                                                                      <  2:4306547/93‑1 (MQ=255)

   aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAgg                                                                                     <  1:4289523/93‑1 (MQ=255)

      gcgcTTGCACGAGTCCGGCGACCTGACCAGCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGtt                                                                                  <  1:3956894/92‑1 (MQ=255)

        gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTaa                                                                                >  1:1678604/1‑93 (MQ=255)

        gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTaa                                                                                <  1:3922215/93‑1 (MQ=255)

        gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTaa                                                                                >  1:2455446/1‑93 (MQ=255)

         cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAc                                                                               >  1:4417709/1‑93 (MQ=255)

         cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGAGGCTATAGCCAGCAAGCGGCCGGAATTTATCGCTGCAGCGTTAAc                                                                               >  1:525375/1‑93 (MQ=255)

          ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGTGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAcc                                                                              <  1:3636207/93‑1 (MQ=255)

          ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAcc                                                                              >  1:4534394/1‑93 (MQ=255)

          ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAcc                                                                              <  1:3466427/93‑1 (MQ=255)

            gCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgc                                                                            <  2:2472682/93‑1 (MQ=255)

            gCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgc                                                                            >  2:4529965/1‑93 (MQ=255)

             cATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcg                                                                           >  2:3577100/1‑93 (MQ=255)

             cATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcg                                                                           <  1:4048015/93‑1 (MQ=255)

              aTGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgc                                                                          <  1:1983400/93‑1 (MQ=255)

               tGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgcc                                                                         >  2:2646285/1‑93 (MQ=255)

                gAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgccg                                                                        >  2:4526038/1‑93 (MQ=255)

                 aaTCCGGCTCGCTGAACACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgccgc                                                                       <  1:21910/93‑1 (MQ=255)

                 aaTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCTGGCGGAATTTATCGCTGCAGGGTTAAccgcgccgc                                                                       <  1:3393823/93‑1 (MQ=255)

                 aaTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgccgc                                                                       >  1:2448465/1‑93 (MQ=255)

                 aaTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAAccgcgccgc                                                                       <  1:1520060/93‑1 (MQ=255)

                  aTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGCTGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGGTAAccggcccgcg                                                                      >  1:2740541/1‑93 (MQ=255)

                      ggCGCGCTGACCCCCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTAt                                                                  <  1:3696606/93‑1 (MQ=255)

                      ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTAt                                                                  >  1:2721847/1‑93 (MQ=255)

                      ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTAt                                                                  >  1:3059210/1‑93 (MQ=255)

                      ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTAt                                                                  <  2:370870/93‑1 (MQ=255)

                         gcgcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCAt                                                               >  1:1307215/1‑93 (MQ=255)

                           gcTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGc                                                             >  1:1353338/1‑93 (MQ=255)

                            cTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                            <  2:1888719/93‑1 (MQ=255)

                             tGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCt                                                           <  1:3407119/93‑1 (MQ=255)

                              gACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGTAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCtt                                                          <  1:2900899/93‑1 (MQ=255)

                                 caccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAg                                                       >  2:2031670/1‑93 (MQ=255)

                                  cccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGt                                                      <  1:1948387/92‑1 (MQ=255)

                                   ccGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTc                                                     <  1:1118887/93‑1 (MQ=255)

                                     gTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGc                                                   >  2:315255/1‑93 (MQ=255)

                                       gCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGccc                                                 <  1:1277567/93‑1 (MQ=255)

                                       gCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGccc                                                 <  1:2913629/93‑1 (MQ=255)

                                       gCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGccc                                                 <  1:3814366/93‑1 (MQ=255)

                                       gCACCCGGAAAATGAGGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTACCCGCGCCGCGCTATCATGCCTTGAGTCGccc                                                 <  2:1713093/93‑1 (MQ=255)

                                         aCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctc                                               <  1:2988105/93‑1 (MQ=255)

                                         aCCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctc                                               >  1:4139438/1‑93 (MQ=255)

                                          cccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCctct                                              <  2:1950156/93‑1 (MQ=255)

                                           ccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTg                                             <  2:3620272/93‑1 (MQ=255)

                                           ccGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTg                                             >  1:1548547/1‑93 (MQ=255)

                                            cGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTgg                                            >  1:1425074/1‑93 (MQ=255)

                                             ggAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGa                                           >  2:4423899/1‑93 (MQ=255)

                                             ggAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGa                                           >  1:3930922/1‑93 (MQ=255)

                                             ggAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGa                                           >  2:3277205/1‑93 (MQ=255)

                                              gAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGaa                                          >  2:4271303/1‑93 (MQ=255)

                                              gAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGaa                                          <  1:1020671/93‑1 (MQ=255)

                                               aaaaTGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAAt                                         <  2:3510814/93‑1 (MQ=255)

                                               aaaaTGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAAt                                         <  1:4702154/93‑1 (MQ=255)

                                                aaaTGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATg                                        <  1:31344/93‑1 (MQ=255)

                                                    gatgCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGaa                                    <  1:363172/93‑1 (MQ=255)

                                                    gaggCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGaa                                    <  2:3847395/90‑1 (MQ=255)

                                                         ggCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGa                               >  2:3912105/1‑93 (MQ=255)

                                                         ggCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGa                               >  2:70740/1‑93 (MQ=255)

                                                          gCTATCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGaa                              <  1:3134627/93‑1 (MQ=255)

                                                            tATCGCCAGCAAGCGGGCGGAATTTATTGCTGCAGGGTTAACCGCGCCGCGATATCATGCCTTGAGTCGCCCTCTGGAATGcgaagccgaagc                            >  2:2227638/1‑93 (MQ=255)

                                                             aTCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGcc                           >  1:625022/1‑93 (MQ=255)

                                                              tCGCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCa                          <  2:2456587/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:3092212/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:4025434/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        >  1:2257965/1‑93 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:1109864/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:817683/93‑1 (MQ=255)

                                                                gCCAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATc                        <  2:1955757/93‑1 (MQ=255)

                                                                  cAGCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGc                      >  2:4160266/1‑93 (MQ=255)

                                                                    gCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                    <  2:974417/93‑1 (MQ=255)

                                                                    gCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCATTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                    >  1:4133486/1‑93 (MQ=255)

                                                                    gCAAGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCATTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                    >  1:4133485/1‑93 (MQ=255)

                                                                      aaGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                  >  2:4592306/1‑93 (MQ=255)

                                                                      aaGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                  >  1:3747159/1‑93 (MQ=255)

                                                                      aaGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                  <  2:622946/93‑1 (MQ=255)

                                                                      aaGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                  >  2:1223880/1‑93 (MQ=255)

                                                                       aGCGGGCGTAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 >  2:1517610/1‑93 (MQ=255)

                                                                       aGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 >  1:4437208/1‑93 (MQ=255)

                                                                       aGCGGGCGGAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCa                 >  2:2723899/1‑93 (MQ=255)

                                                                        gcgggcggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCAt                >  2:3808053/1‑93 (MQ=255)

                                                                           ggcggAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGAt             <  2:2215349/93‑1 (MQ=255)

                                                                             cggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTa           <  1:628439/93‑1 (MQ=255)

                                                                              ggAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTaa          >  1:3655336/1‑93 (MQ=255)

                                                                               gAATTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAc         >  2:869854/1‑93 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        >  1:3391642/1‑93 (MQ=255)

                                                                                aaTTTATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcc        <  1:1877258/93‑1 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAcccg      <  2:1455284/93‑3 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTg      >  2:4543072/1‑93 (MQ=255)

                                                                                  tttATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTg      >  2:4543073/1‑93 (MQ=255)

                                                                                    tATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCTGCATGATTAACCTGGc    <  2:2505159/93‑1 (MQ=255)

                                                                                    tATCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGc    >  2:302867/1‑93 (MQ=255)

                                                                                      tCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGcac  <  1:2985774/93‑1 (MQ=255)

                                                                                      tCGCTGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGcac  >  1:2544537/1‑93 (MQ=255)

ATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCAC  >  NC_000913/3008552‑3008730

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