BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	620,296	A→G	V39A (GTA→GCA)	fepG ←	FepG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	620,296	0	A	G	99.0%	659.6 / NA	200	V39A (GTA→GCA)	fepG	FepG
Reads supporting (aligned to +/- strand): ref base A (1/0); major base G (110/87); minor base T (0/2); total (111/89)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 205 total aligned reads displayed.

ATCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGTGATG  >  NC_000913/620206‑620385

 tCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCt                                                                                        <  1:2992453/93‑1 (MQ=255)

  cAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTg                                                                                       >  2:1553453/1‑93 (MQ=255)

   aaCAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGc                                                                                      >  1:4455261/1‑93 (MQ=255)

      aGCGCCATCAGCACGCGTGGGAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAg                                                                                   >  1:3364388/1‑93 (MQ=255)

         gcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGtt                                                                                >  1:3733373/1‑93 (MQ=255)

          ccATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGttt                                                                               >  2:3536631/1‑93 (MQ=255)

           cATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTc                                                                              <  2:4987105/93‑1 (MQ=255)

           cATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTc                                                                              >  1:2917291/1‑93 (MQ=255)

            aTCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTcc                                                                             >  1:4631167/1‑93 (MQ=255)

            aTCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTcc                                                                             >  2:5711850/1‑93 (MQ=255)

             tCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTccg                                                                            >  2:5308686/1‑92 (MQ=255)

             tCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCa                                                                            <  1:1167296/93‑1 (MQ=255)

             tCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCa                                                                            >  1:4381645/1‑93 (MQ=255)

             tCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCa                                                                            <  1:3745730/93‑1 (MQ=255)

               aGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAgg                                                                          >  1:6086628/1‑92 (MQ=255)

               aGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGc                                                                          >  2:520366/1‑93 (MQ=255)

               aGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGc                                                                          <  1:422420/93‑1 (MQ=255)

               aGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGc                                                                          <  1:2078042/93‑1 (MQ=255)

                    gcgTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAc                                                                     <  1:4414577/93‑1 (MQ=255)

                    gcgTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAc                                                                     <  1:435357/93‑1 (MQ=255)

                     cgTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACg                                                                    <  1:2872149/93‑1 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   <  2:2307872/93‑1 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   <  1:3614216/93‑1 (MQ=255)

                      gTGGTAAACGCCATTAGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGACTGCGAGGCTTCCAGCGAGACgg                                                                   >  2:34391/1‑93 (MQ=255)

                          tAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCAcc                                                               <  2:5621046/93‑1 (MQ=255)

                            aaCGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGc                                                             >  2:2679183/1‑93 (MQ=255)

                              cGCCATTCGGTGCCCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTg                                                           <  2:4732524/93‑1 (MQ=255)

                              cGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATTCCTGCGAGGTTTCCAGCGTGACGGCACCGCTg                                                           <  2:1618482/93‑1 (MQ=255)

                              cGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTg                                                           >  1:261760/1‑93 (MQ=255)

                               gCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGc                                                          <  1:2152462/93‑1 (MQ=255)

                               gCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGc                                                          >  2:3219184/1‑93 (MQ=255)

                                ccATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCg                                                         >  1:4700209/1‑93 (MQ=255)

                                 cATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGt                                                        >  1:1905429/1‑93 (MQ=255)

                                 cATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGt                                                        >  1:2767799/1‑93 (MQ=255)

                                   ttCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTaa                                                      <  2:1113597/93‑1 (MQ=255)

                                    tCGGGGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAAt                                                     <  1:1509150/93‑1 (MQ=255)

                                     cGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATc                                                    <  1:687498/93‑1 (MQ=255)

                                         gACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCa                                                <  1:2369903/93‑1 (MQ=255)

                                         gACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCa                                                >  1:2185558/1‑93 (MQ=255)

                                          accaccATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAg                                               >  1:5097596/1‑93 (MQ=255)

                                            caccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGAt                                             >  1:5918313/1‑93 (MQ=255)

                                             accaTCGTCATACTGCGCGGCGCACCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGAAACGCTGCGTAATCCCCAGATa                                            >  2:2632512/1‑93 (MQ=255)

                                              ccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAc                                           >  2:1355339/1‑93 (MQ=255)

                                              ccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAc                                           >  1:1777277/1‑93 (MQ=255)

                                               caTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAcc                                          >  1:1306422/1‑93 (MQ=255)

                                                 tCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATACCCAGATACCTg                                        >  1:5984296/1‑93 (MQ=255)

                                                 tCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTg                                        >  1:5343112/1‑93 (MQ=255)

                                                  cGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGc                                       <  2:5612506/93‑1 (MQ=255)

                                                   gTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCa                                      >  2:4837052/1‑93 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  2:2571706/1‑93 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  1:103033/1‑93 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  1:5826607/1‑93 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCAGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  1:3202954/1‑93 (MQ=255)

                                                     cATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAc                                    <  1:165572/93‑1 (MQ=255)

                                                      aTACTGCGCGGCGCATCGCCCATCAGCGCGGCGATTGCATGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAcc                                   >  2:3919261/1‑93 (MQ=255)

                                                       tACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACca                                  >  2:1913875/1‑93 (MQ=255)

                                                       tACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGGTTCCAGCGTGACGGCACCGCTGCGTAGTCCCCAGATACCTGCAACca                                  >  2:6075052/1‑93 (MQ=255)

                                                       tACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGGTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACca                                  >  1:1228523/1‑93 (MQ=255)

                                                       tACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGCTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACca                                  >  1:3410606/1‑93 (MQ=255)

                                                        aCTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcac                                 >  2:646420/1‑93 (MQ=255)

                                                        aCTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGGGACGGCACCGCTGCGTAATCCCCCGACACCTGTAACcac                                 >  2:5801016/1‑93 (MQ=255)

                                                          tGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacac                               >  2:5118268/1‑93 (MQ=255)

                                                          tGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacac                               <  2:2113411/93‑1 (MQ=255)

                                                            cgcgGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAg                             >  2:1376764/1‑93 (MQ=255)

                                                            cgcgGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAg                             <  2:3194496/93‑1 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAgg                            <  2:590256/93‑1 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAgg                            >  2:4690584/1‑93 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAgg                            <  2:4056825/93‑1 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAgg                            >  2:773163/1‑93 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAgg                            <  2:2918813/93‑1 (MQ=255)

                                                               ggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCg                          >  2:2254923/1‑93 (MQ=255)

                                                                   cATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGaaa                      <  1:1833182/93‑1 (MQ=255)

                                                                    aTCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAAc                     >  1:6067213/1‑93 (MQ=255)

                                                                     tCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAAcc                    <  1:3714653/93‑1 (MQ=255)

                                                                      cGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACca                   >  1:3329227/1‑93 (MQ=255)

                                                                      cGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACca                   <  1:4487059/93‑1 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 <  2:1053774/93‑1 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 >  2:5057805/1‑93 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 >  1:811225/1‑93 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 <  1:672365/93‑1 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 <  1:2811630/93‑1 (MQ=255)

                                                                         ccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATCCCTGCAACCCCACAGGCGGAAACcagca                <  1:2719258/93‑1 (MQ=255)

                                                                         ccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagca                <  1:3550704/93‑1 (MQ=255)

                                                                         ccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagca                >  1:5509415/1‑93 (MQ=255)

                                                                          cATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagcag               >  1:5148066/1‑93 (MQ=255)

                                                                          cATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagcag               <  1:1243160/93‑1 (MQ=255)

                                                                          cATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagcag               >  1:488610/1‑93 (MQ=255)

                                                                           aTCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCagcagc              >  2:5590243/1‑93 (MQ=255)

                                                                              aGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCaaa           <  2:4170707/93‑1 (MQ=255)

                                                                              aGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCaaa           <  1:4514619/93‑1 (MQ=255)

                                                                               gcgcGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAAc          <  1:2495321/93‑1 (MQ=255)

                                                                                   ggcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGt      <  1:5976704/93‑1 (MQ=255)

                                                                                   ggcgAATGCCTGCGAGGTTTCCAGCGTAACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGt      >  1:960349/1‑93 (MQ=255)

                                                                                    gcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtg     >  2:2438414/1‑93 (MQ=255)

                                                                                    gcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtg     >  1:5254731/1‑93 (MQ=255)

                                                                                    gcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtg     <  1:392909/93‑1 (MQ=255)

                                                                                     cgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtga    <  2:285162/93‑1 (MQ=255)

                                                                                     cgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtga    >  1:3094724/1‑93 (MQ=255)

                                                                                     ccaaTTCCCGCGAAGTTTCCCGCGTGAACGCACCCGTGCCTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtga    <  2:3128919/91‑1 (MQ=255)

                                                                                       aaTGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtgatg  >  2:880320/1‑93 (MQ=255)

ATCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGTGATG  >  NC_000913/620206‑620385

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