BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,580	T→A	I389N (ATC→AAC)	katG →	catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,580	0	T	A	98.8%	927.7 / ‑6.2	262	I389N (ATC→AAC)	katG	catalase
Reads supporting (aligned to +/- strand): ref base T (0/1); major base A (156/100); minor base G (0/3); total (156/106)
Fisher's exact test for biased strand distribution p-value = 6.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

TCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGA  >  NC_000913/4132490‑4132669

tCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAact                                                                                         >  1:4154136/1‑93 (MQ=255)

 cccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactc                                                                                        >  1:5480024/1‑93 (MQ=255)

  ccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGTAGAactct                                                                                       <  1:6363256/93‑1 (MQ=255)

  ccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactct                                                                                       >  2:319818/1‑93 (MQ=255)

   cGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactctc                                                                                      >  2:1263639/1‑93 (MQ=255)

    ggATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcg                                                                                     >  2:1866420/1‑93 (MQ=255)

     gATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgt                                                                                    <  1:2924801/93‑1 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGCCCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   <  1:4623405/93‑1 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   >  1:1595209/1‑93 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   <  2:4954072/93‑1 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   >  2:22095/1‑93 (MQ=255)

        ccGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTcgtcgt                                                                                 >  2:324093/1‑93 (MQ=255)

           tttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTc                                                                              >  1:3820542/1‑93 (MQ=255)

            ttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTcc                                                                             <  1:4440315/93‑1 (MQ=255)

             tGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCt                                                                            <  2:62811/93‑1 (MQ=255)

                tCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCaa                                                                         >  1:6288874/1‑93 (MQ=255)

                 ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAAc                                                                        <  2:3415448/93‑1 (MQ=255)

                  cGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACg                                                                       >  2:5436025/1‑93 (MQ=255)

                  cGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACg                                                                       <  1:5212479/93‑1 (MQ=255)

                   gTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGa                                                                      <  2:3512254/93‑1 (MQ=255)

                    tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCCACGAt                                                                     >  1:4586175/1‑93 (MQ=255)

                    tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGAt                                                                     <  2:4824297/93‑1 (MQ=255)

                      gaagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                   >  1:6158748/1‑93 (MQ=255)

                      gaagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                   <  1:6014378/93‑1 (MQ=255)

                       aagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCg                                                                  <  1:2939326/93‑1 (MQ=255)

                       aagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCg                                                                  >  2:5770925/1‑93 (MQ=255)

                        agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGc                                                                 <  1:5094891/93‑1 (MQ=255)

                        agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGc                                                                 <  2:5521657/93‑1 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                >  2:1411919/1‑93 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                <  1:2066135/93‑1 (MQ=255)

                          aaaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAg                                                               >  2:6383759/1‑93 (MQ=255)

                          aaaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAg                                                               <  1:1172345/93‑1 (MQ=255)

                           aaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAgg                                                              >  1:3575103/1‑93 (MQ=255)

                            aCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGc                                                             <  2:2752064/93‑1 (MQ=255)

                              gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGt                                                           <  1:1120674/93‑1 (MQ=255)

                              gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGt                                                           <  1:1372079/93‑1 (MQ=255)

                              gTAAACCGACAATGCTGGTGACCAACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGt                                                           >  2:1051439/1‑93 (MQ=255)

                                aaaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTc                                                         >  1:2236666/1‑93 (MQ=255)

                                aaaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTc                                                         >  1:273067/1‑93 (MQ=255)

                                aaaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTc                                                         >  2:4186729/1‑93 (MQ=255)

                                 aaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCa                                                        >  1:299165/1‑93 (MQ=255)

                                 aaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCa                                                        <  2:5619565/93‑1 (MQ=255)

                                  aCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCaa                                                       >  2:2135667/1‑93 (MQ=255)

                                    cGACAATGCTGGTGACCGACCTGACGCTGCTTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACg                                                     <  2:5087433/93‑1 (MQ=255)

                                    cGACAATGCTGGTGACCGACCTGACGCTGCG‑TTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGa                                                    <  1:1537529/93‑1 (MQ=255)

                                     gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGa                                                    >  2:386142/1‑93 (MQ=255)

                                     gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGa                                                    >  2:1810977/1‑93 (MQ=255)

                                      aCAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGaa                                                   <  1:3452559/93‑1 (MQ=255)

                                      aCAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGaa                                                   >  2:824029/1‑93 (MQ=255)

                                      aCAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGaa                                                   >  1:4909112/1‑93 (MQ=255)

                                       cAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAg                                                  <  2:1162263/93‑1 (MQ=255)

                                         aTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCCACGATCCGCAGGCGTTCAACGaacct                                                >  2:2343268/1‑90 (MQ=255)

                                         aTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGcc                                                >  1:4739854/1‑93 (MQ=255)

                                          tGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCt                                               >  2:4897154/1‑93 (MQ=255)

                                            cTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCttt                                             >  2:13384/1‑93 (MQ=255)

                                            cTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCttt                                             <  1:3394929/93‑1 (MQ=255)

                                            cTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCttt                                             <  1:3999034/93‑1 (MQ=255)

                                               gTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGcc                                          <  1:3871043/93‑1 (MQ=255)

                                                tcaccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGccc                                         <  2:1556807/91‑1 (MQ=255)

                                                tGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGccc                                         <  2:685134/93‑1 (MQ=255)

                                                  accgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGt                                       >  2:5732164/1‑93 (MQ=255)

                                                  accgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGt                                       >  1:2696040/1‑93 (MQ=255)

                                                   ccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTg                                      >  1:1798274/1‑93 (MQ=255)

                                                     gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    >  1:2549963/1‑93 (MQ=255)

                                                     gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    <  1:4987641/93‑1 (MQ=255)

                                                      accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                   <  2:3069180/93‑1 (MQ=255)

                                                      accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                   >  2:2401731/1‑93 (MQ=255)

                                                      accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                   >  1:720245/1‑93 (MQ=255)

                                                      accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                   >  1:5317851/1‑93 (MQ=255)

                                                       ccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTg                                  <  1:5045962/93‑1 (MQ=255)

                                                         tGACGCTGCGTTTTGAGCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                <  1:1700321/93‑1 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               >  1:2406412/1‑93 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               <  2:3459128/93‑1 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               <  2:2786995/93‑1 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               >  2:4731106/1‑93 (MQ=255)

                                                            cGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCa                             >  2:5942376/1‑93 (MQ=255)

                                                             gCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCaa                            >  2:6440055/1‑93 (MQ=255)

                                                                  gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                       >  1:377370/1‑93 (MQ=255)

                                                                  gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                       >  1:1688272/1‑93 (MQ=255)

                                                                   ttttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAc                      >  1:5634055/1‑93 (MQ=255)

                                                                       gATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcac                  <  2:920936/93‑1 (MQ=255)

                                                                       gATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcac                  >  2:2592218/1‑93 (MQ=255)

                                                                        aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                 >  2:4566911/1‑93 (MQ=255)

                                                                        aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                 >  2:1932216/1‑93 (MQ=255)

                                                                        aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                 >  2:6109382/1‑93 (MQ=255)

                                                                          ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg               >  2:3943051/1‑93 (MQ=255)

                                                                          ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg               >  2:3627395/1‑93 (MQ=255)

                                                                           cTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAggg              <  1:1964033/93‑1 (MQ=255)

                                                                            tGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGa             >  2:5276010/1‑93 (MQ=255)

                                                                            tGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGa             <  2:3359143/93‑1 (MQ=255)

                                                                              aGTTCGAGAAGAACGCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGata           <  2:5218229/93‑1 (MQ=255)

                                                                               gTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGatat          >  1:972600/1‑93 (MQ=255)

                                                                               gTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGatat          >  2:2616396/1‑93 (MQ=255)

                                                                                ttCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATg         <  2:4348340/93‑1 (MQ=255)

                                                                                   gagaAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgggg      <  2:1766316/93‑1 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCg   >  2:4062587/1‑93 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCg   >  1:4338402/1‑93 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCg   <  1:2473212/93‑1 (MQ=255)

                                                                                       agaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGa  <  1:5878934/93‑1 (MQ=255)

                                                                                       agaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGCATATGGGGCCGa  <  1:6184644/93‑1 (MQ=255)

TCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGA  >  NC_000913/4132490‑4132669

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