BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,585,157:1	+C	19.9%	coding (1330/2280 nt)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,156	1	.	C	19.9%	406.0 / 124.5	141	coding (1331/2280 nt)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (46/67); new base C (20/8); total (66/75)
Fisher's exact test for biased strand distribution p-value = 5.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 152 total aligned reads displayed.

CGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGT  >  NC_000913/1585058‑1585255

cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑c                                                                                                      <  2:2795480/100‑1 (MQ=255)

   tcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                                   <  1:3869331/100‑1 (MQ=255)

    cgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagca                                                                                                  >  1:3294802/1‑100 (MQ=255)

    cgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagca                                                                                                  <  1:1687238/100‑1 (MQ=255)

      ctTCCCGATCTGGCATTGCCAGCGCAAAATTGTCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATg                                                                                                <  1:1451993/100‑1 (MQ=255)

      ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                                   >  2:3706650/1‑97 (MQ=255)

      ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATg                                                                                                <  1:2717858/100‑1 (MQ=255)

      ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATg                                                                                                <  1:1681575/100‑1 (MQ=255)

       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTTCCCCGTACATATTGCTTCCATGCTGGCAATTg‑cagc                                                                                                   >  2:4173354/1‑96 (MQ=255)

       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGt                                                                                               <  2:906092/100‑1 (MQ=255)

        tCCCCATCCGGCAGTGCCAGCGCAAAATTGCCCCCCACGGAGATCAATTCTCGAGCCCGCCCCGTACATATTGCTTGCTTGCT‑GCAAGTGCCAGCATgc                                                                                               <  2:2565539/99‑2 (MQ=255)

         tttgcgcTGACAATGCCAGCGCAAAATAGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCAGGACATATTACTTGCATGCCGGCAATTG‑CAGCCTTTcc                                                                                             >  2:235990/7‑100 (MQ=255)

          ccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                                            >  2:4157488/1‑100 (MQ=255)

          ccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaag                                                                                            <  2:3604187/100‑5 (MQ=255)

           cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaagc                                                                                           <  1:1562927/100‑6 (MQ=38)

           cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCAAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaan                                                                                            <  2:3974558/99‑5 (MQ=255)

            gATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAgg                                                                                          <  1:1566716/100‑1 (MQ=255)

               cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGc                                                                                       >  2:2017938/1‑100 (MQ=255)

                tGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCa                                                                                      <  2:3261852/100‑1 (MQ=255)

                 ggCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAt                                                                                     <  1:3581439/100‑1 (MQ=255)

                  gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATg                                                                                    <  1:3179182/100‑1 (MQ=255)

                      tGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGt                                                                                <  2:1170107/100‑1 (MQ=255)

                       gCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTg                                                                               >  2:2643974/1‑100 (MQ=255)

                       gCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTg                                                                               >  2:1541918/1‑100 (MQ=255)

                        ccAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgg                                                                              >  1:499067/1‑100 (MQ=255)

                         cAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                             <  1:982133/100‑1 (MQ=255)

                         cAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                             <  1:3052536/100‑1 (MQ=255)

                          aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAtggggtgggg                                                                            <  2:190470/100‑1 (MQ=255)

                          aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                            <  2:3755707/100‑1 (MQ=255)

                          aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                            <  2:1082448/100‑1 (MQ=255)

                           gcgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGt                                                                           >  1:2892091/1‑100 (MQ=255)

                            cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                                          <  2:1354850/100‑1 (MQ=255)

                            cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                                          <  1:1019055/100‑1 (MQ=255)

                             gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGa                                                                         <  1:3397420/100‑1 (MQ=255)

                             gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGa                                                                         >  1:3169206/1‑100 (MQ=255)

                               aaaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                       <  1:1828904/100‑1 (MQ=255)

                               aaaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTGGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                       <  2:753977/100‑1 (MQ=255)

                                aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGagcg                                                                      >  1:426737/1‑97 (MQ=255)

                                aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                      >  1:3434822/1‑100 (MQ=255)

                                 aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGcc                                                                     <  1:1036251/100‑1 (MQ=255)

                                  aTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCa                                                                    >  2:2449534/1‑100 (MQ=255)

                                    tGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                  <  2:209718/100‑1 (MQ=255)

                                     gCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAg                                                                 <  1:3445133/100‑1 (MQ=255)

                                       cccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                                               >  2:3295622/1‑100 (MQ=255)

                                       cccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGAATGAGGTGGAGTGAAGCCATAgc                                                                >  1:326094/1‑99 (MQ=255)

                                        ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                                              >  2:1060245/1‑100 (MQ=255)

                                         cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCt                                                             >  2:2767987/1‑100 (MQ=255)

                                          ccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTc                                                            >  1:535361/1‑100 (MQ=255)

                                           cATGCAGATCAATGCTCGCGCCTGCCCCGTACATATTGCTTGCATGCTTGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCa                                                           <  1:4032009/100‑1 (MQ=255)

                                            aTGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGCCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                          <  2:2005750/100‑1 (MQ=255)

                                            aTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                          <  2:499066/100‑1 (MQ=255)

                                              gCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAccc                                                        <  2:3322113/100‑1 (MQ=255)

                                                   tCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACg                                                   >  2:3204562/1‑100 (MQ=255)

                                                    cAATGCTCGAGCCTGCCCCGAACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGa                                                  >  2:2556634/1‑100 (MQ=255)

                                                     aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAg                                                 >  1:710964/1‑100 (MQ=255)

                                                        gCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                              <  2:2734044/100‑1 (MQ=255)

                                                        gCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                              <  1:754074/100‑1 (MQ=255)

                                                         cTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAg                                             >  1:3656257/1‑100 (MQ=255)

                                                             aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc                                         <  2:372888/100‑1 (MQ=255)

                                                             aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc                                         <  1:2260534/100‑1 (MQ=255)

                                                             aGCCTGCCCCGTACATAGTGCTGTCATGCGGGCAATGG‑CAGCACGTCCCGGTGCCTGAGGTGGCGTGAAGTCATAGCGCCCACCCAGACGAGCCAGAAAc                                         <  2:3809109/100‑1 (MQ=255)

                                                                cTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc‑‑tct                                    <  1:4067034/100‑1 (MQ=255)

                                                                cTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc‑‑tct                                    <  1:3956017/100‑1 (MQ=255)

                                                                 tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTg                                   >  2:1580592/1‑100 (MQ=255)

                                                                 tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTg                                   <  2:179633/100‑1 (MQ=255)

                                                                  gCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGc                                  <  1:1523020/100‑1 (MQ=255)

                                                                  gCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGc                                  <  2:292041/100‑1 (MQ=255)

                                                                  gCCCCGTACATATTACTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGc                                  <  1:1452656/100‑1 (MQ=255)

                                                                    cccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATAAGGTGGGGTGAAGCCATAGCGCTCACCCCGACGAGCCAGAAAC‑‑TCTGCAg                                >  1:1895902/1‑100 (MQ=255)

                                                                       gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGGCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACg                             <  1:2139319/100‑1 (MQ=255)

                                                                              ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTtctc                      <  2:3540501/100‑1 (MQ=255)

                                                                               tgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCg                     >  2:3154358/1‑100 (MQ=255)

                                                                               cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCg                     >  1:1526733/7‑100 (MQ=38)

                                                                               cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCg                     >  2:2784473/7‑100 (MQ=38)

                                                                               cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCg                     >  1:3413958/7‑100 (MQ=38)

                                                                                gcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACCGTCTGCAGagtttctggc                      <  1:357456/100‑10 (MQ=255)

                                                                                gcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCgg                    <  1:3271819/100‑1 (MQ=255)

                                                                                gcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCgg                    <  1:2221528/100‑1 (MQ=255)

                                                                                ctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCgg                    >  1:1857986/6‑100 (MQ=38)

                                                                                ctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCgg                    >  1:1991210/6‑100 (MQ=38)

                                                                                 tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   >  1:529717/5‑100 (MQ=255)

                                                                                 tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   >  2:4118002/5‑100 (MQ=255)

                                                                                 tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   >  1:373805/5‑100 (MQ=255)

                                                                                 tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   >  1:3315476/5‑100 (MQ=255)

                                                                                 tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   >  1:2087525/5‑100 (MQ=255)

                                                                                 cttgcGTGCTGGCAACTG‑CAGCGTGTCCAGGTGCATCAGGTGGGGTGACGCCATAGCGCTTGCCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGt                   <  2:473167/100‑1 (MQ=255)

                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTg                  <  1:2405164/100‑1 (MQ=255)

                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTg                  <  1:2465382/100‑1 (MQ=255)

                                                                                   gacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGa                 >  2:2456730/3‑100 (MQ=255)

                                                                                    acATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGAt                >  2:1247933/2‑100 (MQ=255)

                                                                                     cATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATa               <  1:2419064/100‑1 (MQ=255)

                                                                                     cATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTAACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATa               >  2:693535/1‑100 (MQ=255)

                                                                                     cATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATa               >  2:2191382/1‑100 (MQ=255)

                                                                                     cATGATGTCAATTT‑CAGCAGGTCCAAGTGCATGAGGGTGGGTGAAGACATAGTGCTCTCCCAGACGAGCCAGAACC‑‑CCTGCAGACGGTTTCTCGGTg                  <  2:681693/97‑1 (MQ=255)

                                                                                             cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGt       >  2:1558296/1‑100 (MQ=255)

                                                                                             cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGt       <  1:2977769/100‑1 (MQ=255)

                                                                                             cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGt       >  2:1279427/1‑100 (MQ=255)

                                                                                              aaTTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACg         >  2:154095/1‑97 (MQ=255)

                                                                                                ttG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACGGACGGTTCg    >  2:2354011/1‑100 (MQ=255)

                                                                                                  g‑cagcaTGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGt  >  2:773131/1‑100 (MQ=255)

CGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAC‑‑TCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGT  >  NC_000913/1585058‑1585255

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: `‑`