BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	592,232	C→T	R57Q (CGA→CAA)	nfrB ←	NfrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	592,232	0	C	T	99.4%	606.0 / NA	175	R57Q (CGA→CAA)	nfrB	NfrB
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (114/60); minor base G (0/1); total (114/61)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

ATATTGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTCGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCTGATG  >  NC_000913/592142‑592321

 tatTGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTtgcg                                                                                        >  2:2123465/1‑93 (MQ=255)

 tatTGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTtgcg                                                                                        <  1:3918378/93‑1 (MQ=255)

 tatTGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTtgcg                                                                                        >  1:4179066/1‑93 (MQ=255)

    tGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGt                                                                                     <  2:1823351/93‑1 (MQ=255)

     gCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTa                                                                                    <  1:2824981/93‑1 (MQ=255)

       cGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAgc                                                                                  >  2:1500597/1‑93 (MQ=255)

        gATGACGCCCGTTTCATTCCACTCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAgcg                                                                                 >  1:844625/1‑93 (MQ=255)

         aTGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAgcgg                                                                                >  2:3729848/1‑93 (MQ=255)

           gACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAgcggcg                                                                              >  1:3746930/1‑93 (MQ=255)

            aCGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGa                                                                             <  2:4321722/93‑1 (MQ=255)

             cGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGAt                                                                            >  2:3173910/1‑93 (MQ=255)

               cccGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATaa                                                                          >  1:3788916/1‑93 (MQ=255)

               cccGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATaa                                                                          <  1:4392057/93‑1 (MQ=255)

               cccGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATaa                                                                          >  2:4536529/1‑93 (MQ=255)

                ccGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATaaa                                                                         >  2:1529122/1‑93 (MQ=255)

                 cGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAac                                                                        >  2:1747840/1‑93 (MQ=255)

                 cGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAac                                                                        >  1:2367363/1‑93 (MQ=255)

                  ggttCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGATAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       >  2:2253984/3‑93 (MQ=255)

                  gTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       >  1:366632/1‑93 (MQ=255)

                  gTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       >  1:3401818/1‑93 (MQ=255)

                  gTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       <  1:2053150/93‑1 (MQ=255)

                  gTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       <  1:4208898/93‑1 (MQ=255)

                  gTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAaca                                                                       >  2:501811/1‑93 (MQ=255)

                   tttCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAacac                                                                      >  2:2341650/1‑93 (MQ=255)

                       aTTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCaa                                                                  >  1:642526/1‑93 (MQ=255)

                         tCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACt                                                                >  2:4406524/1‑93 (MQ=255)

                          ccACGCCGGAACAATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACCCTCAACtt                                                               <  1:3185904/93‑1 (MQ=255)

                             cGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCg                                                            <  1:11499/93‑1 (MQ=255)

                             cGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCg                                                            >  2:2733298/1‑93 (MQ=255)

                               ccGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCgct                                                          >  1:2196711/1‑91 (MQ=255)

                                ccgAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGttt                                                         <  1:1564892/91‑1 (MQ=255)

                                cGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGttt                                                         >  2:397354/1‑93 (MQ=255)

                                cGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGttt                                                         >  2:1401322/1‑93 (MQ=255)

                                cGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGttt                                                         >  2:1917801/1‑93 (MQ=255)

                                 ggAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTa                                                        <  2:3212936/93‑1 (MQ=255)

                                 ggAACCATAATCGCTAAAGATTTTTCATCTGGTTTATACAGTTCGCGATAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTa                                                        <  2:4214682/93‑1 (MQ=255)

                                  gAACTATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTaa                                                       >  2:1930625/1‑93 (MQ=255)

                                  gAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTaa                                                       >  2:4157569/1‑93 (MQ=255)

                                  gAACCATAATCGCGACTGGTCTTACATCTGATTTATACAGTTCGCGATATCGCATTCGCGGGTAGCGGCGATAACCACTCAACTTGCGTTTaa                                                       <  2:91720/93‑1 (MQ=255)

                                    aCCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATg                                                     >  2:1870114/1‑93 (MQ=255)

                                    aCCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATg                                                     >  1:2873564/1‑93 (MQ=255)

                                     ccATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgc                                                    >  1:4640029/1‑93 (MQ=255)

                                     ccATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgc                                                    >  1:956246/1‑93 (MQ=255)

                                      cATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgcg                                                   <  1:448641/93‑1 (MQ=255)

                                       aTAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgcgg                                                  >  2:2825064/1‑93 (MQ=255)

                                         aaTCGCTAACTGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgcggcg                                                >  2:3023018/1‑93 (MQ=255)

                                         aaTCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATgcggcg                                                <  1:185173/93‑1 (MQ=255)

                                          aTCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGt                                               >  1:2052259/1‑93 (MQ=255)

                                              cTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTAccc                                           >  2:2682476/1‑93 (MQ=255)

                                              cTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTAccc                                           >  2:2922878/1‑93 (MQ=255)

                                               tAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCa                                          >  1:551497/1‑93 (MQ=255)

                                               tAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCa                                          >  1:488694/1‑93 (MQ=255)

                                               tAACGGTTTTTCATATGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCa                                          >  1:4388267/1‑93 (MQ=255)

                                                aaCGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAg                                         >  2:2465519/1‑93 (MQ=255)

                                                 aCGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGt                                        <  2:4392416/93‑1 (MQ=255)

                                                  cGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTa                                       <  2:569484/93‑1 (MQ=255)

                                                  cGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTa                                       <  1:2726389/93‑1 (MQ=255)

                                                     tttttCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAgac                                    >  2:4116141/1‑93 (MQ=255)

                                                     tttttCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAgac                                    <  1:3030447/93‑1 (MQ=255)

                                                       tttCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAgacga                                  <  1:579875/93‑1 (MQ=255)

                                                       tttCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAgacga                                  >  2:2598821/1‑93 (MQ=255)

                                                       tttCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAgacga                                  >  1:3559046/1‑93 (MQ=255)

                                                         tCATCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACa                                >  1:2066125/1‑93 (MQ=255)

                                                            tCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCa                             >  1:519279/1‑93 (MQ=255)

                                                            tCTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCa                             >  1:1768591/1‑93 (MQ=255)

                                                             cTGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCaa                            <  1:1358740/93‑1 (MQ=255)

                                                              tGGTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAAt                           >  1:2019463/1‑93 (MQ=255)

                                                               ggTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATa                          <  1:4611343/93‑1 (MQ=255)

                                                               ggTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATa                          >  2:2888886/1‑93 (MQ=255)

                                                               ggTGTATACAGTGCGCGGTACCTCATTCGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTCCCCAGTAGACGACATCAATa                          <  1:589300/93‑1 (MQ=255)

                                                                gTTTATACAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGACAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATaa                         >  1:26824/1‑93 (MQ=255)

                                                                   tataCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATaaaaa                      >  1:1590220/1‑93 (MQ=255)

                                                                    ataCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATaaaaaa                     <  1:2429932/93‑1 (MQ=255)

                                                                     taCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATaaaaaaa                    >  2:976712/1‑93 (MQ=255)

                                                                     taCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATaaaaaaa                    >  2:3561492/1‑93 (MQ=255)

                                                                      aCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAAt                   >  1:3535038/1‑93 (MQ=255)

                                                                      aCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAAt                   >  2:4546339/1‑93 (MQ=255)

                                                                      aCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAAt                   <  2:3259924/93‑1 (MQ=255)

                                                                      aCAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAAt                   >  1:2241310/1‑93 (MQ=255)

                                                                       cAGTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATc                  <  1:3422480/93‑1 (MQ=255)

                                                                       cAGTTCGCGGTAACGCATTGGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATc                  <  1:634811/93‑1 (MQ=255)

                                                                         gTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGt                <  2:868201/93‑1 (MQ=255)

                                                                         gTTCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGt                >  1:4070545/1‑93 (MQ=255)

                                                                         gTTCGCGGTAACTCATTTGCGGGGAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGt                >  1:3695506/1‑93 (MQ=255)

                                                                          ttCGCGGTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTc               <  1:103296/93‑1 (MQ=255)

                                                                               ggTAACTCATTTGCGGGTAGCGGCGATAAACCCTCAACTTGCGTTTAATGCGGCCTTCCCAGTAGACGACATCAATAAAAAAATCGTCCAGcc          <  1:1453603/93‑1 (MQ=255)

                                                                               ggTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGcc          <  2:1982576/93‑1 (MQ=255)

                                                                                gTAACTCATTTGCGGGTAGCGGCGATAAACACTCAGCTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGccc         <  1:2772303/93‑1 (MQ=255)

                                                                                gTAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGccc         <  1:3559782/93‑1 (MQ=255)

                                                                                 tAACTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCg        >  1:2096354/1‑93 (MQ=255)

                                                                                  aaCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGc       >  1:907703/1‑93 (MQ=255)

                                                                                  aaCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGc       >  1:3591297/1‑93 (MQ=255)

                                                                                  aaCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGc       >  1:3105678/1‑93 (MQ=255)

                                                                                   aCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCt      <  2:1176492/93‑1 (MQ=255)

                                                                                   aCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCt      >  1:3897536/1‑93 (MQ=255)

                                                                                   aCTCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCt      <  2:2425683/93‑1 (MQ=255)

                                                                                     tCATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCtga    >  1:3001115/1‑93 (MQ=255)

                                                                                      cATTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCtgat   <  1:1563171/93‑1 (MQ=255)

                                                                                       aTTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTCGACGACATCAATAAAAAAATCGTCCAGCCCGCtgatg  >  1:4149847/1‑93 (MQ=255)

                                                                                       aTTTGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCtgatg  >  2:477566/1‑93 (MQ=255)

ATATTGCCGATGACGCCCGTTTCATTCCACGCCGGAACCATAATCGCTAACGGTTTTTCATCTGGTTTATACAGTTCGCGGTAACTCATTCGCGGGTAGCGGCGATAAACACTCAACTTGCGTTTAATGCGGCGTACCCAGTAGACGACATCAATAAAAAAATCGTCCAGCCCGCTGATG  >  NC_000913/592142‑592321

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