BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,303	Δ1 bp	coding (3819/4281 nt)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,303	0	A	.	81.2%	617.5 / 177.0	245	coding (3819/4281 nt)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base A (2/1); major base . (125/70); minor base G (20/25); total (147/98)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-09

Only 100 of 276 total aligned reads displayed.

TTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAA  >  NC_000913/526222‑526392

tttaaatccactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTgtggtg                                                                                 >  2:5568576‑M2/79‑93 (MQ=255)

tttaaatccactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTgtggtg                                                                                 >  1:4880198‑M2/79‑93 (MQ=255)

  tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGT‑‑ACGG‑GGGGGAGTTTGTGGTGtt                                                                               <  2:4870038/93‑1 (MQ=255)

    aatccactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCt                                                                             >  2:9231514‑M2/75‑93 (MQ=255)

     atccactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCtt                                                                            >  1:6926885‑M2/74‑93 (MQ=255)

       ccactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTc                                                                          >  2:3570313‑M2/72‑93 (MQ=255)

       ccacaacaacaaattgaccctatgggattattgcagactggggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTc                                                                          <  2:2769824‑M2/22‑1 (MQ=255)

        cactacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCa                                                                         >  2:7278096‑M2/71‑93 (MQ=255)

        cACTACAAAAAATTGCCCTTATGGGTTTATTGCAGACTGGGGATGATGCAAGTTGTGGAGCAGGT‑‑CGGG‑GGGGGGGTTTGTGGTGTTCTTTCa                                                                         <  2:8368926/93‑1 (MQ=255)

         actacaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCAc                                                                        >  1:10025382‑M2/70‑93 (MQ=255)

          ctacaacaaattgaccctatgggattattgcagacttgggatgatgccaggtctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACg                                                                       <  2:6998951‑M2/25‑1 (MQ=255)

           tacaacaaattgaccctatgggattattgcagacttgggatgataccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGt                                                                      >  2:6373163‑M2/68‑93 (MQ=255)

            acaacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGta                                                                     >  1:4120207‑M2/67‑93 (MQ=255)

             caacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑GGTTGTGGTGTTCTTTCACGta                                                                     >  2:3795425‑M2/66‑93 (MQ=255)

             caacaaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGtat                                                                    >  1:9614641‑M2/66‑93 (MQ=255)

              aacaaATTGACCTTATGGGATTATTGCAGACTTGGTATGATCCCAGATCTGGAGCATGT‑‑CCGG‑GGGGGGGTTTGTGGTGTTCTTTCACGtata                                                                   <  2:7468873/93‑1 (MQ=255)

              aacaaATTGACCCTATGGGATTATTGCAGACTTTGGATGATGCCAGATCTGGAGCATGT‑‑AC‑GAGGGGGGGTTTTTGGTGTTCTGTCACGtata                                                                   >  1:4250559/1‑93 (MQ=255)

                caaattgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTTTTTCAAGTataat                                                                 >  2:2801555‑M2/63‑93 (MQ=255)

                  acttgcccctatgggattattgcagacttgggatgatgccagcgctggagcatgtccgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAg                                                               <  2:2334602‑M2/33‑1 (MQ=255)

                  aaTTGCCCCTATGGGATTTTGGCAGACTTGGGGAGATGCCAGATGTGGGGGATGTGGGG‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTataata                                                                <  1:5787013/93‑1 (MQ=255)

                   attgaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAgg                                                              >  2:6153155‑M2/60‑93 (MQ=255)

                      gaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑G‑TTGTGGTGTTCTTTCACGTATAATAGGAcc                                                           >  2:4542705‑M2/57‑93 (MQ=255)

                      gaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGAcc                                                           <  1:910205‑M2/37‑1 (MQ=255)

                      gaccctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGAcc                                                           <  2:5297383‑M2/37‑1 (MQ=255)

                       accctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑G‑TTGTGGTGTTTTTTCACGTATAATAGGACCa                                                          >  2:2402793‑M2/56‑93 (MQ=255)

                         cctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                                        >  1:3648831‑M2/54‑93 (MQ=255)

                          ctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAaga                                                       >  1:3201204‑M2/53‑92 (MQ=255)

                          ctatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                                       >  1:83992‑M2/53‑93 (MQ=255)

                          ctatgggattatggcggacttgggatgatgcaagatctggagcatgtccgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                                       <  1:5876400‑M2/41‑1 (MQ=255)

                           tatgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                      >  2:6717048‑M2/52‑93 (MQ=255)

                             tgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑GTTTGTGGTGTTCTTTCACGTTTAATAGGACCAAGTaa                                                     >  1:3306453‑M2/50‑93 (MQ=255)

                             tgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑G‑TTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                                    >  2:3121368‑M2/50‑93 (MQ=255)

                             tgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTGTAATAGGACCAAGTaaa                                                    >  2:8647793‑M2/50‑93 (MQ=255)

                             tgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                                    >  1:9561004‑M2/50‑93 (MQ=255)

                             tgggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGGGTTCTTTCACGTATAATAGGACCAAGTaaa                                                    >  1:8744976‑M2/50‑93 (MQ=255)

                              gggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAt                                                   <  1:6693831‑M2/45‑1 (MQ=255)

                               gtattattgcagatttgggatgatgacagatctggagcatgtgcgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                                  <  1:237516‑M2/46‑1 (MQ=255)

                               ggattattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                                  >  2:4715002‑M2/48‑93 (MQ=255)

                                  ttattgcaggcttgggatgatgccagatctggagcatttacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGa                                               <  2:8206070‑M2/49‑1 (MQ=255)

                                  ttattgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGGGGGTTTGTGGTGTTTTTTTACGTTTATAAGGACCCAGTAAAttt                                                 >  2:189873‑M2/45‑93 (MQ=255)

                                    attgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                             <  1:2995415‑M2/51‑1 (MQ=255)

                                    attgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                             <  2:584701‑M2/51‑1 (MQ=255)

                                    attgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                             <  2:8495944‑M2/51‑1 (MQ=255)

                                     ttgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                            >  1:7122643‑M2/42‑93 (MQ=255)

                                     ttgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                            >  2:1026143‑M2/42‑93 (MQ=255)

                                     ttgcagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                            >  1:2545828‑M2/42‑93 (MQ=255)

                                     ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                                           >  1:511525/1‑93 (MQ=255)

                                        gagacttgggatggtgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                         <  1:6171496‑M2/55‑1 (MQ=255)

                                        cagacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                         <  1:6052060‑M2/55‑1 (MQ=255)

                                         agactggggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                                        <  1:7828‑M2/56‑1 (MQ=255)

                                          gacttgggatgatgccagatctggagcatgtacgggG‑GGGGG‑GTTTGGGGGGTTCTTTCAAGTATAAAAGGACCAAGTAAATTTGATAGTacg                                        >  2:3441150‑M2/37‑92 (MQ=255)

                                          gacttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                                       >  1:5832963‑M2/37‑93 (MQ=255)

                                           agttgggatgatgccagatctggagcatgtgcgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGc                                      <  1:3801111‑M2/58‑1 (MQ=255)

                                            ctggggatgatgccagagttggagcatgtgcgggG‑GGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGc                                      <  2:6852475‑M2/59‑1 (MQ=255)

                                            cTTGGGATGATGCCAGATCTGGAGGATGT‑‑AC‑G‑GGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                                   <  1:1742110/93‑1 (MQ=255)

                                            cTTGGGATGATGCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                                    >  2:9818523/1‑93 (MQ=255)

                                            cTTGGGATGATGCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                                    <  1:4475990/93‑1 (MQ=255)

                                             ttgggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                                    >  1:2146807‑M2/34‑93 (MQ=255)

                                             ttgggatgatgccagatcgggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                                    <  1:144171‑M2/60‑1 (MQ=255)

                                                ggatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                                 >  1:1351507‑M2/31‑93 (MQ=255)

                                                 gatgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                                >  1:8240376‑M2/30‑93 (MQ=255)

                                                   tgatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                              >  1:9953430‑M2/28‑93 (MQ=255)

                                                    tgctccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctgc                             >  2:4780114‑M2/27‑93 (MQ=255)

                                                    gatgccagatctggagcatgtacgggG‑GGGGG‑G‑TTGTGGGGTTCTTTCACGTTTAATAGGACCAAGTAAATTTGATAGTACTGCAGatgctgg                             >  2:5125003‑M2/27‑92 (MQ=255)

                                                    gatgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctgc                             >  1:9243758‑M2/27‑93 (MQ=255)

                                                    gatgCCAGATCTGGAGCATGT‑‑AC‑GCGGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctgc                             >  2:9651888/1‑93 (MQ=255)

                                                    gatgCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                            <  1:5630350/93‑1 (MQ=255)

                                                     atgccagatctggagcatgtacgggG‑GGGGG‑GTTTGGGGGGTTCTTTCACGGATAATAGAACCAAGTAAATTTGATAGTACTGCAGAtgctgc                             >  2:4266160‑M2/26‑93 (MQ=255)

                                                     atgccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                            <  1:1841552‑M2/68‑1 (MQ=255)

                                                        ccagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                         >  2:2670233‑M2/23‑93 (MQ=255)

                                                         cagatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        >  1:7397791‑M2/22‑93 (MQ=255)

                                                         cagatcgggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        <  2:7555512‑M2/72‑1 (MQ=255)

                                                          agatctggagcatgtacgggG‑GGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                        >  1:7421638‑M2/21‑93 (MQ=255)

                                                          agatctggagcatgtacgggG‑GGGGG‑GTTTGTGGGGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTACAGATGCTGcgttcg                        >  1:6348164‑M2/21‑91 (MQ=255)

                                                          agatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa                       <  1:9816112‑M2/73‑1 (MQ=255)

                                                          agatctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa                       >  2:2590490‑M2/21‑93 (MQ=255)

                                                           gttctggagcatgtgcgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt                      <  2:7214006‑M2/74‑1 (MQ=255)

                                                           agatcgggagcatgtagggG‑GGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa                       <  1:4749975‑M2/74‑1 (MQ=255)

                                                             tctggagcatgtacgggG‑GGGGG‑GTTTGTGGTGTTCTTTCACGGATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTagagg                     >  2:2517976‑M2/18‑91 (MQ=255)

                                                             tctggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGc                    >  2:6069843‑M2/18‑93 (MQ=255)

                                                              cTGGAGCATGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtt                  <  1:5361089/93‑1 (MQ=255)

                                                               tggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtt                  >  2:1928004‑M2/16‑93 (MQ=255)

                                                               tggagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtt                  <  2:1806481‑M2/78‑1 (MQ=255)

                                                               tGGAGCATGT‑‑CC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt                 <  2:9366906/93‑1 (MQ=255)

                                                                 gagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTg                >  2:1686573‑M2/14‑93 (MQ=255)

                                                                 gagcatgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTg                >  2:7729240‑M2/14‑93 (MQ=255)

                                                                     atgtacgggG‑GGGG‑‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaag            >  1:3501540‑M2/10‑93 (MQ=255)

                                                                      tGT‑‑AC‑G‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaa          >  2:4015327/1‑93 (MQ=255)

                                                                        gttgggG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga           <  2:4450521‑M2/87‑1 (MQ=255)

                                                                         ttcggG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaa          <  2:2498742‑M2/88‑1 (MQ=255)

                                                                          tgcgG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:8956983‑M2/89‑1 (MQ=255)

                                                                          tacgG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         >  2:5462112‑M2/5‑93 (MQ=255)

                                                                          tacgG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:6542042‑M2/89‑1 (MQ=255)

                                                                          tacgG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:8781242‑M2/89‑1 (MQ=255)

                                                                          tacgG‑GGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaagaaa         <  2:6311495‑M2/89‑1 (MQ=255)

                                                                              g‑gggggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGcg     >  1:1117357/1‑92 (MQ=255)

                                                                              g‑gggggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCa     <  2:3877487/93‑1 (MQ=255)

                                                                              g‑gggggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCa     <  1:4652771/93‑1 (MQ=255)

                                                                                 ggggggTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGa   >  2:1176264/1‑93 (MQ=255)

                                                                                  ggggGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGaa  >  1:2722498/1‑93 (MQ=255)

TTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGT‑‑AC‑G‑GGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAA  >  NC_000913/526222‑526392

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