BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,918	A→G	C87R (TGC→CGC)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,918	0	A	G	100.0%	410.7 / NA	116	C87R (TGC→CGC)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (53/63); total (53/63)

Only 100 of 116 total aligned reads displayed.

ATACTGCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATT  >  NC_000913/4441772‑4442065

aTACTGCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACt                                                                                                                                                 >  2:689679/1‑151 (MQ=255)

  aCTGCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTc                                                                                                                                               >  2:1121653/1‑151 (MQ=255)

  aCTGCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTc                                                                                                                                               >  2:458297/1‑151 (MQ=255)

     gCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCg                                                                                                                                            <  4:298361/151‑1 (MQ=255)

         gCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATaa                                                                                                                                        <  1:45608/151‑1 (MQ=255)

          ccGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATa                                                                                                                                         >  1:291435/1‑149 (MQ=255)

           cGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGt                                                                                                                                      >  3:998497/1‑151 (MQ=255)

           cGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGt                                                                                                                                      >  2:35186/1‑151 (MQ=255)

            gTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGTATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCTCCGGAGCGCACTTCCCGATAAGTc                                                                                                                                     <  3:869855/151‑1 (MQ=255)

            gTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTc                                                                                                                                     <  2:132900/151‑1 (MQ=255)

               gTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTTTCACCGGAGCGCACTTCCCGATAAGTCGGa                                                                                                                                  <  3:1203480/151‑1 (MQ=255)

               gTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCgg                                                                                                                                   >  2:338096/1‑150 (MQ=255)

               gTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGa                                                                                                                                  <  3:427181/151‑1 (MQ=255)

                 cATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGAtt                                                                                                                                <  4:409143/151‑1 (MQ=255)

                       ccTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGt                                                                                                                          <  3:504631/151‑1 (MQ=255)

                          gACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGta                                                                                                                         <  3:686158/149‑1 (MQ=255)

                           aCGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGtat                                                                                                                        >  3:247384/1‑149 (MQ=255)

                           aCGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGtat                                                                                                                        >  1:362906/1‑149 (MQ=255)

                           aCGCATTCCCTGGGCGGGATAGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAg                                                                                                                      <  3:862823/151‑1 (MQ=255)

                            cGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgc                                                                                                                     <  2:438510/151‑1 (MQ=255)

                            cGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgc                                                                                                                     >  4:417208/1‑151 (MQ=255)

                             gCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgcc                                                                                                                    >  3:1152546/1‑151 (MQ=255)

                             gCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgcc                                                                                                                    >  3:1194266/1‑151 (MQ=255)

                               aTTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgccgc                                                                                                                  <  4:408457/151‑1 (MQ=255)

                                ttCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAgccgcc                                                                                                                 >  3:890419/1‑151 (MQ=255)

                                 tCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGccgccg                                                                                                                >  3:86261/1‑151 (MQ=255)

                                  cccTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGccgccg                                                                                                                <  2:795385/150‑1 (MQ=255)

                                     tGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGtat                                                                                                            >  3:531087/1‑151 (MQ=255)

                                      gggCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGtata                                                                                                           <  2:1273222/151‑1 (MQ=255)

                                      gggCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGta                                                                                                             <  3:257318/149‑1 (MQ=255)

                                       ggCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAg                                                                                                          <  2:921319/151‑1 (MQ=255)

                                         cGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGcc                                                                                                        >  2:531465/1‑151 (MQ=255)

                                         cGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGcc                                                                                                        >  3:147513/1‑151 (MQ=255)

                                            gATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgc                                                                                                     <  3:543823/151‑1 (MQ=255)

                                            gATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgc                                                                                                     >  4:309979/1‑151 (MQ=255)

                                               cGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgcggc                                                                                                  >  2:1288588/1‑151 (MQ=255)

                                               cGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgcggc                                                                                                  <  4:422911/151‑1 (MQ=255)

                                                gTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgcggcg                                                                                                 >  2:221442/1‑151 (MQ=255)

                                                 tGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGcggcgg                                                                                                >  3:953291/1‑151 (MQ=255)

                                                 tGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGcggcgg                                                                                                >  3:1165785/1‑151 (MQ=255)

                                                  gATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGACGAGGGAAGGATCGTAACCAATGCGTACCGCTTCGGCATGACCCGTACCACCGGAGCGCACTGCCCGATAAGTGGGATTTGGCGTATAGCCGCCGGTATAGCCTgcggcg                                                                                                 <  4:195074/149‑1 (MQ=255)

                                                     tctcCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTgcggcg                                                                                                 >  3:9933/1‑146 (MQ=255)

                                                         ccAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTa                                                                                         >  3:492093/1‑150 (MQ=255)

                                                          cAAAATACCTGTAGCAACTGCTCACAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTGCGGCATGGCCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTaaa                                                                                       <  1:397512/151‑1 (MQ=255)

                                                            aaaTACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACg                                                                                     <  2:417440/151‑1 (MQ=255)

                                                             aaTACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGc                                                                                    <  3:555649/151‑1 (MQ=255)

                                                             aaTACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGc                                                                                    >  1:207463/1‑151 (MQ=255)

                                                              aTACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGcc                                                                                   >  2:927416/1‑151 (MQ=255)

                                                                 ccTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCggg                                                                                >  2:402671/1‑151 (MQ=255)

                                                                   tGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGt                                                                               >  1:124508/1‑150 (MQ=255)

                                                                    gTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGGAAACGCCGGGt                                                                               >  3:1138281/1‑149 (MQ=255)

                                                                     tAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAAc                                                                            >  2:588234/1‑151 (MQ=255)

                                                                      aGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTaa                                                                             >  4:71629/1‑149 (MQ=255)

                                                                        cAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTg                                                                          <  3:222285/150‑1 (MQ=255)

                                                                             gCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGa                                                                     <  3:552229/150‑1 (MQ=255)

                                                                             gCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGACAGa                                                                     <  2:1237396/150‑1 (MQ=255)

                                                                              cTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGaa                                                                    >  3:32296/1‑150 (MQ=255)

                                                                                cATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACa                                                                  <  2:623030/150‑1 (MQ=255)

                                                                                cATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGCCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAAc                                                                   <  2:167637/149‑1 (MQ=255)

                                                                                 aTAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAAATGCCAGAAc                                                                   >  1:414451/1‑148 (MQ=255)

                                                                                       gatgaCGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTcc                                                          <  3:401870/151‑1 (MQ=255)

                                                                                       gatgaCGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGGAACTGCCAGAAAAGACGCTcc                                                          >  3:79537/1‑151 (MQ=255)

                                                                                        atgaCGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGCATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTcc                                                          >  1:414395/1‑150 (MQ=255)

                                                                                         tgaCGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCcac                                                        >  2:341719/1‑151 (MQ=255)

                                                                                           aCGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCcaca                                                       <  1:176258/150‑1 (MQ=255)

                                                                                            cGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATATGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACAcc                                                     >  1:84204/1‑151 (MQ=255)

                                                                                            cGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCacac                                                      <  4:124507/150‑1 (MQ=255)

                                                                                             ggAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAccc                                                    >  4:135231/1‑151 (MQ=255)

                                                                                                aGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAg                                                 <  2:327762/151‑1 (MQ=255)

                                                                                                  gATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAAAGCCGGGTAAATCCCAGAACAGAAGCTCCACACCCCAGa                                                >  2:888538/1‑150 (MQ=255)

                                                                                                   aTCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGaaa                                              <  1:388872/151‑1 (MQ=255)

                                                                                                   aTCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGaaa                                              <  4:476422/151‑1 (MQ=255)

                                                                                                   aTCGTAAAAAATGCTTACCGTTTCGGCATTACCCGTATCACCGGAGCGCATTTCCCGATAAGTCGGATTTGCCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGaa                                               <  3:458332/150‑1 (MQ=255)

                                                                                                     cGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAaaca                                            <  2:32289/151‑1 (MQ=255)

                                                                                                            aaTGCGTACCGCTTCGGCAGGCCCCGTATCCCCGGAGCGCACTTCCCGATAAGCCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTACCTGCCAGAACAGACGCTCCACACCCCAGAAACACCCCATc                                     <  3:251410/151‑1 (MQ=255)

                                                                                                             aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATc                                     >  1:279784/1‑150 (MQ=255)

                                                                                                             aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCg                                    <  1:371082/151‑1 (MQ=255)

                                                                                                             aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCg                                    <  2:953957/151‑1 (MQ=255)

                                                                                                             aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAAAGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCg                                    >  2:1274358/1‑151 (MQ=255)

                                                                                                              tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCg                                    >  4:283490/1‑150 (MQ=255)

                                                                                                               gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCa                                  <  4:255897/151‑1 (MQ=255)

                                                                                                               gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCa                                  <  1:150230/151‑1 (MQ=255)

                                                                                                                 gcaCCGCTTCGGCATGACCCGTATCACCGGAGCGCACTCCCCGATAAGGCGGATTTGGGGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCaaa                                <  3:1118446/149‑1 (MQ=255)

                                                                                                                  tACCGCTTCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCaaa                                >  2:1290967/1‑150 (MQ=255)

                                                                                                                         tCGGCATGACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGGAAACGCCGGGTAAATGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAAAGGCa                         >  2:1103182/1‑150 (MQ=255)

                                                                                                                                gACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                  <  3:1288656/150‑1 (MQ=255)

                                                                                                                                gACCCGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                 <  2:1082614/151‑1 (MQ=255)

                                                                                                                                  cccGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                 <  3:35193/149‑1 (MQ=255)

                                                                                                                                    cGTATCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg             <  1:193599/151‑1 (MQ=255)

                                                                                                                                       aTCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc           <  3:531499/150‑1 (MQ=255)

                                                                                                                                        tCACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCCGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg             >  3:1158765/1‑147 (MQ=255)

                                                                                                                                         cACCGGAGCGCCCTCCCCGATAAGGCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAgg        <  3:369191/151‑1 (MQ=255)

                                                                                                                                         cACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc           <  3:689704/148‑1 (MQ=255)

                                                                                                                                         cACCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg         <  2:953253/150‑1 (MQ=255)

                                                                                                                                          aCCGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAgg        >  1:445877/1‑150 (MQ=255)

                                                                                                                                           ccGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt       <  4:231253/150‑1 (MQ=255)

                                                                                                                                            cGGAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTAc     >  2:169165/1‑151 (MQ=255)

                                                                                                                                             ggAGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACa    <  4:120144/151‑1 (MQ=255)

                                                                                                                                              gAGCGCACTTCCCGATAAGTCGGATTTGGCGTAAAGCCGCCGGTATACCCTGCGGCGGGGGGGTAAAAGCCGGGTAACAGCCAGAACAAACGCCCCACACCCCCTAAAAAACCCATTTCAACAACGGGAATCTACAGTCCATCacggaa     >  2:218123/1‑144 (MQ=255)

                                                                                                                                               aGCGCACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAtt  >  2:311190/1‑151 (MQ=255)

ATACTGCGTGCCGTGGTCATTGCCCTGACGCATTCCCTGGGCGGGATCGTGATTCTCCCAAAATACCTGTAGCAACTGCTCATAGCTGATGACGGAAGGATCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATT  >  NC_000913/4441772‑4442065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`