BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	484,881	T→C	5.5%	I247V (ATC→GTC) ‡	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	484,881	0	T	C	5.5%	1142.0 / 35.3	329	I247V (ATC→GTC) ‡	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base T (152/158); new base C (13/5); total (165/164)
Fisher's exact test for biased strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 335 total aligned reads displayed.

TGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACC  >  NC_000913/484735‑485028

    cgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                               <  4:699167/149‑1 (MQ=255)

    cgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                               <  1:1855360/149‑1 (MQ=255)

         gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                                 <  1:904131/142‑1 (MQ=255)

         gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                                 <  1:404514/142‑1 (MQ=255)

         gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACtt                                                                                                                                          >  3:969043/1‑149 (MQ=255)

         gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACtt                                                                                                                                          <  4:1852153/149‑1 (MQ=255)

             aTACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGc                                                                                                                                      <  2:602535/149‑1 (MQ=255)

               aCCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCtt                                                                                                                                    <  1:1278022/149‑1 (MQ=255)

                 ccGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTg                                                                                                                                  <  3:868989/149‑1 (MQ=255)

                  cGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                                                                                 <  2:869798/149‑1 (MQ=255)

                   ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                                                                                <  2:1792186/149‑1 (MQ=255)

                       gcagAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtc                                                                                                                            <  4:1114739/149‑1 (MQ=255)

                          gagTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                               <  1:1049119/127‑1 (MQ=255)

                             tgtgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGt                                                                                                                      >  3:205377/1‑149 (MQ=255)

                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                                                                                <  1:1796312/135‑1 (MQ=255)

                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCa                                                                                                                  >  4:522540/1‑149 (MQ=255)

                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCa                                                                                                                  >  3:1772712/1‑149 (MQ=255)

                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCa                                                                                                                  >  1:614156/1‑149 (MQ=255)

                                      ggTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTg                                                                                                             <  3:522540/149‑1 (MQ=255)

                                           gggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                               >  3:1812044/1‑110 (MQ=255)

                                            ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTc                                                                                                       >  2:193561/1‑149 (MQ=255)

                                            ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAAAGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTc                                                                                                       >  1:1491710/1‑149 (MQ=255)

                                              aaGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGc                                                                                                     >  4:324049/1‑149 (MQ=255)

                                               aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCNTGCCATTCGcc                                                                                                    <  4:1772712/149‑1 (MQ=255)

                                                          aGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtct                                                                                                                           >  4:259067/1‑115 (MQ=255)

                                                          aGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtct                                                                                                                           <  3:259067/115‑1 (MQ=255)

                                                            ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTACCATTCGCCAGTTc                                                                                               >  3:1164248/1‑141 (MQ=255)

                                                             gTGATAGACCCAGTGGTCTGATCAACGGTAAAGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTctg                                                                                      <  2:1491710/149‑3 (MQ=255)

                                                                aTAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgcg                                                                                   >  2:1597785/1‑149 (MQ=255)

                                                                aTAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgcg                                                                                   >  3:1712768/1‑149 (MQ=255)

                                                                 tAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTc                                                                                        <  1:983689/143‑1 (MQ=255)

                                                                   gACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAg                                                                                <  3:1357183/149‑1 (MQ=255)

                                                                    aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCatg                                                                               >  2:1136042/1‑147 (MQ=255)

                                                                    aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCa                                                                                 <  2:1348427/147‑1 (MQ=255)

                                                                    aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCa                                                                                 >  3:217681/1‑147 (MQ=255)

                                                                    aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCa                                                                                 >  1:1348427/1‑147 (MQ=255)

                                                                      ccAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGtt                                                                                                                     >  3:34804/1‑109 (MQ=255)

                                                                      ccAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGtt                                                                                                                     <  4:34804/109‑1 (MQ=255)

                                                                        aGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtc                                                                                                                            >  2:399423/1‑100 (MQ=255)

                                                                          tGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCg                                                                         <  1:1136042/149‑1 (MQ=255)

                                                                           ggTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGt                                                                        <  3:367444/149‑1 (MQ=255)

                                                                             tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                                                                                 >  2:1484399/1‑90 (MQ=255)

                                                                             tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTc                                                                                        >  3:1479366/1‑131 (MQ=255)

                                                                               tGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGc                                                                                                     >  2:534868/1‑116 (MQ=255)

                                                                                gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAGCTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                   <  3:1629449/149‑1 (MQ=255)

                                                                                gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCa                                                                                       <  1:800214/129‑1 (MQ=255)

                                                                                gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                   <  2:42947/149‑1 (MQ=255)

                                                                                 aTCAACGGTAACGTCAGAGAATTCCTGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                   <  4:1588075/148‑1 (MQ=255)

                                                                                 aTCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                   >  3:1588075/1‑148 (MQ=255)

                                                                                 aTCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                   <  4:909774/148‑1 (MQ=255)

                                                                                  tCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGa                                                                 <  4:661901/149‑1 (MQ=255)

                                                                                    aaCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACt                                                               <  3:1418294/149‑1 (MQ=255)

                                                                                     aCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTc                                                                                        <  2:1234199/123‑1 (MQ=255)

                                                                                       ggTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg                                                            >  2:1160414/1‑149 (MQ=255)

                                                                                       ggTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGATCACCAGTGacg                                                                                                                                             <  4:379823/68‑2 (MQ=38)

                                                                                        gTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGt                                                           <  1:1781693/149‑1 (MQ=255)

                                                                                              tCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca                                                       <  3:463491/147‑1 (MQ=255)

                                                                                              tCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca                                                       >  4:463491/1‑147 (MQ=255)

                                                                                                agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                                                                      <  1:511640/130‑1 (MQ=255)

                                                                                                agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   >  4:820816/1‑149 (MQ=255)

                                                                                                agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   >  1:1295280/1‑149 (MQ=255)

                                                                                                agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   >  1:540284/1‑149 (MQ=255)

                                                                                                 gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg                                                            >  3:689793/1‑139 (MQ=255)

                                                                                                 gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAAc                                                  <  4:1326013/149‑1 (MQ=255)

                                                                                                  agaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACg                                                 >  1:1405621/1‑149 (MQ=255)

                                                                                                    aaTTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTa                                               <  4:1712768/149‑1 (MQ=255)

                                                                                                      ttCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGa                                             >  2:878325/1‑149 (MQ=255)

                                                                                                       tCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGAt                                            <  4:1663504/149‑1 (MQ=255)

                                                                                                        ccAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   <  1:870824/141‑1 (MQ=255)

                                                                                                        ccAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATc                                           >  2:1292745/1‑149 (MQ=255)

                                                                                                         cAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCg                                          >  4:212773/1‑149 (MQ=255)

                                                                                                           gCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGACCGGa                                        <  1:1597785/149‑1 (MQ=255)

                                                                                                              tACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAg                                                                                                  >  2:1718940/1‑88 (MQ=255)

                                                                                                               aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                    <  1:1144673/149‑1 (MQ=255)

                                                                                                               aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                    <  1:1444091/149‑1 (MQ=255)

                                                                                                               aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                    <  2:952858/149‑1 (MQ=255)

                                                                                                                  gTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGAt                                            >  4:745762/1‑138 (MQ=255)

                                                                                                                    ccTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTAGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtgc                               <  4:1373113/149‑1 (MQ=255)

                                                                                                                     cTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtgct                              <  2:1405621/149‑1 (MQ=255)

                                                                                                                        cGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca                                                       >  3:101364/1‑121 (MQ=255)

                                                                                                                        cGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca                                                       <  4:101364/121‑1 (MQ=255)

                                                                                                                         gggAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGAATCAAGTtgc                               >  1:1746588/1‑144 (MQ=255)

                                                                                                                             aCTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCCAGTTGCTGCACGGtt                      <  2:1852469/149‑1 (MQ=255)

                                                                                                                               ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAg                                                                            >  1:1716057/1‑93 (MQ=255)

                                                                                                                               ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                    >  4:1194957/1‑149 (MQ=255)

                                                                                                                               ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                    <  3:1194957/149‑1 (MQ=255)

                                                                                                                                       gTCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtg                                <  3:129378/129‑1 (MQ=255)

                                                                                                                                       gTCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCAGTTTCAGGCGCAGGAAGTCGCTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtg                                >  4:129378/1‑129 (MQ=255)

                                                                                                                                        tCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGt           >  2:1107926/1‑149 (MQ=255)

                                                                                                                                        tCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGt           <  2:782852/149‑1 (MQ=255)

                                                                                                                                        tCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAg                                              <  1:327907/114‑1 (MQ=255)

                                                                                                                                        tCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGt           >  3:785184/1‑149 (MQ=255)

                                                                                                                                         cACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCTCAGTc          <  2:1774914/149‑1 (MQ=255)

                                                                                                                                            tGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGcc       >  1:1229214/1‑149 (MQ=255)

                                                                                                                                             ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCt      <  3:324049/149‑1 (MQ=255)

                                                                                                                                                 aTCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATTAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGAcc  >  3:56155/1‑149 (MQ=255)

                                                                                                                                                 aTCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCa                           <  1:1473519/124‑1 (MQ=255)

                                                                                                                                                 aTCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGAcc  >  1:212956/1‑149 (MQ=255)

                                                                                                                                                  tCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   >  2:1444368/1‑99 (MQ=255)

                                                                                                                                                  tCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                   <  4:1457352/99‑1 (MQ=255)

TGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGACC  >  NC_000913/484735‑485028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`