BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,507,420	T→C	T18A (ACG→GCG)	b1437 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,507,420	0	T	C	99.3%	419.6 / ‑4.5	137	T18A (ACG→GCG)	b1437	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (55/81); total (56/81)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

Only 100 of 143 total aligned reads displayed.

ACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGA  >  NC_000913/1507330‑1507509

aCGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCa                                                                                         <  2:1502050/93‑1 (MQ=255)

    tGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCAcccc                                                                                     <  1:369428/93‑1 (MQ=255)

       cTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGc                                                                                  >  2:4509871/1‑93 (MQ=255)

         cGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGa                                                                                <  1:4614083/93‑1 (MQ=255)

          gAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGAt                                                                               <  1:1311168/93‑1 (MQ=255)

          gAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGAt                                                                               <  1:1559661/93‑1 (MQ=255)

          gAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGAt                                                                               <  1:771834/93‑1 (MQ=255)

          gAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGAt                                                                               <  1:2705358/93‑1 (MQ=255)

             tctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCTGCGCCACCCCGGCGATgag                                                                            <  2:4576142/93‑1 (MQ=255)

              ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                           <  2:1393055/93‑1 (MQ=255)

              ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                           <  1:4271428/93‑1 (MQ=255)

              ctcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATgaga                                                                           <  1:2178323/93‑1 (MQ=255)

               tcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAt                                                                          <  2:984235/93‑1 (MQ=255)

               tcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAt                                                                          >  2:2593761/1‑93 (MQ=255)

               tcAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAt                                                                          >  1:857753/1‑93 (MQ=255)

                cAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGAtt                                                                         <  2:1181358/93‑1 (MQ=255)

                  gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaa                                                                       >  1:2583148/1‑93 (MQ=255)

                  gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaa                                                                       >  1:3473434/1‑93 (MQ=255)

                  gggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGACGCGCCACCCCTGCGATGAGATTaa                                                                       >  1:1213355/1‑93 (MQ=255)

                   ggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaaa                                                                      <  1:764845/93‑1 (MQ=255)

                   ggCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTaaa                                                                      <  1:3579782/93‑1 (MQ=255)

                      gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                   >  2:4007900/1‑93 (MQ=255)

                      gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                   >  2:3782117/1‑93 (MQ=255)

                      gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                   <  2:134836/93‑1 (MQ=255)

                      gTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGaa                                                                   <  1:4413432/93‑1 (MQ=255)

                       tCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAAc                                                                  >  2:3380994/1‑93 (MQ=255)

                        cGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAAcc                                                                 <  2:1630812/93‑1 (MQ=255)

                         gATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCa                                                                >  2:2058099/1‑93 (MQ=255)

                          aTGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAt                                                               >  1:1058552/1‑93 (MQ=255)

                          aTGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAt                                                               >  1:2782832/1‑93 (MQ=255)

                           tGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCAtt                                                              <  1:168368/93‑1 (MQ=255)

                             tAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGc                                                            <  2:4470490/93‑1 (MQ=255)

                               gCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGt                                                          >  2:2800826/1‑93 (MQ=255)

                                  aaTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTaaa                                                       >  1:953223/1‑93 (MQ=255)

                                   aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                      <  1:1257870/93‑1 (MQ=255)

                                   aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                      <  2:4494036/93‑1 (MQ=255)

                                   aTGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAg                                                      <  2:2116563/93‑1 (MQ=255)

                                    tGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGc                                                     <  1:3489339/93‑1 (MQ=255)

                                    tGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGc                                                     <  1:600180/93‑1 (MQ=255)

                                    tGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGc                                                     >  1:4169934/1‑93 (MQ=255)

                                       cagcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAAt                                                  <  2:4063377/93‑1 (MQ=255)

                                        agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                 <  2:379958/93‑1 (MQ=255)

                                        agcaACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATc                                                 <  2:3383842/93‑1 (MQ=255)

                                          caACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCt                                               <  2:4113814/93‑1 (MQ=255)

                                            aCCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGa                                             >  1:1651310/1‑93 (MQ=255)

                                              cATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaaa                                           <  1:3497822/93‑1 (MQ=255)

                                              cATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGaaa                                           >  2:3945445/1‑93 (MQ=255)

                                               aTTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaac                                          >  1:4277182/1‑93 (MQ=255)

                                               aTTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaac                                          >  2:4004492/1‑93 (MQ=255)

                                                tttGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAaaca                                         >  1:781516/1‑93 (MQ=255)

                                                   gAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACt                                      <  1:3573371/93‑1 (MQ=255)

                                                    aaaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTc                                     <  1:2246675/93‑1 (MQ=255)

                                                     aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                    >  2:315366/1‑93 (MQ=255)

                                                     aaCTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCg                                    >  2:703841/1‑93 (MQ=255)

                                                       cTCAGGTTTCATGGGAGGCGCATTGTCATGCCGCGCCACCACTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                  <  2:1207750/93‑1 (MQ=255)

                                                       cTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGt                                  <  1:2056425/93‑1 (MQ=255)

                                                        tCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGtt                                 <  1:1856072/93‑1 (MQ=255)

                                                        tCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGtt                                 <  2:2950976/93‑1 (MQ=255)

                                                          aGGTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTg                               >  2:2895732/1‑93 (MQ=255)

                                                           ggTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGa                              >  2:3003192/1‑93 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             <  2:4216453/93‑1 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             <  2:3950999/93‑1 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             <  1:4193393/93‑1 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             >  1:2117289/1‑93 (MQ=255)

                                                            gTTTCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAg                             >  1:3690718/1‑93 (MQ=255)

                                                             tttCATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGt                            <  2:3835669/93‑1 (MQ=255)

                                                                cATGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTaa                         <  1:4624915/93‑1 (MQ=255)

                                                                  tGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATc                       <  1:278966/93‑1 (MQ=255)

                                                                  tGGGAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTAAAACAACTCGGTTTGAGTTAATc                       <  1:1711565/93‑1 (MQ=255)

                                                                  tGGGAGGCGCAGTGTCATGCCGCGCCACCCCGGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATc                       <  1:2043887/93‑1 (MQ=255)

                                                                   gggAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCg                      >  1:3070826/1‑93 (MQ=255)

                                                                   gggAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAACCCTGAAACAACTCGGTTTGAGTTAATCg                      <  2:4276919/93‑1 (MQ=255)

                                                                     gAGGCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGcc                    <  2:263839/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  <  1:3821134/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCaa                  <  1:1490155/93‑1 (MQ=255)

                                                                       ggCGCAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGATTTGAGTTAATCGCCaa                  <  2:1874279/93‑1 (MQ=255)

                                                                        gcgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAt                 <  1:518750/93‑1 (MQ=255)

                                                                        gcgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAAt                 <  1:275092/93‑1 (MQ=255)

                                                                         cgcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACATCTCGGTTTGAGTTAATCGCCAAtt                <  1:3999467/93‑1 (MQ=255)

                                                                          gcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTa               <  1:2836411/93‑1 (MQ=255)

                                                                          gcAGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTa               >  2:1359519/1‑93 (MQ=255)

                                                                            aGTGTCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaa             <  1:1843529/93‑1 (MQ=255)

                                                                             gtgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaa            >  2:519027/1‑93 (MQ=255)

                                                                              tgtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTaaaaa           >  2:4398326/1‑93 (MQ=255)

                                                                               gtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAg          <  1:2770873/93‑1 (MQ=255)

                                                                               gtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAg          <  1:4048216/93‑1 (MQ=255)

                                                                               gtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAg          >  2:2814668/1‑93 (MQ=255)

                                                                               gtCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAg          <  1:1310682/93‑1 (MQ=255)

                                                                                tCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAgg         >  2:3163285/1‑93 (MQ=255)

                                                                                tCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAgg         >  1:4125155/1‑93 (MQ=255)

                                                                                tCATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAgg         >  2:768674/1‑93 (MQ=255)

                                                                                 cATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGt        <  2:1347193/93‑1 (MQ=255)

                                                                                 cATGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGt        >  2:976991/1‑93 (MQ=255)

                                                                                  aTGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGtt       >  1:2526320/1‑93 (MQ=255)

                                                                                   tGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTa      >  1:4018533/1‑93 (MQ=255)

                                                                                   tGCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTa      <  2:740716/93‑1 (MQ=255)

                                                                                    gCCGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTaa     >  2:1895367/1‑93 (MQ=255)

                                                                                     ccGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAAt    <  1:2237637/93‑1 (MQ=255)

                                                                                     ccGCGCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAAt    <  1:2117368/93‑1 (MQ=255)

                                                                                      cgcgCCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATg   >  2:4053486/1‑93 (MQ=255)

ACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTCATGCCGCGTCACCCCTGCGATGAGATTAAAGAACCATTGCGTAAAGCAATCCTGAAACAACTCGGTTTGAGTTAATCGCCAATTAAAAAGGTTAATGA  >  NC_000913/1507330‑1507509

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