BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,123,585	+G	coding (953/2286 nt)	glcF ←	glycolate oxidase iron‑sulfur subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,123,582	1	.	G	94.0%	205.4 / 7.3	83	coding (956/2286 nt)	glcF	glycolate oxidase iron‑sulfur subunit
Reads supporting (aligned to +/- strand): ref base . (2/3); new base G (45/33); total (47/36)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑G‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAT  >  NC_000913/3123507‑3123652

gcgggggTAAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGc                                                                            <  2:4571597/73‑1 (MQ=255)

cccGGGGTTAACCCCGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCCGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGc                                                                            <  1:1181026/75‑1 (MQ=255)

aCCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGc                                                                            <  2:791083/76‑1 (MQ=255)

 ccGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g                                                                          <  2:4422980/76‑1 (MQ=255)

 ccGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g                                                                          >  1:3581506/1‑76 (MQ=255)

 ccGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g                                                                          >  2:1756529/1‑76 (MQ=255)

 ccGGGGATAAACACCCCCCAAGGGTGGAGCGGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g                                                                          <  1:45845/76‑1 (MQ=255)

  cGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑g                                                                        <  2:554103/76‑1 (MQ=255)

  cGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑g                                                                        <  2:2471038/76‑1 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       >  2:715400/1‑76 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       >  2:2314112/1‑76 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       <  2:2909410/76‑1 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       >  2:770954/1‑76 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       <  1:4121969/76‑1 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       >  2:3980997/1‑76 (MQ=255)

   ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCCGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gg                                                                       >  1:1851978/1‑76 (MQ=255)

   ggggTTAAACACGCCGCAAGGGACGAGCTGCTGTTTAAGCTGGTGGTGATTGCGGAATAAA‑CGAG‑CA‑CATAGg‑g‑gg                                                                       >  1:931244/1‑76 (MQ=255)

     ggTTAAACACGCCGCAGGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gggc                                                                     <  2:4218450/76‑3 (MQ=255)

     ggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gggc                                                                     >  1:4085664/1‑74 (MQ=255)

     ggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑g‑gggc                                                                     >  1:1270878/1‑74 (MQ=255)

      gTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGT‑g‑gggca                                                                    >  1:2469970/1‑73 (MQ=255)

          aaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GCTACA‑C‑GGC‑aaaagc                                                               >  2:4649971/1‑73 (MQ=255)

          aaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gc                                                               <  1:1514216/76‑1 (MQ=255)

          aaCACCCCGCAAGGGGCGAGCGGCTGTTTAAGGTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gc                                                               <  2:875299/76‑1 (MQ=255)

             ccgccgcAGGGGTCGAGCTGCTTTTTAAGTTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gccac                                                            <  1:3547733/75‑1 (MQ=255)

             acGCCGCAAGGGCCGAGCTGCTGTTTCAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gccac                                                            >  1:2756781/1‑76 (MQ=255)

              cgccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gccacc                                                           <  1:879831/76‑1 (MQ=255)

               gccgcAGGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gccacca                                                          <  1:4804911/76‑1 (MQ=255)

               gccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑Gccacca                                                          <  2:2543044/76‑1 (MQ=255)

                ccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCAt                                                         <  1:1432518/76‑1 (MQ=255)

                    agggggCGGGCTGCTGTTGAGGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTcc                                                     <  2:3225277/74‑1 (MQ=255)

                    aaGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGTGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGG‑CAA‑GCCACCATCTcc                                                     >  1:4857385/1‑76 (MQ=255)

                    aaGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTcc                                                     <  2:2631373/76‑1 (MQ=255)

                     aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCg                                                    >  1:4262715/1‑76 (MQ=255)

                     aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCg                                                    <  1:3882947/76‑1 (MQ=255)

                      gggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCgg                                                   >  1:3504086/1‑76 (MQ=255)

                       ggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGc                                                  >  2:4519784/1‑76 (MQ=255)

                       ggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGc                                                  >  1:1268032/1‑76 (MQ=255)

                        gtggAGCTGCTGTTTAGGGTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCa                                                 <  1:375069/73‑1 (MQ=255)

                         tCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCAc                                                >  2:2918189/1‑76 (MQ=255)

                         tCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCAc                                                <  1:1185945/76‑1 (MQ=255)

                          cGAGCTGCTGTTTAAGCTGCTGGTGATAACGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCCCCCGCAcc                                               >  2:479473/1‑75 (MQ=255)

                          cGAGCGGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACt                                               <  1:2832191/76‑1 (MQ=255)

                           gAGCTGCTGTTTAGGGTACCGGGGATAGCGGAAGAAG‑GGAG‑CC‑GATGGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTa                                              <  2:4522140/76‑1 (MQ=255)

                           gAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCCTCTCCGGCACTa                                              >  2:432660/1‑76 (MQ=255)

                           gAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCCCTa                                              >  2:4579034/1‑76 (MQ=255)

                           gAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCAcca                                              >  2:3995262/1‑74 (MQ=255)

                            aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCCACATCCCCCGCCCtca                                             >  2:4866647/1‑74 (MQ=255)

                            aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTaa                                             >  2:3452448/1‑76 (MQ=255)

                            aGCTGCTGTTTAAGCTGCTGGCGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCCTCTCCGGCACTaa                                             >  2:2631146/1‑76 (MQ=255)

                             gctgctTTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTaaa                                            <  1:2355117/76‑1 (MQ=255)

                             gctgctGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTaaa                                            <  1:2247983/76‑1 (MQ=255)

                                  tgTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGAGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCggg                                       >  1:1153207/1‑76 (MQ=255)

                                   gTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGt                                      >  1:732112/1‑76 (MQ=255)

                                   gTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGt                                      <  1:1958439/76‑1 (MQ=255)

                                    tttAAGCTGCTGTTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTc                                     <  1:4685481/76‑1 (MQ=255)

                                    tttAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑GCGGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTc                                     >  2:1851236/1‑76 (MQ=255)

                                     taaaGCTGCTGGTGAAGGGGGAAAAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCg                                    <  2:262851/74‑1 (MQ=255)

                                       aaGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCa                                  >  1:3410408/1‑76 (MQ=255)

                                       aaGCTGCGGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCa                                  <  1:770954/76‑1 (MQ=255)

                                        agatgcggGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                 <  2:4608057/73‑1 (MQ=255)

                                        aGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                 >  2:4424531/1‑76 (MQ=255)

                                        aGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                 <  1:3847899/76‑1 (MQ=255)

                                         ggtgcggGTGATAGCGGAATAGA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATg                                <  2:1930958/74‑1 (MQ=255)

                                         gctgctGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATg                                >  2:1849776/1‑76 (MQ=255)

                                          ctgctgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGa                               >  1:1548313/1‑76 (MQ=255)

                                           tgctgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGAc                              >  1:2836010/1‑76 (MQ=255)

                                           tgctgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGAc                              >  2:3544678/1‑76 (MQ=255)

                                            gcgggTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCCCTAAAGCGGGTCGCATGAcc                             <  2:2929968/73‑1 (MQ=255)

                                             tgggTGATAGCGGAATAAG‑GGCG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCCTCTCCGGCACTAAAGCGGGTCGCATGACCg                            <  2:2124579/74‑1 (MQ=255)

                                             ctgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                            >  1:3779757/1‑76 (MQ=255)

                                             ctgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                            <  1:4719366/76‑1 (MQ=255)

                                              tgGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGc                           >  1:771262/1‑76 (MQ=255)

                                              tgGTGATAGCGGAAAAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGc                           <  1:3757392/76‑1 (MQ=255)

                                                gTGATAGCGAAAAAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCa                         <  1:1349216/76‑1 (MQ=255)

                                                 tGATAGCGGAATAAA‑GGAG‑CC‑GATAGAGG‑GGC‑AAA‑GCCACCATCTCCGGCCCTAAAGCGGGTCGCATGACCCCCAg                        >  2:396939/1‑76 (MQ=255)

                                                   aTAGCGGAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGAACGCCAGCg                      >  1:794836/1‑76 (MQ=255)

                                                      ggggAAGAAA‑GGGG‑GCGGATAGG‑G‑GGCGAAA‑GCCGCCATCGCCGGCACTAAAGCGGGTCGCATGCCCGCCAGCGtt                    <  2:4113163/74‑1 (MQ=255)

                                                      gCGAAATGAA‑GGAG‑CG‑GATGGC‑G‑GGC‑AAACGCCACCACCCCCGGCACTAAAGCGGGTCGAAGGCCCGCCAGCGTTg                   <  1:1167752/76‑1 (MQ=255)

                                                      agggAAAAAAGGGAG‑CC‑GAGAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                    <  1:642187/74‑1 (MQ=255)

                                                        ggAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCCGCGTTGCg                 >  1:196015/1‑76 (MQ=255)

                                                        ggAATAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCg                 >  2:4363993/1‑76 (MQ=255)

                                                        ggAATAAA‑GGAG‑CC‑GAGAGCGG‑GGC‑AAA‑GCAACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGGTGCg                 >  1:2092508/1‑76 (MQ=255)

                                                          aaaaaa‑GGAGCCC‑GATGGC‑G‑GGGAAAA‑GCCCCCATCTCCGGCACAAAAGCGGGTCGCATGACCGCCAGCGTTGCgg                <  1:4026884/73‑1 (MQ=255)

                                                          aaTAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCggg               >  2:1695094/1‑76 (MQ=255)

                                                           aTAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgc              >  1:3067607/1‑76 (MQ=255)

                                                            tAAA‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCCTGCACGCCCGCGCTGCGGgcg             >  1:4806786/1‑76 (MQ=255)

                                                             aaa‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAg                                           >  1:4051263/1‑45 (MQ=255)

                                                             aaa‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAg                                           <  2:4051263/45‑1 (MQ=255)

                                                             aaa‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCACGACCGCCAGCGCTGCGGgcga            >  2:1118512/1‑76 (MQ=255)

                                                              ga‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCCGCGTTGCGGgcgag           >  1:4701433/2‑75 (MQ=255)

                                                              aa‑GGAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGCAGCATGACCGCCCGCGTTGCGGgcgag           >  1:4711785/1‑75 (MQ=255)

                                                                 ggAG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgc         >  1:3022204/1‑76 (MQ=255)

                                                                   aG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcga       >  1:3115583/1‑76 (MQ=255)

                                                                   aG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcga       >  1:2143361/1‑76 (MQ=255)

                                                                   aG‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcga       >  1:1175736/1‑76 (MQ=255)

                                                                    g‑CC‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGGTGCCGgcgatgcgat      >  1:3699294/1‑76 (MQ=255)

                                                                    accc‑GAGAGCGG‑GGA‑AAA‑CCCCCCATCCCGGGCACTAAGGCGGGCCGCATGACCGCCAGCGTTGCGGgcgatgcga       <  1:4118149/74‑1 (MQ=255)

                                                                       c‑GATAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGGTGCGGTgg    >  2:1203587/1‑76 (MQ=255)

                                                                          aTAGCGG‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAt  >  2:4005974/1‑76 (MQ=255)

ACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAA‑GGAG‑CC‑GATAGC‑G‑GGC‑AAA‑GCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAT  >  NC_000913/3123507‑3123652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 5 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`