BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	585,303	T→C	N111D (AAC→GAC)	ompT ←	DLP12 prophage; outer membrane protease VII; outer membrane protein 3b

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	585,303	0	T	C	100.0%	441.3 / NA	115	N111D (AAC→GAC)	ompT	DLP12 prophage; outer membrane protease VII; outer membrane protein 3b
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (48/67); total (48/67)

Only 100 of 115 total aligned reads displayed.

ATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCG  >  NC_000913/585157‑585448

ataACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTg                                                                                                                                               <  1:1199597/151‑1 (MQ=255)

   aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTgg                                                                                                                                              >  2:894963/1‑149 (MQ=255)

    cGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAat                                                                                                                                           >  1:940161/1‑151 (MQ=255)

      gCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAatcc                                                                                                                                         >  2:278532/1‑151 (MQ=255)

          tCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAatccatcc                                                                                                                                     >  2:124387/1‑151 (MQ=255)

          tCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAatccatc                                                                                                                                      <  1:261146/150‑1 (MQ=255)

          tCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAatccatc                                                                                                                                      <  3:131855/150‑1 (MQ=255)

           ccTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGGCCAGGTTCCGGGGTCACTGGAAtccatcca                                                                                                                                    <  4:330543/151‑1 (MQ=255)

              gATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGt                                                                                                                                  <  3:182143/150‑1 (MQ=255)

              gATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTc                                                                                                                                 <  3:385920/151‑1 (MQ=255)

              gATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTc                                                                                                                                 <  2:1278238/151‑1 (MQ=255)

                   ccGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTTCTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGAt                                                                                                                            <  2:983095/151‑1 (MQ=255)

                     ggCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCg                                                                                                                          >  2:981408/1‑151 (MQ=255)

                       ccATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGa                                                                                                                         <  2:1195605/150‑1 (MQ=255)

                        cATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGAcc                                                                                                                       <  1:711648/151‑1 (MQ=255)

                           gAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCata                                                                                                                    <  4:364643/151‑1 (MQ=255)

                           gAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCat                                                                                                                     >  2:399977/1‑150 (MQ=255)

                                  aGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCAc                                                                                                             >  3:40998/1‑151 (MQ=255)

                                     cGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTc                                                                                                          <  3:348286/151‑1 (MQ=255)

                                     cGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTc                                                                                                          <  1:316248/151‑1 (MQ=255)

                                      ggTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCg                                                                                                         >  2:1159198/1‑151 (MQ=255)

                                      ggTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCg                                                                                                         <  4:74855/151‑1 (MQ=255)

                                       gTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCgg                                                                                                        <  2:980875/151‑1 (MQ=255)

                                       gTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCgg                                                                                                        <  2:1270052/151‑1 (MQ=255)

                                          aTTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTg                                                                                                     <  2:140131/151‑1 (MQ=255)

                                            tGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCTCTGGAACCCATCCAGTCCGGATCGACCATATTACCACCTCGGCTg                                                                                                     <  2:760107/149‑1 (MQ=255)

                                              ggTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCga                                                                                                 <  2:266635/151‑1 (MQ=255)

                                              ggTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCga                                                                                                 <  1:973509/151‑1 (MQ=255)

                                                   gTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGtt                                                                                            <  2:297/151‑1 (MQ=255)

                                                       aggagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTcc                                                                                        >  2:915148/1‑151 (MQ=255)

                                                         gagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTcc                                                                                        >  1:729379/1‑149 (MQ=255)

                                                          agccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagc                                                                                     >  1:539370/1‑151 (MQ=255)

                                                            ccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcca                                                                                   <  1:26635/151‑1 (MQ=255)

                                                                 cTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCagcagc                                                                              >  2:1172251/1‑151 (MQ=255)

                                                                  tttGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcc                                                                             >  1:10907/1‑151 (MQ=255)

                                                                  tttGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcc                                                                             >  1:607232/1‑151 (MQ=255)

                                                                    tGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGccc                                                                            <  1:453592/150‑1 (MQ=255)

                                                                      atatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCg                                                                          >  1:99026/1‑150 (MQ=255)

                                                                        atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgata                                                                       <  2:397865/151‑1 (MQ=255)

                                                                        atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                        >  2:412846/1‑150 (MQ=255)

                                                                        atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                        >  3:152519/1‑150 (MQ=255)

                                                                          tCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                        >  2:223842/1‑148 (MQ=255)

                                                                            aGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagat                                                                    >  2:1165526/1‑150 (MQ=255)

                                                                              aTCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATc                                                                 >  1:572847/1‑151 (MQ=255)

                                                                                cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                <  2:1087048/150‑1 (MQ=255)

                                                                                cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTg                                                               <  2:515829/151‑1 (MQ=255)

                                                                                 aaaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgg                                                              <  4:396811/151‑1 (MQ=255)

                                                                                 aaaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgg                                                              <  1:230037/151‑1 (MQ=255)

                                                                                 aaaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgg                                                              <  1:830433/151‑1 (MQ=255)

                                                                                  aaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                             <  1:15157/151‑1 (MQ=255)

                                                                                  aaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                             >  1:1150798/1‑151 (MQ=255)

                                                                                  aaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                             >  2:1230597/1‑151 (MQ=255)

                                                                                            cATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATc                                                   <  3:453797/151‑1 (MQ=255)

                                                                                               aaTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCa                                                <  1:651957/151‑1 (MQ=255)

                                                                                               aaTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCa                                                <  4:251319/151‑1 (MQ=255)

                                                                                               aaTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAAGccc                                                 <  2:831371/150‑1 (MQ=255)

                                                                                                aTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCa                                                >  2:1196854/1‑150 (MQ=255)

                                                                                                   gAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaatta                                            <  1:1056647/151‑1 (MQ=255)

                                                                                                   gAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaatta                                            <  2:384875/151‑1 (MQ=255)

                                                                                                   gAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaatta                                            <  2:10907/151‑1 (MQ=255)

                                                                                                    aGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaa                                           <  4:455687/151‑1 (MQ=255)

                                                                                                    aGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaa                                           <  3:195432/151‑1 (MQ=255)

                                                                                                    aGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaatta                                            <  2:798748/150‑1 (MQ=255)

                                                                                                     gTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCATTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaat                                          <  1:381485/151‑1 (MQ=255)

                                                                                                     gTTGTGTATCAGGGTGTCTACGCTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGACTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaactaat                                          <  2:1249583/151‑1 (MQ=255)

                                                                                                      ttGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaatt                                         >  3:18154/1‑151 (MQ=255)

                                                                                                       tgtgTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTg                                        >  2:1136371/1‑151 (MQ=255)

                                                                                                       tgtgTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTg                                        >  2:947725/1‑151 (MQ=255)

                                                                                                        gtgtATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTg                                        <  2:215561/150‑1 (MQ=255)

                                                                                                        gtgtATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGc                                       >  2:376994/1‑151 (MQ=255)

                                                                                                             tCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCtt                                  >  2:858715/1‑151 (MQ=255)

                                                                                                               aGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTta                                <  4:147899/151‑1 (MQ=255)

                                                                                                               aGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTta                                <  3:1867/151‑1 (MQ=255)

                                                                                                                gggTGTCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaa                               <  2:148897/151‑1 (MQ=255)

                                                                                                                  gtgtCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaat                              >  2:657741/1‑150 (MQ=255)

                                                                                                                    gtCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTaataat                           >  3:206044/1‑151 (MQ=255)

                                                                                                                    gtCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTaataat                           >  2:969477/1‑151 (MQ=255)

                                                                                                                     tCTACTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAAtt                          >  2:745345/1‑151 (MQ=255)

                                                                                                                        aCTTTCATCCGTCCAGGTTCCGGGTTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTg                         >  2:515081/1‑149 (MQ=255)

                                                                                                                        aCTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTg                         <  2:869347/149‑1 (MQ=255)

                                                                                                                        aCTTTCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCa                       <  2:1161999/151‑1 (MQ=255)

                                                                                                                          tttCATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGc                     >  2:241746/1‑151 (MQ=255)

                                                                                                                             taTCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGtt                  >  2:1272317/2‑151 (MQ=255)

                                                                                                                             cATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGtt                  <  2:1219721/151‑1 (MQ=255)

                                                                                                                             cATCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGt                   <  1:1253749/150‑1 (MQ=255)

                                                                                                                              aTCCGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTa                 <  4:204041/151‑1 (MQ=255)

                                                                                                                                 cGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTg              <  2:939977/151‑1 (MQ=255)

                                                                                                                                 cGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTg              <  2:366417/151‑1 (MQ=255)

                                                                                                                                 cGTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGAGCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTAtt               <  2:1272032/150‑1 (MQ=255)

                                                                                                                                  gTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTg              >  1:907137/1‑150 (MQ=255)

                                                                                                                                  gTCCAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa             <  2:700952/151‑1 (MQ=255)

                                                                                                                                    ccAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAt           >  2:570272/1‑151 (MQ=255)

                                                                                                                                     cAGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAtt          <  1:516993/151‑1 (MQ=255)

                                                                                                                                      aGGTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAttt         <  2:216142/151‑1 (MQ=255)

                                                                                                                                       ggTTCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTc        <  3:288695/151‑1 (MQ=255)

                                                                                                                                          tCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTcc       >  1:483021/1‑149 (MQ=255)

                                                                                                                                          tCCGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTcc       >  1:320444/1‑149 (MQ=255)

                                                                                                                                            cGGGGTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTc   <  4:407634/151‑1 (MQ=255)

                                                                                                                                             ggggTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTca  <  2:174125/151‑2 (MQ=255)

                                                                                                                                             ggggTCACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCg  >  1:104814/1‑151 (MQ=255)

ATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCG  >  NC_000913/585157‑585448

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