BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	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.3%	417.7 / 1.0	117	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (59/56); total (61/56)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 131 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                <  2:530909/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                              <  2:1290943/100‑1 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                             >  1:526689/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                             >  2:1859922/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                             <  2:2366426/100‑1 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                             >  2:2872605/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                          >  2:1428902/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                          <  1:988303/100‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGTAAAGGAGAAGACTATGTTGAAAT‑cccccc                                                                                         >  1:2741895/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccccc                                                                                         >  2:134026/1‑100 (MQ=255)

        tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccct                                                                                       <  1:1072247/100‑3 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccctgtt                                                                                    >  2:96340/1‑95 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTc                                                                                    >  1:667299/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTc                                                                                    >  1:1723819/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                  >  2:2794520/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                  <  1:2859843/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTgg                                                                                 >  1:2167710/1‑100 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTgg                                                                                 >  2:839288/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGaa                                                                               <  1:2081881/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCGGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGTaa                                                                               <  2:446017/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGa                                                                          >  1:2903626/1‑100 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACt                                                                        <  1:186862/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTa                                                                       <  1:2131652/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTa                                                                       >  2:2898943/1‑100 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTa                                                                       >  1:859030/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    >  1:1489156/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  2:1924013/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  1:2370694/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    >  1:65561/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                    >  2:2191241/1‑100 (MQ=255)

                              tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTg                                                                  <  2:2870811/100‑1 (MQ=255)

                              tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTg                                                                  <  1:722329/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                >  2:314636/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               >  2:338620/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               >  1:1907522/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                              >  2:2269079/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             <  2:2821721/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             >  1:949162/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             >  2:1782922/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                            <  1:2739476/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                            <  2:949165/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTGTTCGTATCCTGCTGCCATGGCAAAGGAGAAGATTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGGGc                                                          <  1:2621366/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                         <  2:2797025/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                         <  2:552557/100‑1 (MQ=255)

                                           ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCGTGTTTGGTTCAGTAt                                                     >  1:2758395/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   >  1:2809039/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  >  1:2265366/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  <  1:219902/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  2:427111/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  1:1197960/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                <  2:2587101/100‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                               >  2:1805694/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             >  1:2297346/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             <  1:2956140/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             >  2:1364445/1‑100 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCACTATTGCtgttg                                             >  1:2577086/1‑99 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                          <  2:1723826/99‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        >  1:1464887/1‑100 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        <  2:80032/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  1:93733/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  1:353804/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  2:617930/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGCGCAGTATTGCTGTTGTTGATTcc                                     <  1:94518/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  1:1859915/98‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                   >  1:423869/1‑100 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                   <  2:2758411/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                 >  1:1527508/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                 >  1:2018123/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGGAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                >  1:115003/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                              <  2:2577100/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  1:1320833/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             >  2:1832392/1‑100 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  1:109478/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            >  1:255234/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                           <  2:859033/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                           >  2:14329/1‑100 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGGAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                         <  2:115003/100‑1 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                        >  2:1004892/1‑100 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                        >  1:1953662/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                     >  1:2364000/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                     >  2:1251082/1‑100 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                   <  2:2120523/100‑1 (MQ=255)

                                                                               aGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                 <  2:1313919/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGGTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               >  1:2809004/1‑100 (MQ=255)

                                                                                  agaagaCTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             >  1:187602/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           <  1:525231/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           >  1:171947/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           <  1:1059075/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           >  2:970169/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt          >  1:805523/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTACTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         >  1:2152922/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        >  2:2320128/1‑100 (MQ=255)

                                                                                         tATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa       <  2:2741911/100‑1 (MQ=255)

                                                                                         tATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa       >  1:1562832/1‑100 (MQ=255)

                                                                                         tATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa       >  2:2294768/1‑100 (MQ=255)

                                                                                           tGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGCTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGGCGAAc     >  1:326918/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    >  2:395013/1‑100 (MQ=255)

                                                                                            gTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    <  1:1251078/100‑1 (MQ=255)

                                                                                             ttGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt   >  2:2048508/1‑100 (MQ=255)

                                                                                              tGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg  >  2:2974577/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280

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