BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)_6→7	100%	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	91.3%	430.7 / 36.4	138	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (5/7); new base C (67/58); total (73/65)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

   cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c                                                                                              <  1:533079/100‑1 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑c                                                                                            <  2:2819557/100‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑cc                                                                                           >  1:2120170/1‑100 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑ccc                                                                                          <  2:1135673/100‑2 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCGTCGTATCCTGCTGCCATTGCAAAGGAGAATACTATGTTGAAAT‑c‑cc‑c‑ccc                                                                                          <  2:2381282/100‑2 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑ccctgt                                                                                       >  2:2906159/1‑96 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑ccctgt                                                                                       >  2:2219481/1‑96 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑ccctgtt                                                                                      >  1:1461907/1‑95 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc‑c‑ccctgtt                                                                                      >  1:2460689/1‑95 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTc                                                                                     >  1:2089796/1‑100 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTTTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTg                                                                                   <  1:3094670/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCt                                                                                    >  2:143624/1‑99 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGCAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTg                                                                                   <  1:3730682/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGGAATCC‑CC‑C‑CC‑TGTTCTgg                                                                                  >  2:2678352/1‑100 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTgg                                                                                  >  2:435116/1‑100 (MQ=255)

                ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGaa                                                                                <  2:618637/100‑1 (MQ=255)

                   tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAAt                                                                             >  2:2595439/1‑100 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATg                                                                            <  2:3869774/100‑1 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATg                                                                            <  2:1676690/100‑1 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAAt                                                                             <  2:3930670/98‑1 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGa                                                                           >  2:283278/1‑100 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGAATATGTTGAAAT‑CACC‑C‑CC‑TGTTCTGGAAAATGa                                                                           >  1:3439334/1‑100 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTACAATTGCAAAAAAGAAGACTATGTTGAAAT‑C‑CC‑CTCC‑TGTTCTGGAAAATGa                                                                           >  1:1015985/1‑100 (MQ=255)

                       ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACt                                                                         >  2:277172/1‑100 (MQ=255)

                       ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCAGGAAAATGACt                                                                         <  2:1781145/100‑1 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGcc                                                                     >  1:2879099/1‑100 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCAGGAAAATGACTAGc                                                                      >  1:3867558/1‑99 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCTC‑CC‑TGTTCTGGAAAATGACTAGcc                                                                     <  2:1247564/100‑1 (MQ=255)

                             tctcTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTg                                                                   <  2:3270117/97‑1 (MQ=255)

                             tAGCTTCAAATTCTTCCCACATAGTCTTTGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTg                                                                   <  1:2170684/100‑1 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:1953705/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTCTGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCATGttt                                                                >  2:2746230/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:3741006/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGACATTGCCAAGGAGCAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:774913/100‑1 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:3443899/1‑99 (MQ=255)

                                   cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  1:2890390/100‑1 (MQ=255)

                                   cAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:813435/1‑100 (MQ=255)

                                    aaaTTCTTCCCATATAGTGTTCGTATCTTGTTGTTTTTTCAAAGGAGAAGAGTATGGTGAAAT‑C‑CC‑C‑CCATGGTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:1108902/100‑1 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:311393/100‑1 (MQ=255)

                                      aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  1:2621235/1‑100 (MQ=255)

                                       ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:2029291/100‑1 (MQ=255)

                                        tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        >  2:3648793/1‑100 (MQ=255)

                                          ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATATTGAAAT‑C‑CCAC‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:1057256/1‑100 (MQ=255)

                                            cctacaTAGTCTTCGTATACTGCTGCCATTGCAAATGAGAAGAGTATGTTGAAATAC‑CC‑G‑CC‑TGTTCTGGAAAATGACTAGCCTATTTGTTGCaaaatt                                                     >  2:1349704/4‑94 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:283284/100‑1 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  <  1:277177/100‑1 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:1980164/1‑100 (MQ=255)

                                                caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:1206295/100‑1 (MQ=255)

                                                caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  1:3883861/100‑1 (MQ=255)

                                                caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  1:2594074/99‑1 (MQ=255)

                                                 aTAGTCTTTGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGTAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               <  2:3338402/100‑1 (MQ=255)

                                                 aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  1:1937096/99‑1 (MQ=255)

                                                     ttttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:3761569/98‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATTTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1535568/99‑1 (MQ=255)

                                                     tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:1544547/100‑1 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:746211/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:3911428/1‑100 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:666795/1‑100 (MQ=255)

                                                         cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  2:1810789/100‑1 (MQ=255)

                                                         cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       >  2:251510/1‑100 (MQ=255)

                                                            aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:1645133/1‑100 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:558667/1‑100 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  1:3348440/1‑100 (MQ=255)

                                                                tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  1:3261963/1‑100 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:2846756/100‑1 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCTC‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:2289009/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:666787/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAAAAAATTATGTTTAAAT‑C‑CCTC‑CC‑TGTTCTTGAAAATGACTAGCCTGTTTGGTGCATTATTGCTGTTGTTGATTCCGATAATGc                              <  2:3161212/99‑1 (MQ=255)

                                                                   ctccATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:3381921/98‑1 (MQ=255)

                                                                     ccATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGAAGTTGTTGATTCCGATAATGCTGa                           >  1:3216789/1‑100 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  1:3360636/1‑100 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  2:2372655/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATgagatcggaa                      >  2:3464279/1‑91 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:505286/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:3702323/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGCTTGGTGCAGTATTGCTGTAGTTGATTACGATAATGCTGATTCg                       >  1:837418/1‑99 (MQ=255)

                                                                           cAAATGAGAAGACTATTTTGAAAT‑C‑CC‑CTCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      <  1:1666302/99‑1 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  2:2935308/1‑98 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:340036/1‑100 (MQ=255)

                                                                           cAAAGAAGAAGACTATGTTGAAAT‑C‑CCGC‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:3540620/98‑1 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAttn                        <  2:1225967/96‑2 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:1169152/100‑1 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTACAGTATTGCTGTTGATGATTCCGATAATGCTGATTCGGCa                    >  2:3305419/1‑100 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  1:1338472/1‑100 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGGTTGGTGCAGGATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:3558423/1‑100 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  1:3820355/1‑100 (MQ=255)

                                                                                  gaagaCTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAAAGACTAACCTGTTTGGTGCAGTATTGCTGTTGTTGACTCCGATAATGCTGATTCGGCAGGTGAt              <  1:2746254/100‑1 (MQ=255)

                                                                                       cTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGTCTGTTTGGTGCAGTATTGCTGGTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:3443938/99‑1 (MQ=255)

                                                                                       cTATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  2:2684597/100‑1 (MQ=255)

                                                                                        tATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        >  1:3520586/1‑100 (MQ=255)

                                                                                        tATGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        >  1:2453157/1‑100 (MQ=255)

                                                                                         aTGTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  1:2994629/1‑100 (MQ=255)

                                                                                           gTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2187273/1‑100 (MQ=255)

                                                                                           gTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1696994/1‑100 (MQ=255)

                                                                                           gTTGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  2:2994602/100‑1 (MQ=255)

                                                                                             tGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  2:1607725/1‑100 (MQ=255)

                                                                                             tGAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   <  1:883792/100‑1 (MQ=255)

                                                                                              gAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTgn  <  2:3787074/100‑2 (MQ=255)

                                                                                              gAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTga  >  2:883781/1‑99 (MQ=255)

                                                                                              gAAATCC‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:3629709/100‑1 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CC‑C‑CC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

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