BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,046,045	T→C	V339V (GTA→GTG)	gcvP ←	glycine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,046,045	0	T	C	99.6%	721.7 / ‑5.0	221	V339V (GTA→GTG)	gcvP	glycine decarboxylase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (91/129); total (92/129)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 228 total aligned reads displayed.

ATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCC  >  NC_000913/3045955‑3046134

atCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCAc                                                                                         >  2:1400066/1‑93 (MQ=255)

 acGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCAcc                                                                                        <  2:563268/92‑1 (MQ=255)

  cGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCt                                                                                       <  2:4095943/93‑1 (MQ=255)

   ggTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTg                                                                                      >  2:81379/1‑93 (MQ=255)

   ggTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTg                                                                                      >  2:6548108/1‑93 (MQ=255)

   ggTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTg                                                                                      <  1:644272/93‑1 (MQ=255)

   ggTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTg                                                                                      >  2:5711160/1‑93 (MQ=255)

      cAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGa                                                                                   >  2:3556895/1‑93 (MQ=255)

      cAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGa                                                                                   <  1:3356813/93‑1 (MQ=255)

       aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCACCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGaa                                                                                  <  2:6315074/93‑1 (MQ=255)

       aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGaa                                                                                  <  2:1210821/93‑1 (MQ=255)

       aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGaa                                                                                  >  2:3765670/1‑93 (MQ=255)

       aGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGaa                                                                                  <  1:2726512/93‑1 (MQ=255)

        gACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAg                                                                                 >  2:4498958/1‑93 (MQ=255)

         aCGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGt                                                                                <  2:6056741/93‑1 (MQ=255)

            gTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACa                                                                             <  2:5894679/93‑1 (MQ=255)

             tGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACaa                                                                            <  1:6301709/93‑1 (MQ=255)

              gAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACaaa                                                                           >  1:256385/1‑93 (MQ=255)

              gAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACaaa                                                                           <  1:6337528/93‑1 (MQ=255)

               aaTGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAAt                                                                          <  1:2538193/93‑1 (MQ=255)

                aTGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATg                                                                         <  2:2000923/93‑1 (MQ=255)

                 tGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGt                                                                        >  2:1581278/1‑93 (MQ=255)

                 tGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGt                                                                        >  2:2468364/1‑93 (MQ=255)

                 tGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGt                                                                        <  1:5311356/93‑1 (MQ=255)

                   cGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTg                                                                      >  2:1938923/1‑93 (MQ=255)

                   cGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTg                                                                      >  1:4608853/1‑93 (MQ=255)

                    ggTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTgg                                                                     <  2:6380767/93‑1 (MQ=255)

                     gTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGa                                                                    <  2:3216572/93‑1 (MQ=255)

                       tAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGt                                                                  <  1:168116/93‑1 (MQ=255)

                        aGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGtt                                                                 <  2:614601/93‑1 (MQ=255)

                        aGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGtt                                                                 <  2:2695713/93‑1 (MQ=255)

                        aGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGtt                                                                 <  2:3868113/93‑1 (MQ=255)

                         gCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTc                                                                >  1:451308/1‑93 (MQ=255)

                         gCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTc                                                                <  1:4197656/93‑1 (MQ=255)

                          cGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCg                                                               >  1:3774566/1‑93 (MQ=255)

                           gATACGTTTCAGGCCACCCGGGCCGTGATAAGCGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGc                                                              <  2:1442078/93‑1 (MQ=255)

                           gATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGc                                                              >  1:5730998/1‑93 (MQ=255)

                           gATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGc                                                              >  2:2144071/1‑93 (MQ=255)

                            aTACGTTTCAGGCCACCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                             <  1:3900413/93‑1 (MQ=255)

                            aTACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                             <  1:6282405/93‑1 (MQ=255)

                            aTACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                             >  2:2806526/1‑93 (MQ=255)

                            aTACGTTTCAGGCCAACCGGGCCGAGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCt                                                             <  1:3433146/93‑1 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGCGTTCGCtt                                                            <  1:6573059/93‑1 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                            <  2:4302830/93‑1 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                            >  1:95898/1‑93 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                            <  2:1769927/93‑1 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                            <  1:3749874/93‑1 (MQ=255)

                             tACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCtt                                                            <  2:1128947/93‑1 (MQ=255)

                               cGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTtc                                                          <  2:8661/93‑1 (MQ=255)

                                gTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTtct                                                         >  2:4789417/1‑93 (MQ=255)

                                   tCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAc                                                      <  1:3724500/93‑1 (MQ=255)

                                    cAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcg                                                     <  1:1345873/93‑1 (MQ=255)

                                     aGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcgg                                                    >  1:5094571/1‑93 (MQ=255)

                                     aGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcgg                                                    >  2:4607318/1‑93 (MQ=255)

                                     aGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcgg                                                    >  2:3642801/1‑93 (MQ=255)

                                      ggCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcggc                                                   >  2:4264095/1‑93 (MQ=255)

                                      ggCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCAcggc                                                   <  2:3904816/93‑1 (MQ=255)

                                         cAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCTCCTGGGAAGTACAAAAGTTGGAGTTCGCTTTCTCACGGCGGa                                                >  2:6734753/1‑93 (MQ=255)

                                          aaCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGat                                               >  1:149195/1‑93 (MQ=255)

                                           aCCGGGCGGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGata                                              <  1:3959379/93‑1 (MQ=255)

                                           aCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCTTTTTCTCACGGCGGata                                              <  2:3999159/93‑1 (MQ=255)

                                           aCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGata                                              >  2:2229812/1‑93 (MQ=255)

                                           aCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGata                                              >  1:668568/1‑93 (MQ=255)

                                           aCCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCAACTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGAGGata                                              >  1:2981329/1‑93 (MQ=255)

                                            ccGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGatat                                             <  1:6240828/93‑1 (MQ=255)

                                               ggCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGtt                                          >  1:2260765/1‑93 (MQ=255)

                                                   gTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTc                                      >  2:101155/1‑93 (MQ=255)

                                                    tGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATTGTGCTcg                                     >  2:4674797/1‑93 (MQ=255)

                                                    tGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTcg                                     >  2:2493316/1‑93 (MQ=255)

                                                       tAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGa                                  >  2:3140857/1‑93 (MQ=255)

                                                         aaCGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTTGTCGCGAgc                                >  1:1606662/1‑92 (MQ=255)

                                                         aaCGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGt                                >  1:1326554/1‑93 (MQ=255)

                                                          aCGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTc                               <  2:5147034/93‑1 (MQ=255)

                                                          aCGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTc                               <  1:5396912/93‑1 (MQ=255)

                                                           cGGCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCt                              <  2:208166/93‑1 (MQ=255)

                                                            ggCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGTATATGTTGCTCGCGAGTCTg                             >  1:6652000/1‑93 (MQ=255)

                                                             gCATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGGTGCTCGCGGGTCTGc                            >  2:4229222/1‑93 (MQ=255)

                                                              cATACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCa                           >  2:5593576/1‑93 (MQ=255)

                                                               aTACAGGCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCAt                          <  1:1063660/93‑1 (MQ=255)

                                                                    ggCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCa                     >  2:5327149/1‑93 (MQ=255)

                                                                    ggCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCa                     >  2:2372531/1‑93 (MQ=255)

                                                                     gCTGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCAt                    <  1:5625636/93‑1 (MQ=255)

                                                                       tGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgc                  >  1:5165930/1‑93 (MQ=255)

                                                                       tGGCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgc                  >  2:2056128/1‑93 (MQ=255)

                                                                        ggCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgcg                 >  2:5174391/1‑93 (MQ=255)

                                                                         gCGATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATgcgc                <  2:4732350/93‑1 (MQ=255)

                                                                           gATGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAg              >  2:6363367/1‑93 (MQ=255)

                                                                            aTGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgc             <  1:6507357/93‑1 (MQ=255)

                                                                             tGTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcg            >  2:6618181/1‑93 (MQ=255)

                                                                              gTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc           >  1:2437523/1‑93 (MQ=255)

                                                                              gTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc           <  2:2360967/93‑1 (MQ=255)

                                                                              gTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc           >  1:1234065/1‑93 (MQ=255)

                                                                              gTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc           >  2:3066790/1‑93 (MQ=255)

                                                                              gTTTGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAgcgc           <  1:3571052/93‑1 (MQ=255)

                                                                                 tGCCAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGt        >  1:3496782/1‑93 (MQ=255)

                                                                                   ccAGCAGCACCTGGTAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAt      <  1:2022713/93‑1 (MQ=255)

                                                                                   ccAGCAGCACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTAt      <  1:1818114/93‑1 (MQ=255)

                                                                                     agcagcACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTg    <  1:6377938/93‑1 (MQ=255)

                                                                                     agcagcACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTg    >  1:90330/1‑93 (MQ=255)

                                                                                     agcagcACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTg    >  2:1052827/1‑93 (MQ=255)

ATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCAGTACCTGGGAAGTACAAATGTTGGAGTTCGCTTTCTCACGGCGGATATGTTGCTCGCGAGTCTGCATCGCCATGCGCAGCGCGGTATTGCC  >  NC_000913/3045955‑3046134

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