BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	707,634	A→G	A26A (GCA→GCG)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	707,634	0	A	G	99.6%	833.9 / ‑5.9	266	A26A (GCA→GCG)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (134/131); total (134/132)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 274 total aligned reads displayed.

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         <  1:4985752/93‑1 (MQ=255)

 gagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATg                                                                                        >  1:507104/1‑93 (MQ=255)

 gagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATg                                                                                        <  2:2116812/93‑1 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGGTGGCTTTCGGCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGcc                                                                                      <  1:18334/93‑1 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGcc                                                                                      <  2:6507114/93‑1 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGcc                                                                                      >  1:7694697/1‑93 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGcc                                                                                      >  2:5607592/1‑93 (MQ=255)

    gATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCg                                                                                     >  1:4857526/1‑93 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGa                                                                                   <  1:3920206/93‑1 (MQ=255)

       tAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAg                                                                                  <  2:3732894/93‑1 (MQ=255)

        aaCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAgg                                                                                 >  1:2057937/1‑93 (MQ=255)

         aCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGc                                                                                <  1:5653772/93‑1 (MQ=255)

          cccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCt                                                                               <  2:4550311/93‑1 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              >  1:5827493/1‑93 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              >  1:203670/1‑93 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              >  2:1724244/1‑93 (MQ=255)

            cATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcg                                                                             >  2:3233785/1‑93 (MQ=255)

            cATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcg                                                                             >  2:4044010/1‑93 (MQ=255)

            cATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGATGGCTATCGACGGAGTTACAGCAATGTCGGGCTTTATGGGGATGCCGGAGGCTcc                                                                             >  2:1124117/1‑92 (MQ=255)

                gggaCGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGcc                                                                         <  1:3744757/90‑1 (MQ=255)

                cGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGcc                                                                         >  2:3472121/1‑93 (MQ=255)

                cGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGcc                                                                         >  2:2585567/1‑93 (MQ=255)

                 gTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCg                                                                        >  2:668346/1‑93 (MQ=255)

                 gTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCg                                                                        >  1:717006/1‑93 (MQ=255)

                 gTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGAGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCg                                                                        >  2:6559307/1‑93 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGa                                                                       <  2:2452470/93‑1 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCGTTATGCAGATGCTGGGGGCTCGCGTCgg                                                                       >  2:4757188/1‑92 (MQ=255)

                  gaCGTTTAGGGGCAAACGTAGGACGCTGGCGCTGGCTAACGCCGGTGTTACAGCAATGTCGGGCGTTATGGCGATGCCGGAGGCTCGCGCCGa                                                                       >  2:6032958/2‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACCGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGCAg                                                                     >  1:3303072/1‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAg                                                                     >  1:6109895/1‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAg                                                                     <  2:6844039/93‑1 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    >  1:1093337/1‑93 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    >  1:3651356/1‑93 (MQ=255)

                       ttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAAGAt                                                                  >  2:5732095/1‑93 (MQ=255)

                        tAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGagt                                                                 >  2:630106/1‑91 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                >  2:7658544/1‑93 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                >  2:2951013/1‑93 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                >  1:2443538/1‑93 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                >  1:6493796/1‑93 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               >  1:3450898/1‑93 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               >  2:301409/1‑93 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  1:2910301/93‑1 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  1:7496446/93‑1 (MQ=255)

                             gCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcg                                                            <  2:1238902/93‑1 (MQ=255)

                              cAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATcga                                                           <  2:6207822/93‑1 (MQ=255)

                                 aCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgacga                                                        <  1:2753954/93‑1 (MQ=255)

                                  cGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAt                                                       >  2:6310281/1‑93 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      <  2:1188937/93‑1 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      <  1:676502/93‑1 (MQ=255)

                                    tagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTc                                                     <  2:3904775/93‑1 (MQ=255)

                                      gtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCaa                                                   >  1:7437460/1‑93 (MQ=255)

                                       taCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAAc                                                  <  1:4865423/93‑1 (MQ=255)

                                           ctggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTa                                              >  1:7811989/1‑93 (MQ=255)

                                             ggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAc                                            <  2:5468017/93‑1 (MQ=255)

                                             ggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAc                                            <  2:7229363/93‑1 (MQ=255)

                                             ggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAc                                            <  1:1529153/93‑1 (MQ=255)

                                             ggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAc                                            <  2:87048/93‑1 (MQ=255)

                                              gcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAcc                                           <  1:1154669/93‑1 (MQ=255)

                                                gctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcgg                                         >  1:58819/1‑93 (MQ=255)

                                                 ctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggc                                        <  2:4560522/93‑1 (MQ=255)

                                                 ctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCCCCGAAGGATTCATCGACGATTCAACCTTAACcggc                                        <  1:3234562/93‑1 (MQ=255)

                                                  tggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcg                                       <  1:5762356/93‑1 (MQ=255)

                                                  tggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcg                                       <  1:1064795/93‑1 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      <  1:6655439/93‑1 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      >  1:4434564/1‑93 (MQ=255)

                                                    gcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGt                                     >  2:5355770/1‑93 (MQ=255)

                                                     cTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTa                                    <  2:5259218/93‑1 (MQ=255)

                                                     cTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTa                                    <  2:222150/93‑1 (MQ=255)

                                                      tATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTAt                                   >  2:5785978/1‑93 (MQ=255)

                                                       aTCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATc                                  <  2:2637262/93‑1 (MQ=255)

                                                       aTCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATc                                  <  2:7870046/93‑1 (MQ=255)

                                                        tCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCt                                 <  1:4513717/93‑1 (MQ=255)

                                                        tCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCt                                 >  1:1245108/1‑93 (MQ=255)

                                                          gCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAt                               <  1:2312606/93‑1 (MQ=255)

                                                              gtgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACt                           <  1:1340519/93‑1 (MQ=255)

                                                              gtgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACt                           <  2:5857985/93‑1 (MQ=255)

                                                               tgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTg                          >  1:947954/1‑93 (MQ=255)

                                                               tgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTg                          >  1:2683479/1‑93 (MQ=255)

                                                                gtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGGGCCGAAGGATTCATCGACGAATCAACCTTAACCGGCGGGATTTATTACTgg                         >  2:2277945/1‑93 (MQ=255)

                                                                gtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTgg                         <  2:1500159/93‑1 (MQ=255)

                                                                  tACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCa                       <  1:5303941/93‑1 (MQ=255)

                                                                  tACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCa                       <  1:3999938/93‑1 (MQ=255)

                                                                  tACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCa                       <  1:4363291/93‑1 (MQ=255)

                                                                    cAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGc                     >  1:5315722/1‑93 (MQ=255)

                                                                     aGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCg                    >  1:6931023/1‑93 (MQ=255)

                                                                     aGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCg                    <  2:7103303/93‑1 (MQ=255)

                                                                     aGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCg                    >  2:1101748/1‑93 (MQ=255)

                                                                        aaTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTTGCAGCGTGa                 >  2:1314815/1‑93 (MQ=255)

                                                                         aTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                >  1:4614733/1‑93 (MQ=255)

                                                                         aTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                <  1:7549970/93‑1 (MQ=255)

                                                                          tGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGTAGCCTGAAc               >  2:2039946/1‑93 (MQ=255)

                                                                          tGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAc               <  2:7451902/93‑1 (MQ=255)

                                                                            tCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcgc             <  2:2400371/93‑1 (MQ=255)

                                                                               ggCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAc          <  1:635455/93‑1 (MQ=255)

                                                                               ggCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAc          <  2:4244114/93‑1 (MQ=255)

                                                                                gCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGGTTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAcc         >  1:6798351/1‑93 (MQ=255)

                                                                                gCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAcc         >  1:6905797/1‑93 (MQ=255)

                                                                                    tATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTaa     >  1:7347093/1‑93 (MQ=255)

                                                                                     aTGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTaaa    <  1:699507/93‑1 (MQ=255)

                                                                                      tGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   >  1:2996581/1‑93 (MQ=255)

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

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