BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,379,796	C→T	43.9%	R283H (CGC→CAC)	menF ←	isochorismate synthase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,379,796	0	C	T	43.9%	79.3 / 196.4	167	R283H (CGC→CAC)	menF	isochorismate synthase 2
Reads supporting (aligned to +/- strand): ref base C (60/33); new base T (36/37); total (97/70)
Fisher's exact test for biased strand distribution p-value = 5.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 173 total aligned reads displayed.

ACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACG  >  NC_000913/2379697‑2379893

aCTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcgc                                                                                                   >  2:3189786/1‑100 (MQ=255)

 cTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTcgcg                                                                                                  <  1:185445/100‑1 (MQ=255)

    tACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTg                                                                                               <  2:3462790/100‑1 (MQ=255)

      cgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCGCCAGCATGTTCTCGTGCTGGt                                                                                             >  2:3675430/1‑100 (MQ=255)

      cgcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCAAGTTCTCGCGCTgg                                                                                              >  1:2422589/1‑99 (MQ=255)

       gcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGtt                                                                                            >  2:131100/1‑100 (MQ=255)

       gcagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGtt                                                                                            >  1:1632145/1‑100 (MQ=255)

        cagacgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGttt                                                                                           >  2:2564824/1‑100 (MQ=255)

           acgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTa                                                                                        <  2:4220672/100‑1 (MQ=255)

           acgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTa                                                                                        <  1:1437401/100‑1 (MQ=255)

            cgcagTACCTGCGGCGGTGAAACATCCAGCGTCTTGGTATCTGCCTGTTATCTTTGACAGATTTCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTa                                                                                        <  1:3652302/99‑1 (MQ=255)

            cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTat                                                                                       <  1:2147456/100‑1 (MQ=255)

            cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTa                                                                                        >  1:2189233/1‑99 (MQ=255)

            cgcagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTat                                                                                       >  2:3591645/1‑100 (MQ=255)

              cagTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTatca                                                                                     >  2:682487/1‑100 (MQ=255)

                 tACCTGAGGCGGTAAAAAATCCAGCGTCTGGGTACCTGCCTGGAATCGTTGACAGATCTCGTCGACCACCAGCATGTTGTCGCGCTGGTTTTTATCATcc                                                                                  >  1:4058411/1‑100 (MQ=255)

                  aCCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCg                                                                                 >  2:1714416/1‑100 (MQ=255)

                  aCCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCg                                                                                 >  2:1114546/1‑100 (MQ=255)

                    cggcggcgGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGcc                                                                               <  1:971137/98‑1 (MQ=255)

                    cTGCGGCGGTAAGACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGcc                                                                               >  2:2804838/1‑100 (MQ=255)

                    cTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGcc                                                                               <  2:3834142/100‑1 (MQ=255)

                     tGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCa                                                                              >  2:700258/1‑100 (MQ=255)

                      gcggcgGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCAt                                                                             >  1:2673612/1‑100 (MQ=255)

                         gcggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAg                                                                          <  1:3189774/100‑1 (MQ=255)

                          tggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGc                                                                         <  2:2848111/99‑1 (MQ=255)

                          cggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGc                                                                         <  2:388226/100‑1 (MQ=255)

                           ggTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCTCTGGTTTTTATCATCCGCCATCAGcc                                                                        <  2:2597981/100‑1 (MQ=255)

                            gTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCa                                                                       >  2:1364548/1‑100 (MQ=255)

                             tAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCATCCAc                                                                      >  1:2760234/1‑100 (MQ=255)

                              aaaaCATCCAGCGTCTGGGTATCGTCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCAct                                                                     <  1:1384717/100‑1 (MQ=255)

                              aaaaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTTACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCAct                                                                     <  1:3121695/100‑1 (MQ=255)

                              aaaaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCAct                                                                     <  1:2859967/100‑1 (MQ=255)

                              aaaaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCAct                                                                     <  2:2811946/100‑1 (MQ=255)

                              aaaaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCAct                                                                     >  1:1026675/1‑100 (MQ=255)

                                aaCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCActct                                                                   <  2:2610759/100‑1 (MQ=255)

                                 aCATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGCTTTTATCATCCGCCATCAGCCACtctc                                                                  >  2:337427/1‑100 (MQ=255)

                                    tCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTa                                                               <  2:2685299/100‑1 (MQ=255)

                                    tCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTa                                                               <  2:4161012/100‑1 (MQ=255)

                                    tCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTa                                                               >  1:2798106/1‑100 (MQ=255)

                                    tCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTa                                                               <  1:3369/100‑1 (MQ=255)

                                     ccAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTaa                                                              <  1:2401909/100‑1 (MQ=255)

                                      cAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAc                                                             >  2:2405772/1‑100 (MQ=255)

                                       aGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAAct                                                            <  1:840582/100‑1 (MQ=255)

                                           tCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAActgctg                                                        <  2:3355267/100‑1 (MQ=255)

                                           tCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAActgctg                                                        >  2:2338478/1‑100 (MQ=255)

                                            cTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCGCCACCAGCATGTTCGCGGGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACtgctgc                                                       >  1:4124841/1‑100 (MQ=255)

                                            cTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACtgctgc                                                       >  2:3920483/1‑100 (MQ=255)

                                             tGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTgcg                                                      >  2:431401/1‑100 (MQ=255)

                                              gggTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTgcgc                                                     >  2:765684/1‑100 (MQ=255)

                                               ggtcctctccTCTAATCGTTGACAGTTATCTCCCCCCACCATCATCTTCTCGCGCTGCTTTCTCTCATCCGCCATCAGCCACTCTCCTAActg                                                           <  1:3985048/85‑1 (MQ=255)

                                               ggTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGcc                                                    <  1:3261516/100‑1 (MQ=255)

                                                 tATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTg                                                  >  1:2153467/1‑100 (MQ=255)

                                                  aTCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTTGTTTTTATCATCCGCCATCAGACACGCTCCTAACTGCTGCGCCTGc                                                 <  1:1912188/100‑1 (MQ=255)

                                                  aTCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGc                                                 >  2:2992347/1‑100 (MQ=255)

                                                   tCGGCCTGTACTCGTTGACAGATATCTTCCACCGCCAGCATATTCTCGCGCTGGTTTTTCTCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCt                                                <  1:3143008/100‑1 (MQ=255)

                                                   tCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCt                                                >  1:538420/1‑100 (MQ=255)

                                                     ggCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTa                                              <  1:1102081/100‑1 (MQ=255)

                                                       ccTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatc                                            <  2:516159/100‑1 (MQ=255)

                                                         tGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatcat                                          <  1:3120598/100‑1 (MQ=255)

                                                         tGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTatcat                                          <  2:2353900/100‑1 (MQ=255)

                                                           tAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTAtcatca                                        <  2:1632149/100‑1 (MQ=255)

                                                            aaTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAg                                       <  1:3650349/100‑1 (MQ=255)

                                                             aTCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAgg                                      <  2:651757/100‑1 (MQ=255)

                                                              tCGTTGACAGATTTCTGCCACCACGAGCGTGTTGTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGa                                     <  1:3237101/100‑1 (MQ=255)

                                                                 ttGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGatt                                   >  1:697693/1‑99 (MQ=255)

                                                                  tGACAGATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGattat                                 <  2:1997064/100‑1 (MQ=255)

                                                                  tGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGattat                                 >  1:1849816/1‑100 (MQ=255)

                                                                     cAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTg                              >  1:1168993/1‑100 (MQ=255)

                                                                      aGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGc                             <  1:1342125/100‑1 (MQ=255)

                                                                      aGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGc                             <  2:3049346/100‑1 (MQ=255)

                                                                       gATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCt                            >  2:1076120/1‑100 (MQ=255)

                                                                       gATATCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCt                            >  2:4087836/1‑100 (MQ=255)

                                                                        atatCTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTa                           <  2:3268406/100‑1 (MQ=255)

                                                                            cTTCCACCACCAGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGt                       >  1:1041458/1‑100 (MQ=255)

                                                                            cTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGt                       >  2:1458037/1‑100 (MQ=255)

                                                                               ccaccaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTcc                    <  1:3767056/100‑1 (MQ=255)

                                                                                caccacctGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTccc                   >  2:2240242/1‑100 (MQ=255)

                                                                                caccaccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTccc                   >  2:304761/1‑100 (MQ=255)

                                                                                caccaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCCACTGTTccc                   >  2:3791646/1‑100 (MQ=255)

                                                                                 accaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCg                  >  1:1243248/1‑100 (MQ=255)

                                                                                  ccaccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGc                 >  2:1538553/1‑100 (MQ=255)

                                                                                   caccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGcc                >  2:2502349/1‑100 (MQ=255)

                                                                                   caccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGcc                >  2:505010/1‑100 (MQ=255)

                                                                                   caccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGcc                <  2:626058/100‑1 (MQ=255)

                                                                                   caccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGcc                <  2:538422/100‑1 (MQ=255)

                                                                                   caccaGCATGTTCACGTGCTGGTTTTTATCATCCGCCATCACCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCCACTCTTccc                   >  1:1209793/1‑97 (MQ=255)

                                                                                     ccaGCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAg              <  1:2580576/100‑1 (MQ=255)

                                                                                     ccaGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAg              >  2:198942/1‑100 (MQ=255)

                                                                                       aGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAgcg            >  1:4028370/1‑100 (MQ=255)

                                                                                       aGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAgcg            >  2:2272098/1‑100 (MQ=255)

                                                                                        gCATGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAgcgc           >  1:2294970/1‑100 (MQ=255)

                                                                                         cATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCt          <  1:1145849/100‑1 (MQ=255)

                                                                                           tGTTCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTc        >  2:3358014/1‑100 (MQ=255)

                                                                                           tGTTCTCGCGCTTGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTc        <  1:3585870/100‑1 (MQ=255)

                                                                                           tGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTc        <  2:3534822/100‑1 (MQ=255)

                                                                                            gTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCa       >  2:2105790/1‑100 (MQ=255)

                                                                                             ttCTCGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAg      <  2:2354975/100‑1 (MQ=255)

                                                                                             ttCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAg      >  2:888377/1‑100 (MQ=255)

                                                                                               ctcGTGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTa    >  2:4115026/1‑100 (MQ=255)

                                                                                                 cgcgCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACg  >  2:4242660/1‑100 (MQ=255)

ACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAGATATCTTCCACCACCAGCATGTTCTCGCGCTGGTTTTTATCATCCGCCATCAGCCACTCTCCTAACTGCTGCGCCTGCTTATCATCAGGATTATTTGCTACTGTTCCCGCCAGCGCTTCAGTACG  >  NC_000913/2379697‑2379893

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