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	98.9%	329.9 / ‑3.7	91	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (54/36); total (54/37)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑t                                                                                                 <  1:1903286/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑c                                                                                               >  1:1696162/1‑100 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑c                                                                                               >  1:1577742/1‑100 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cccc                                                                                            <  2:921720/100‑1 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                           >  2:1854169/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTGCCCCCATAGTCTTCGTAGCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccc                                                                                           <  1:1750972/100‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑cccccc                                                                                          >  1:1004679/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCGAAGGAGAAGACTATGTTGAAA‑T‑ccccccc                                                                                         <  2:2333557/100‑2 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccccc                                                                                         >  2:490918/1‑97 (MQ=255)

         cGAATTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccccctg                                                                                       <  2:2166637/100‑4 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATACTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccccccctgtt                                                                                     >  2:1121108/1‑95 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTc                                                                                     >  2:396783/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTg                                                                                   <  1:14810/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTg                                                                                   >  1:1924761/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTACCCGCATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTg                                                                                   <  2:443975/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaa                                                                                >  2:1268648/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaa                                                                                >  1:1406101/1‑100 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaaa                                                                               <  2:1347138/100‑1 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaaaa                                                                              >  1:687637/1‑100 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGaaaa                                                                              <  1:248276/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCGCATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAAt                                                                             <  2:1821883/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATg                                                                            >  2:538605/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGa                                                                           <  1:1212924/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGAc                                                                          <  1:1768452/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGAc                                                                          <  1:1861465/100‑1 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAg                                                                       >  1:181205/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAg                                                                       >  1:510902/1‑100 (MQ=255)

                           acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGc                                                                      <  2:1471126/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:801405/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  2:801402/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:119790/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:221881/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:2242829/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:2308254/1‑100 (MQ=255)

                                  ttCAAATTCTTCGCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:1857619/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAATGACAAGACTATTTAGAAA‑TCCCCCCCCGTTCTGGAAAATGACTAGcc                                                                     >  2:1605661/1‑94 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:1960010/1‑100 (MQ=255)

                                   tCAAATTCTTTCCACATAATCTTTGTATCCTGCTGCCATTGCAAAGGGGCAGACTATGTTTAAA‑TCCCCCCCTTTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:347964/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:396786/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:821274/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:2350809/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1530835/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  2:1859634/1‑100 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  1:762645/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                      >  1:1770105/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:2313162/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:1463913/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:525457/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:335862/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:1099344/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             >  2:255925/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:2053383/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  1:679300/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  2:1389569/1‑100 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       >  2:1944637/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:1458370/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:1579189/100‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  2:221880/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  1:1000395/1‑100 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:1103772/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:254809/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:1787159/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  2:523326/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:2160331/1‑100 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:302850/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAGA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:2094326/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:920239/100‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGCCTAGCGTGTTTGGTGCAGTATTGCAGTTGTTGATTCCGATAATGCTGATTc                        >  2:1106197/1‑98 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGATCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:48738/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:324238/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:2259478/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:197643/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:410537/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  2:1701696/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  1:156825/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:1854182/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  1:1671482/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:1107454/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:1110580/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:1105545/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:647153/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:57138/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:57006/1‑100 (MQ=255)

                                                                                         tATGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        <  1:175678/100‑1 (MQ=255)

                                                                                          aTGTTGGAA‑TCCCCCCCGGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  2:468158/100‑1 (MQ=255)

                                                                                          aTGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  2:2269558/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:856551/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:603593/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:79767/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:636488/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCGTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  2:575306/100‑1 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:39923/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:1956912/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:2333135/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:252715/1‑100 (MQ=255)

                                                                                            gTTGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1683290/1‑100 (MQ=255)

                                                                                             ttGAAA‑TCCCCCCCTGTTCTGGAAAATGACTATCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2233983/1‑100 (MQ=255)

                                                                                             ttGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2166900/1‑100 (MQ=255)

                                                                                              tGAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  1:1522988/1‑100 (MQ=255)

                                                                                               gAAA‑TCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:265219/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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