BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,061,601	T→C	C194R (TGC→CGC)	torD →	TorD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,061,601	0	T	C	99.5%	685.8 / ‑4.8	214	C194R (TGC→CGC)	torD	TorD
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (97/116); total (97/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 227 total aligned reads displayed.

TGGTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCG  >  NC_000913/1061511‑1061690

   ttACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGac                                                                                      <  2:5245419/93‑1 (MQ=255)

   ttACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGac                                                                                      >  1:2151056/1‑93 (MQ=255)

    tACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGacc                                                                                     <  2:6076740/93‑1 (MQ=255)

    tACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGacc                                                                                     <  2:576617/93‑1 (MQ=255)

      ccAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGaccac                                                                                   >  2:3872892/1‑93 (MQ=255)

       cAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGaccacc                                                                                  >  1:3938319/1‑93 (MQ=255)

       cAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGaccacc                                                                                  >  1:4025541/1‑93 (MQ=255)

         gagTTTGTTGCGCGTTGTCTTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCaa                                                                                <  1:2698279/93‑1 (MQ=255)

           gTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCaaaa                                                                              >  1:2481516/1‑93 (MQ=255)

            tttGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAAc                                                                             >  2:1667607/1‑93 (MQ=255)

            tttGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAAc                                                                             >  2:6670126/1‑93 (MQ=255)

              tgttgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAg                                                                           >  1:2298639/1‑93 (MQ=255)

               gttgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGa                                                                          <  1:1629189/93‑1 (MQ=255)

                ttgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGAt                                                                         <  2:235119/93‑1 (MQ=255)

                 tgCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATa                                                                        <  1:2238593/93‑1 (MQ=255)

                  gcgcGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATaa                                                                       <  2:4467079/93‑1 (MQ=255)

                    gcgTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcg                                                                     >  2:5095230/1‑93 (MQ=255)

                     cgTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgt                                                                    <  2:2485691/93‑1 (MQ=255)

                     cgTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgt                                                                    <  1:3070907/93‑1 (MQ=255)

                      gTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgtc                                                                   >  1:1978167/1‑93 (MQ=255)

                      gTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgtc                                                                   <  1:279547/93‑1 (MQ=255)

                      gTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGCGGAGCGCGACCACCAAAACAGATAAcgtc                                                                   >  1:6192136/1‑93 (MQ=255)

                        tGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAAcgtcgt                                                                 >  2:2586906/1‑93 (MQ=255)

                        tGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACAACCAAAACAGATAAcgtcgt                                                                 >  2:2896101/1‑93 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  2:4134821/93‑1 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  1:4953639/93‑1 (MQ=255)

                         gtcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGtt                                                                <  2:5819596/93‑1 (MQ=255)

                          tcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGttt                                                               >  1:5536877/1‑93 (MQ=255)

                          tcgtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGttt                                                               >  2:5570641/1‑93 (MQ=255)

                            gtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTtgt                                                             >  1:2806565/1‑93 (MQ=255)

                            gtcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTtgt                                                             <  2:2201811/93‑1 (MQ=255)

                             tcAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTtgtg                                                            >  1:1929535/1‑93 (MQ=255)

                              cAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgc                                                           <  2:819850/93‑1 (MQ=255)

                               aGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcg                                                          >  1:1052053/1‑93 (MQ=255)

                               aGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcg                                                          <  2:4315250/93‑1 (MQ=255)

                                gTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTgcgc                                                         >  2:2391628/1‑93 (MQ=255)

                                 tATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGcc                                                        <  2:2682151/93‑1 (MQ=255)

                                 tATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGcc                                                        <  2:3968052/93‑1 (MQ=255)

                                 tATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGCTAGTGGAGCGCGACCACCAAAACAGATAACGGCGTTTGTGCGcc                                                        >  2:5120852/1‑93 (MQ=255)

                                  aTGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCt                                                       <  2:1584897/93‑1 (MQ=255)

                                   tGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCCAAACAGATAACGTCGTTTGTGCGCCTg                                                      >  1:5484942/1‑93 (MQ=255)

                                   tGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTg                                                      <  1:6697224/93‑1 (MQ=255)

                                   tGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTg                                                      >  2:1341376/1‑93 (MQ=255)

                                    gACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGa                                                     <  1:3244189/93‑1 (MQ=255)

                                     aCAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaa                                                    >  1:1377707/1‑93 (MQ=255)

                                       aGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaaaa                                                  >  2:3435740/1‑93 (MQ=255)

                                       aGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaaaa                                                  <  2:4164935/93‑1 (MQ=255)

                                       aGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGaaaa                                                  <  2:1755357/93‑1 (MQ=255)

                                          tttGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAc                                               >  1:6326867/1‑93 (MQ=255)

                                          tttGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAc                                               <  1:6257271/93‑1 (MQ=255)

                                           ttGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcg                                              <  1:127661/93‑1 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgcg                                            >  1:255577/1‑93 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgcg                                            <  1:5749020/93‑1 (MQ=255)

                                             ggTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGAcgcg                                            >  2:4036851/1‑93 (MQ=255)

                                               tttttACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGtt                                          <  1:5242449/93‑1 (MQ=255)

                                                 tttACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTa                                        <  1:5616906/93‑1 (MQ=255)

                                                   tACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGc                                      <  2:1007645/93‑1 (MQ=255)

                                                   tACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGc                                      >  1:827718/1‑93 (MQ=255)

                                                    aCGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCg                                     >  2:3528601/1‑93 (MQ=255)

                                                     cgaggCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGt                                    <  2:2165913/90‑1 (MQ=255)

                                                      gcgGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTc                                   >  2:605789/1‑93 (MQ=255)

                                                      gcgGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTc                                   <  1:5153945/93‑1 (MQ=255)

                                                        ggCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGc                                 >  2:192962/1‑93 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTGGCATc                              >  1:1847300/1‑93 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATc                              <  2:3050409/93‑1 (MQ=255)

                                                           aCTAAGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATc                              <  2:5161233/93‑1 (MQ=255)

                                                              aaGCCAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAgg                           >  2:3473491/1‑93 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  1:5093975/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  2:4958000/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  1:3302933/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        <  2:1107800/93‑1 (MQ=255)

                                                                 ccAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAt                        >  2:5991886/1‑93 (MQ=255)

                                                                  cAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAtt                       <  1:5899579/93‑1 (MQ=255)

                                                                  cAGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCAtt                       <  2:2055068/93‑1 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      <  2:5687777/93‑1 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      <  1:1621008/93‑1 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      >  1:4882328/1‑93 (MQ=255)

                                                                   aGTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTa                      <  1:6010472/93‑1 (MQ=255)

                                                                    gTTATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTAt                     >  2:3248786/1‑93 (MQ=255)

                                                                     ttATTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATg                    >  2:6021751/1‑93 (MQ=255)

                                                                       aTTGCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCCGGCATTATgcc                  >  1:1476028/1‑91 (MQ=255)

                                                                          gCTGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCa               <  2:4521623/93‑1 (MQ=255)

                                                                            tGGTGTTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGt             <  1:353387/93‑1 (MQ=255)

                                                                               tgtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGc          <  2:4800927/93‑1 (MQ=255)

                                                                               tgtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGc          <  1:451315/93‑1 (MQ=255)

                                                                               tgtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGc          >  2:5602878/1‑93 (MQ=255)

                                                                               tgtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTACGGCGCAGTTGc          <  2:3751112/93‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  1:867144/93‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         >  2:6357496/1‑93 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:2782594/93‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:3362682/93‑1 (MQ=255)

                                                                                gtTAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGcc         <  2:40879/93‑1 (MQ=255)

                                                                                 ttAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCg        >  2:5912876/1‑93 (MQ=255)

                                                                                 ttAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCg        <  1:4478913/93‑1 (MQ=255)

                                                                                 ttAGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCg        <  1:4173405/93‑1 (MQ=255)

                                                                                   aGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGa      <  2:3005403/93‑1 (MQ=255)

                                                                                   aGTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGa      <  1:2416967/93‑1 (MQ=255)

                                                                                    gTGGAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGAt     <  1:3611380/93‑1 (MQ=255)

                                                                                     tGGAGCGCCACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGCCATTATGGGGCAGTTGCCGGATg    <  1:4844174/93‑1 (MQ=39)

                                                                                      ggAGCGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATgc   <  1:3996521/93‑1 (MQ=255)

TGGTTACCAGAGTTTGTTGCGCGTTGTCGTCAGTATGACAGCTTTGGTTTTTACGCGGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCG  >  NC_000913/1061511‑1061690

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: `‑`