BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,498	T→C	V37A (GTG→GCG)	yahI →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,498	0	T	C	99.4%	538.6 / ‑6.2	175	V37A (GTG→GCG)	yahI	putative kinase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (81/92); total (83/92)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 188 total aligned reads displayed.

TTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGTGCTGGAAATGCTGGCTTCCGATTACGACATTGTGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGATTG  >  NC_000913/339408‑339587

ttGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGgcgc                                                                                         >  2:766352/1‑93 (MQ=255)

  gCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTg                                                                                       <  1:1793122/93‑1 (MQ=255)

   ccATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTgg                                                                                      >  2:2907793/1‑93 (MQ=255)

     aTTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaa                                                                                    >  1:4494984/1‑93 (MQ=255)

     aTTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaa                                                                                    >  2:2590340/1‑93 (MQ=255)

      ttGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGTGCTGGaaa                                                                                   >  1:2837173/1‑93 (MQ=255)

      ttGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaaa                                                                                   <  2:3814469/93‑1 (MQ=255)

      ttGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaaa                                                                                   <  2:4006709/93‑1 (MQ=255)

      tgggtgcCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGaaa                                                                                   <  1:530750/91‑1 (MQ=255)

       tggtggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAAt                                                                                  <  2:172129/93‑1 (MQ=255)

       tggtggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAAt                                                                                  <  1:194360/93‑1 (MQ=255)

         gtggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGc                                                                                <  1:4507097/93‑1 (MQ=255)

          tggCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCt                                                                               <  2:4635267/93‑1 (MQ=255)

              aaCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCt                                                                           >  2:1418811/1‑93 (MQ=255)

               aCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCtt                                                                          >  2:3510947/1‑93 (MQ=255)

               aCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCtt                                                                          >  1:2113764/1‑93 (MQ=255)

               aCAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCtt                                                                          <  2:2857688/93‑1 (MQ=255)

                cAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTc                                                                         <  1:1675571/93‑1 (MQ=255)

                 aGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTcc                                                                        <  1:2494751/93‑1 (MQ=255)

                  gCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAACTGCTGGCTTCCg                                                                       >  2:386731/1‑93 (MQ=255)

                  gCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCg                                                                       <  1:3019281/93‑1 (MQ=255)

                    aTTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAt                                                                     >  1:386076/1‑93 (MQ=255)

                    aTTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAt                                                                     >  2:4510243/1‑93 (MQ=255)

                     ttATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGAtt                                                                    <  2:363161/93‑1 (MQ=255)

                         cAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGa                                                                <  1:4582258/93‑1 (MQ=255)

                         cAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGa                                                                <  1:3468388/93‑1 (MQ=255)

                         cAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGa                                                                <  2:474458/93‑1 (MQ=255)

                          aaaGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGAc                                                               <  1:2744175/93‑1 (MQ=255)

                          aaaGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGAc                                                               <  1:2202373/93‑1 (MQ=255)

                           aaGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACa                                                              <  1:1246030/93‑1 (MQ=255)

                           aaGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACa                                                              <  1:739018/93‑1 (MQ=255)

                           aaGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACa                                                              >  2:2793152/1‑93 (MQ=255)

                            aGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACAt                                                             >  1:3553655/1‑93 (MQ=255)

                            aGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACAt                                                             >  2:2635948/1‑93 (MQ=255)

                            aGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACAt                                                             <  1:176436/93‑1 (MQ=255)

                            aGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACAt                                                             >  2:4404546/1‑93 (MQ=255)

                              aTAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtg                                                           >  1:822893/1‑93 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGGGCTGGAAATGATGGCTTCCGATTACGACATtgg                                                          >  1:496364/1‑92 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                          <  2:4579852/93‑2 (MQ=255)

                               tAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATtgc                                                          >  2:330532/1‑92 (MQ=255)

                                 aCGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGGAATGCTGGCTTCCGGTTACGACATTgcgc                                                        >  1:1271296/1‑93 (MQ=255)

                                 aCGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTgcgc                                                        <  2:2241411/93‑1 (MQ=255)

                                     cagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAc                                                    <  2:1761431/93‑1 (MQ=255)

                                     cagccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAc                                                    >  1:467042/1‑93 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                   >  2:3139082/1‑93 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                   >  2:2413228/1‑93 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                   >  1:3486166/1‑93 (MQ=255)

                                      agccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAcc                                                   <  2:4460238/93‑1 (MQ=255)

                                       gccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAccc                                                  >  2:1863806/1‑93 (MQ=255)

                                       gccagTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGAccc                                                  >  1:1464431/1‑93 (MQ=255)

                                          agTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACg                                               <  1:3824365/93‑1 (MQ=255)

                                          agTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACg                                               >  2:1624449/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGAGATTGCGCTGACCCACgg                                              >  1:3904786/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACgg                                              >  1:1513854/1‑93 (MQ=255)

                                           gTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATNACGACATTGCGCTGACCCANgg                                              >  1:3258876/1‑93 (MQ=255)

                                              gATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCaa                                           <  1:3313661/93‑1 (MQ=255)

                                               aTTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAAc                                          <  2:57361/93‑1 (MQ=255)

                                               aTTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAAc                                          >  1:1348506/1‑93 (MQ=255)

                                                 tGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACgg                                        <  1:3238796/93‑1 (MQ=255)

                                                   aGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGc                                      <  2:694457/93‑1 (MQ=255)

                                                    gCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGcc                                     <  1:4210136/93‑1 (MQ=255)

                                                       tCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCa                                  >  1:3632575/1‑93 (MQ=255)

                                                       tCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCCGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCa                                  <  1:1028169/93‑1 (MQ=255)

                                                        cAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAg                                 <  1:1318782/93‑1 (MQ=255)

                                                         aggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                                >  2:1338336/1‑93 (MQ=255)

                                                         aggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                                <  1:4102663/93‑1 (MQ=255)

                                                         aggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                                <  2:4416256/93‑1 (MQ=255)

                                                         aggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                                <  1:2878179/93‑1 (MQ=255)

                                                         aggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAgg                                >  1:1073076/1‑93 (MQ=255)

                                                          ggcggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGt                               <  2:3850490/93‑1 (MQ=255)

                                                            cggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCg                             <  2:4176133/93‑1 (MQ=255)

                                                            cggaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCg                             <  1:2340258/93‑1 (MQ=255)

                                                              gaggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCggg                           <  1:1133263/93‑1 (MQ=255)

                                                               aggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGc                          >  1:1068156/1‑93 (MQ=255)

                                                                ggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCt                         <  1:160980/93‑1 (MQ=255)

                                                                ggcggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCt                         <  2:4000552/93‑1 (MQ=255)

                                                                  cggTGAAAGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTgg                       <  1:2414042/93‑1 (MQ=255)

                                                                       aaaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTa                  >  2:2943213/1‑93 (MQ=255)

                                                                       aaaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTa                  <  2:2472940/93‑1 (MQ=255)

                                                                       aaaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTa                  <  2:1087362/93‑1 (MQ=255)

                                                                        aaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGTCTGGATTTAc                 <  2:3048783/93‑1 (MQ=255)

                                                                        aaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTAc                 <  2:3582932/93‑1 (MQ=255)

                                                                        aaGCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTAc                 >  2:3742176/1‑93 (MQ=255)

                                                                          gCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGc               >  1:40174/1‑93 (MQ=255)

                                                                          gCCGTCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGc               >  1:3624628/1‑93 (MQ=255)

                                                                          gCCCTCGCCGATACGGCGCTGGAAATGCTGTCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCTCAGGTCGGGCTGGATTTACGc               <  2:2078101/93‑1 (MQ=255)

                                                                              tCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTg           <  1:3956670/93‑1 (MQ=255)

                                                                              tCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTg           >  2:92960/1‑93 (MQ=255)

                                                                              tCGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTg           >  2:3943788/1‑93 (MQ=255)

                                                                               cGCCGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGc          <  1:3426388/93‑1 (MQ=255)

                                                                                 ccGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCgg        >  1:2596534/1‑93 (MQ=255)

                                                                                 ccGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCgg        >  1:1506861/1‑93 (MQ=255)

                                                                                  cGATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGga       <  1:427832/93‑1 (MQ=255)

                                                                                   gATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgag      >  1:2062001/1‑93 (MQ=255)

                                                                                   gATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgag      >  2:1901049/1‑93 (MQ=255)

                                                                                   gATACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgag      >  2:970887/1‑93 (MQ=255)

                                                                                    aTACGGCGCTGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgaga     <  1:1621066/93‑1 (MQ=255)

                                                                                    aTACGGCGCGGGAAATGCTGGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGgaga     <  1:3217128/93‑1 (MQ=255)

                                                                                      aCGGCGCTGGAAATGCTGGCGTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGAtt   <  2:4253814/93‑1 (MQ=255)

                                                                                       cGGCGCTGGAAATGCTCGCTTCCGATTACGACATTGCGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGCTTTACGCCGTGCGGAGAGTg  >  2:614765/1‑93 (MQ=38)

TTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCCGTCGCCGATACGGTGCTGGAAATGCTGGCTTCCGATTACGACATTGTGCTGACCCACGGCAACGGGCCGCAGGTCGGGCTGGATTTACGCCGTGCGGAGATTG  >  NC_000913/339408‑339587

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