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	96.3%	460.9 / 9.1	136	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/3); new base C (80/51); total (82/54)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAt                                                                                                  <  1:513938/100‑1 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑cc                                                                                               >  1:3297213/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACCTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGACGACTATGTTG‑ACAT‑cc‑ccc                                                                                           >  2:3242711/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑cc‑ccc                                                                                           <  2:1459259/100‑1 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑cc‑ccccctgt                                                                                      >  1:3142264/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑cc‑ccccctgt                                                                                      >  1:2207832/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑cc‑ccccctgt                                                                                      >  1:1768995/1‑96 (MQ=255)

              tcagtgTCACTTCACATAGCTTCAAACTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTg                                                                                   <  2:2898088/95‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTg                                                                                   <  1:478204/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTgg                                                                                  <  2:2113154/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTgg                                                                                  >  2:3300830/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGa                                                                                 >  1:1663894/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGa                                                                                 >  1:723063/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGaa                                                                                >  2:1491517/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑CC‑CCCCTGTTCTGGaaa                                                                               >  2:2737013/1‑100 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGaaa                                                                               >  2:802254/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGGGAAGACTATGTTG‑AAATCCC‑CCCCTGTTTTGGTAAATg                                                                            >  1:2310606/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATg                                                                            <  1:3332796/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGa                                                                           >  2:1425426/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGa                                                                           >  1:3600576/1‑100 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACt                                                                         <  2:2505673/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTa                                                                        >  2:689772/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAg                                                                       >  1:3086147/1‑100 (MQ=255)

                           acaTAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGc                                                                      <  1:829900/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     >  2:834211/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:1822884/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:1746243/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 <  1:3471994/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:1805717/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:186750/1‑100 (MQ=255)

                                 ctttgaaTTCGTCCCCCATTGTCTTCGTATCCTGCTGCCAGTGCAAAGGAGATGACTATGTTGTAGTT‑CC‑CCCCTGTTCTGGAAAATGCCTAGCCTGttt                                                                <  2:3428364/95‑1 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  1:1818936/100‑1 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:2862038/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:3725761/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTg                                                                   >  2:917545/1‑95 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:2826999/100‑1 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  2:3697600/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1050386/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:153686/100‑1 (MQ=255)

                                         tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:1187506/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTATGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:1646073/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAATGATAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:665367/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:2462556/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTTGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:94739/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:1363587/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  1:2163350/1‑100 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             >  1:3136997/1‑100 (MQ=255)

                                                     ntCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCCGTTTGGTGCAGTATTGCtgttgt                                             >  1:3386989/2‑99 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:1739210/99‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2983841/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  1:43814/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:3486518/1‑99 (MQ=255)

                                                       cTTCGTATCAAGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTATTCTGGAAAATGAATAGCCAGTTTGGTGCAGTATTGCTGttgttg                                           <  2:3798162/99‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  1:2037867/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGACTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  2:850838/1‑100 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:2829653/1‑100 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:2281451/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  2:1155420/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  2:3147768/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTATTGCTGTTGTTGATTCCg                                     >  1:1316402/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:674583/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  2:2099223/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  1:1250101/1‑100 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AATTCCC‑CCCCTGTTCTTGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:3578594/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGGAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  2:2331370/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:1490577/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 <  1:3697628/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  2:2374680/99‑1 (MQ=255)

                                                                     gcCATTGCAGACGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGCAAATGACTAGCCTGTTTGGTGCAGTATTCCCGTTTTTGATTcctctacag                               >  2:3166008/1‑90 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:2642408/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:277440/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1391258/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  1:854267/1‑100 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          <  1:1567775/100‑1 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  1:602294/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  2:251406/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:1094070/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCCGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1631420/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:2596373/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:841551/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2603697/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  2:1190542/100‑1 (MQ=255)

                                                                               aGGAGAAGACTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  2:1570646/100‑1 (MQ=255)

                                                                                  agaagaCTATGTTG‑AAATCCC‑CCCCTCTTCTGGAAAATGACTAACCTATTGGGTCCGGTGTTTGTGGTGTTGATTGCGATAATGGTGATTCGGCAGGTGa               >  2:2633378/1‑100 (MQ=255)

                                                                                    aagaCTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt             <  2:2800664/100‑1 (MQ=255)

                                                                                     agaCTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:1096783/1‑100 (MQ=255)

                                                                                      gaCTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:1770093/99‑1 (MQ=255)

                                                                                      gaCTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:3711986/1‑100 (MQ=255)

                                                                                        cTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  2:2067304/1‑100 (MQ=255)

                                                                                        cTATGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:2016133/1‑100 (MQ=255)

                                                                                          aTGTTG‑AAAT‑CCGTCTCTGTTTTGGAAAATTATTAGTCTGTTTGGTGCAGTATTGTTTTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  2:2211344/100‑1 (MQ=255)

                                                                                           tGTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:229608/1‑100 (MQ=255)

                                                                                            gTTG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:2958458/1‑100 (MQ=255)

                                                                                            gTTG‑AAATCCC‑CCCCTGTTCTGGAAAAGGACTACCCTGTTTGGTGCAGAATTGCTGTTGTTGCTTCCGATCGTGCTGATTCGGCAGGTGATTGTCGAAc      >  2:221358/1‑99 (MQ=255)

                                                                                            gTTG‑AAATCCC‑CCCCTGTTCAGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:3591724/1‑100 (MQ=255)

                                                                                             ttG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGCTAATGCTGCTTCGGCAGGTGATTTTCTTACGt    >  2:527629/1‑100 (MQ=255)

                                                                                             ttG‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:538106/1‑100 (MQ=255)

                                                                                             tccgatct‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    <  1:247056/93‑1 (MQ=255)

                                                                                               g‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:140241/1‑100 (MQ=255)

                                                                                               g‑AAATCCC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:2779938/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTG‑AAAT‑CC‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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