BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,322,882	C→T	T146T (ACG→ACA)	phnB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,322,882	0	C	T	99.7%	1331.6 / ‑5.6	370	T146T (ACG→ACA)	phnB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (170/199); total (171/199)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 413 total aligned reads displayed.

CGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTG  >  NC_000913/4322792‑4322971

cGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTgtgt                                                                                          <  1:5235913/93‑1 (MQ=255)

cGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTgtgt                                                                                          <  1:9413776/93‑1 (MQ=255)

 gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgtt                                                                                         <  1:7919797/93‑1 (MQ=255)

 gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgtt                                                                                         <  2:3341326/93‑1 (MQ=255)

   tttGGTTAATTAAAAGTGACTAACCGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttgg                                                                                       >  2:8746767/1‑93 (MQ=255)

   tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttgg                                                                                       >  1:7661189/1‑93 (MQ=255)

   tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttgg                                                                                       >  2:2449573/1‑93 (MQ=255)

   tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttgg                                                                                       >  1:6089661/1‑93 (MQ=255)

     tGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttggtt                                                                                     <  2:9374756/93‑1 (MQ=255)

     tGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTgttggtt                                                                                     >  1:48925/1‑93 (MQ=255)

      ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                    >  1:4381357/1‑93 (MQ=255)

       gTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGggtgt                                                                                   >  2:518229/1‑93 (MQ=255)

        ttaattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGgttgtt                                                                                  >  2:1765065/1‑93 (MQ=255)

            ttaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGttt                                                                              >  1:9282455/1‑93 (MQ=255)

            ttaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGttt                                                                              >  2:5336994/1‑93 (MQ=255)

            ttaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGttt                                                                              >  2:5640262/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCGTCCCGGCGGGGTGTTTTGTTGGTTGGTTTTTga                                                                            >  1:2551804/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTga                                                                            >  2:3598007/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTga                                                                            >  1:9760384/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTga                                                                            <  1:2749407/93‑1 (MQ=255)

               aaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTgac                                                                           <  2:5991551/93‑1 (MQ=255)

                aaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTTTGTTGGGTGTTGttttttg                                                                          >  2:3688297/1‑89 (MQ=255)

                  gTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTgacgac                                                                        <  1:4074566/93‑1 (MQ=255)

                   tGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACa                                                                       <  1:1778472/93‑1 (MQ=255)

                   tGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACa                                                                       >  1:6743430/1‑93 (MQ=255)

                      cTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTa                                                                    >  2:4673552/1‑93 (MQ=255)

                       tAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTaa                                                                   <  2:3980760/93‑1 (MQ=255)

                       tAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTACCGGCGGGTTATTGTGTTGGTTGTTGTATGACGACATTaa                                                                   >  1:6032100/1‑93 (MQ=255)

                        aaCAGATGAAGAGTTAACGTGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGCTGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                  <  2:4981884/93‑1 (MQ=255)

                        aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                  <  1:299342/93‑1 (MQ=255)

                        aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                  >  1:8909283/1‑93 (MQ=255)

                        aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                  >  2:659246/1‑93 (MQ=255)

                         aCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatc                                                                 >  1:8053339/1‑93 (MQ=255)

                         aCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatc                                                                 <  2:8552850/93‑1 (MQ=255)

                         aCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTACCGGCGGGCTAATGAGTTGGGTGTTGATTGGCGACATTAatt                                                                 >  2:1305045/1‑92 (MQ=255)

                          cAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatca                                                                >  2:6355418/1‑93 (MQ=255)

                           aGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatcat                                                               >  2:1648831/1‑93 (MQ=255)

                           aGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatcat                                                               <  2:5420159/93‑1 (MQ=255)

                            gATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatcatc                                                              >  1:6954190/1‑93 (MQ=255)

                            gATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGCGTTGGTTGTTGTTTGACGACATTAatcatc                                                              >  1:1374184/1‑93 (MQ=255)

                             aTGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATcc                                                             >  2:9566057/1‑93 (MQ=255)

                             aTGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATcc                                                             >  1:2391696/1‑93 (MQ=255)

                              tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCa                                                            <  1:515593/93‑1 (MQ=255)

                                aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACg                                                          <  1:1074041/93‑1 (MQ=255)

                                aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGAAATTAATCATCCACg                                                          >  2:2103667/1‑93 (MQ=255)

                                 agagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACgg                                                         <  1:9194304/93‑1 (MQ=255)

                                 agagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGGTTGACGACATTAATCATCCACgg                                                         >  2:765300/1‑93 (MQ=255)

                                  gagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGt                                                        >  1:8182972/1‑93 (MQ=255)

                                    gTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTTACGACATTAATCATCCACGGTAc                                                      >  1:3952219/1‑93 (MQ=255)

                                    gTTAACGGGAAGTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTTTTGGTTGTTGTTTGACGACATTAATCATCCACGGTAc                                                      <  1:8307403/93‑1 (MQ=255)

                                     ttAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACg                                                     >  1:1160841/1‑93 (MQ=255)

                                      tAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑c                                                   <  1:961303/93‑1 (MQ=255)

                                         cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCaa                                                >  1:7018103/1‑93 (MQ=255)

                                          gggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCaaa                                               >  2:6159605/1‑93 (MQ=255)

                                          gggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCaaa                                               <  2:3141930/93‑1 (MQ=255)

                                            gAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAAtt                                             <  1:5443115/93‑1 (MQ=255)

                                              aTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTa                                           >  2:9864898/1‑93 (MQ=255)

                                              aTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTa                                           <  1:8490754/93‑1 (MQ=255)

                                              aTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTTTGTTGGTTGTTGTTTGACGACATTAATTCTCCACGGTACGCCCAAAttt                                            >  1:1906431/1‑93 (MQ=255)

                                                tCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATc                                         <  2:6934029/93‑1 (MQ=255)

                                                    tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTg                                     >  1:5810234/1‑93 (MQ=255)

                                                     gACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGa                                    <  1:1002797/93‑1 (MQ=255)

                                                     gACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGa                                    <  2:5273339/93‑1 (MQ=255)

                                                     gACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGa                                    >  2:8447749/1‑93 (MQ=255)

                                                      aCAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGAc                                   <  2:2385467/93‑1 (MQ=255)

                                                       cAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACt                                  <  2:7815282/93‑1 (MQ=255)

                                                         gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACttt                                <  2:2726610/93‑1 (MQ=255)

                                                         gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACttt                                <  2:8876332/93‑1 (MQ=255)

                                                          tGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTg                               >  2:4575286/1‑93 (MQ=255)

                                                           gCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGc                              <  2:5886875/93‑1 (MQ=255)

                                                             gggAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGGACG‑CCAAATTTATCGGTGACTTTGCCa                            >  1:4140666/1‑93 (MQ=255)

                                                              ggAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCaa                           >  1:3136134/1‑93 (MQ=255)

                                                               gAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCaaa                          <  2:1219167/93‑1 (MQ=255)

                                                                aGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAg                         <  1:729980/93‑1 (MQ=255)

                                                                 gggcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGc                        >  2:1406476/1‑93 (MQ=255)

                                                                  ggcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGcc                       >  2:7973510/1‑93 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  2:4897690/93‑1 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  2:2938999/93‑1 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  1:1838711/93‑1 (MQ=255)

                                                                      ggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATgc                   >  2:6748577/1‑93 (MQ=255)

                                                                          tCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc               >  2:7528172/1‑93 (MQ=255)

                                                                           cccGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCa              >  1:9965016/1‑93 (MQ=255)

                                                                           cccGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCa              <  1:4587084/93‑1 (MQ=255)

                                                                             cggcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGa            >  2:3870917/1‑93 (MQ=255)

                                                                               gcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa          >  1:1771751/1‑93 (MQ=255)

                                                                                cggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAg         >  2:5739894/1‑93 (MQ=255)

                                                                                  ggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGtt       >  2:8141580/1‑93 (MQ=255)

                                                                                  ggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGtt       >  2:2595649/1‑93 (MQ=255)

                                                                                   gTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGttt      >  1:8787808/1‑93 (MQ=255)

                                                                                    ttgttGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTc     >  1:6499373/4‑93 (MQ=255)

                                                                                    ttATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTc     <  2:4147357/93‑1 (MQ=255)

                                                                                     tATTGTTTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  2:5381422/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    >  1:6453036/1‑93 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  1:2044970/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  1:1621636/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  1:3696486/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  2:9086259/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  2:4355617/93‑1 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGCCTTTGCCAAAGCCCTACGCCCAGAACGTTTTCt   >  2:7314990/1‑93 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt   <  1:6189676/93‑1 (MQ=255)

CGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACG‑CCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTG  >  NC_000913/4322792‑4322971

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