BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,601,652	T→A	5.6%	L147Q (CTG→CAG)	osmY →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,601,652	0	T	A	5.6%	709.0 / 24.8	232	L147Q (CTG→CAG)	osmY	periplasmic protein
Reads supporting (aligned to +/- strand): ref base T (105/114); new base A (6/7); total (111/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 242 total aligned reads displayed.

AGGCCGAAGAGGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGTTAAAA  >  CP009273/4601505‑4601798

aggccGATGAGGCAGTGAAAGTGGCTAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGATACCGCCACCACCAGTTAAATCAAAGCCAAActgct                                                                                                                                                    <  2:1243660/148‑1 (MQ=255)

      aaGAGGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGa                                                                                                                                              <  2:905996/148‑1 (MQ=255)

      aaGAGGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCGGGCAGCGACAAACGGCACGTTCGCGACGCTAAAGACGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGa                                                                                                                                              <  2:305722/148‑1 (MQ=255)

        gagGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTTAAATCAAAGCCAAACTGCTGGCGGACg                                                                                                                                            <  1:2161495/148‑1 (MQ=255)

        gagGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACg                                                                                                                                            >  1:1858238/1‑148 (MQ=255)

           gCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGata                                                                                                                                         <  2:2519874/148‑1 (MQ=255)

           gCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGata                                                                                                                                         >  2:2116598/1‑148 (MQ=255)

           gCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGata                                                                                                                                         <  1:1927115/148‑1 (MQ=255)

            cAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGatat                                                                                                                                        >  2:1605821/1‑148 (MQ=255)

               tGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGt                                                                                                                                     >  1:3481226/1‑148 (MQ=255)

                gtaaGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGGGACACCGCCACCACCAGTGAAATCTAAGCCAAACTGCAGGCGGACGATATCGTc                                                                                                                                    <  2:2452369/146‑1 (MQ=255)

                 aaaGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTcc                                                                                                                                   >  1:766859/1‑148 (MQ=255)

                 aaaGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTcc                                                                                                                                   <  2:2257447/148‑1 (MQ=255)

                 aaaGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTcc                                                                                                                                   >  2:1915329/1‑148 (MQ=255)

                  taggtgCGAAAGGCGTTGAAGGGGTTACCTCTGTCAGAGAAAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTTAAGGGCTACGCGGGTGACACCGAAACCACCACTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTccg                                                                                                                                  <  2:2120625/143‑2 (MQ=255)

                  aaGGGGCGAAAGGCTTTGACGGGGGGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCGCGGTGACGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTccc                                                                                                                                  <  1:949909/148‑1 (MQ=255)

                    gTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCtt                                                                                                                                <  1:441926/148‑1 (MQ=255)

                     tGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTc                                                                                                                               <  2:3405615/148‑1 (MQ=255)

                      ggCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTcc                                                                                                                              >  1:879150/1‑148 (MQ=255)

                       gagAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTccc                                                                                                                             >  2:730900/3‑148 (MQ=255)

                          aaaGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTc                                                                                                                          >  1:2884645/1‑148 (MQ=255)

                           aaGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCa                                                                                                                         <  1:2532535/148‑1 (MQ=255)

                             ggCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCAtg                                                                                                                       <  1:2187304/148‑1 (MQ=255)

                              gCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCAtgt                                                                                                                      >  2:2900133/1‑148 (MQ=255)

                                  tGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGaaa                                                                                                                  >  2:2023300/1‑148 (MQ=255)

                                   gAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAg                                                                                                                 <  1:1308137/148‑1 (MQ=255)

                                       gggTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGa                                                                                                             >  2:819598/1‑148 (MQ=255)

                                       gggTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCAGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGa                                                                                                             <  2:2434061/148‑1 (MQ=255)

                                         gTGACCTCTGTCAGCGAAAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAACGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGaaa                                                                                                           <  1:1825928/148‑1 (MQ=255)

                                            aCCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATTTGAAAGTTGAAacca                                                                                                        >  1:2411086/1‑148 (MQ=255)

                                            aCCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAacca                                                                                                        >  2:1034340/1‑148 (MQ=255)

                                              ctctGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAaccacc                                                                                                      >  2:968334/1‑148 (MQ=255)

                                              ctctGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAaccacc                                                                                                      >  2:3473989/1‑148 (MQ=255)

                                              ctctGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAaccacc                                                                                                      <  1:2612445/148‑1 (MQ=255)

                                                     aGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCg                                                                                               <  2:436134/148‑1 (MQ=255)

                                                      gCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGt                                                                                              <  2:2109667/148‑1 (MQ=255)

                                                      gCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGt                                                                                              >  1:1016552/1‑148 (MQ=255)

                                                       cGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCAGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTg                                                                                             >  2:2659349/1‑148 (MQ=255)

                                                         aCAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGt                                                                                           >  1:2992949/1‑148 (MQ=255)

                                                          cAAACTGCACGTTCGCGACGCTAAAGACGGCTGGGTGATGGGCTACGGGGGTGACACCGCCACCACCAGTGAAATCAACGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGtt                                                                                          <  1:672226/148‑1 (MQ=255)

                                                          cAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGtt                                                                                          <  1:1915329/148‑1 (MQ=255)

                                                             aCTGCTCGTTCGCGACTCTATAGTAGGCTCGGTGTAGGTCTCCGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCGAAACTGGTGGCTGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACGGACGGCGTGGttaat                                                                                       <  2:2925425/148‑4 (MQ=255)

                                                             aCTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAg                                                                                       >  1:3381496/1‑148 (MQ=255)

                                                              cTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGc                                                                                      >  1:1217264/1‑148 (MQ=255)

                                                                 cACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGctct                                                                                   <  2:1697395/148‑1 (MQ=255)

                                                                 cACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGctct                                                                                   >  2:2705687/1‑148 (MQ=255)

                                                                     ttCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCgg                                                                               >  1:940426/1‑148 (MQ=255)

                                                                     ttCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCAGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCgg                                                                               <  2:575580/148‑1 (MQ=255)

                                                                              cTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGa                                                                      <  1:3552566/148‑1 (MQ=255)

                                                                              cTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGa                                                                      <  1:1490934/148‑1 (MQ=255)

                                                                               tAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGag                                                                     >  1:2420416/1‑147 (MQ=255)

                                                                                  agaagGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATtct                                                                  >  1:3267852/1‑148 (MQ=255)

                                                                                   gaagGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATtctc                                                                 >  2:1491469/1‑148 (MQ=255)

                                                                                   gaagGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATtctc                                                                 >  2:1870175/1‑148 (MQ=255)

                                                                                    aagGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCa                                                                >  1:1274294/1‑148 (MQ=255)

                                                                                      ggCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAgg                                                              >  1:578428/1‑148 (MQ=255)

                                                                                       gCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCGAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGc                                                             >  1:1557520/1‑148 (MQ=255)

                                                                                       gCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGc                                                             >  1:188351/1‑148 (MQ=255)

                                                                                        cTCGGTGATGGGATACGCGGGTGACACCGCCACCACCAGTTAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGca                                                            <  1:2405296/148‑1 (MQ=255)

                                                                                         tCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGcac                                                           <  1:2837964/148‑1 (MQ=255)

                                                                                            gTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACaaa                                                        <  2:1352304/148‑1 (MQ=255)

                                                                                             tGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAg                                                       >  1:2782331/1‑148 (MQ=255)

                                                                                              gAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGt                                                      >  2:404964/1‑148 (MQ=255)

                                                                                                aGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGa                                                    >  2:3451266/1‑148 (MQ=255)

                                                                                                aGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGa                                                    <  2:523072/148‑1 (MQ=255)

                                                                                                aGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGa                                                    >  1:3074881/1‑148 (MQ=255)

                                                                                                 gggCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGAc                                                   <  1:1138512/148‑1 (MQ=255)

                                                                                                  ggCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGAcc                                                  >  1:1604498/1‑148 (MQ=255)

                                                                                                  ggCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCAGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGAGTCTCAGGCACACAGTGAcc                                                  >  2:2968436/1‑148 (MQ=255)

                                                                                                    cTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGt                                                <  2:1316313/148‑1 (MQ=255)

                                                                                                     tACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTg                                               <  1:1154512/148‑1 (MQ=255)

                                                                                                      aCGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGc                                              <  2:1961397/148‑1 (MQ=255)

                                                                                                            gTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAg                                        >  2:2751469/1‑148 (MQ=255)

                                                                                                             tGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGt                                       <  1:1232200/148‑1 (MQ=255)

                                                                                                             tGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCGGCAAGt                                       >  1:1007329/1‑148 (MQ=255)

                                                                                                             tGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCAGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGt                                       >  2:1203201/1‑148 (MQ=255)

                                                                                                               acacCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTAt                                     <  2:2068319/148‑1 (MQ=255)

                                                                                                               acacCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTAt                                     >  1:3422628/1‑148 (MQ=255)

                                                                                                                   cGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTTCTGAAAGTATCGcc                                 >  1:1075721/1‑148 (MQ=255)

                                                                                                                       accaccaGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAg                             <  2:2290330/148‑1 (MQ=255)

                                                                                                                       accaccaGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAg                             >  2:3186229/1‑148 (MQ=255)

                                                                                                                       accaccaGTGAAATCAAAGCCAAACTGCTGGCGGACGAGATCGTCCCGTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAg                             <  2:498019/148‑1 (MQ=255)

                                                                                                                        ccaccaGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGc                            <  2:3218027/148‑1 (MQ=255)

                                                                                                                         caccaGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCg                           <  1:729763/148‑1 (MQ=255)

                                                                                                                         caccaGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCg                           <  2:2623068/148‑1 (MQ=255)

                                                                                                                          accaGTGAAAACAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCgg                          <  1:2698455/148‑1 (MQ=255)

                                                                                                                               tGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGAt                     <  1:1001394/148‑1 (MQ=255)

                                                                                                                                 aaaTCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATgg                   <  1:2017921/148‑1 (MQ=255)

                                                                                                                                 aaaTCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATgg                   >  2:67691/1‑148 (MQ=255)

                                                                                                                                 aaaTCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATgg                   <  1:90899/148‑1 (MQ=255)

                                                                                                                                   aTCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCCCAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGgtg                 >  1:1444391/1‑148 (MQ=255)

                                                                                                                                    tCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGgtgt                <  2:1613702/148‑1 (MQ=255)

                                                                                                                                    tCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGgtgt                <  1:1694692/148‑1 (MQ=255)

                                                                                                                                        aGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGaaa            >  2:158851/1‑148 (MQ=255)

                                                                                                                                            aaaCTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCg        <  2:1153489/148‑1 (MQ=255)

                                                                                                                                             aaCTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGt       >  1:3126802/1‑148 (MQ=255)

                                                                                                                                              aCTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGCCCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGtt      >  2:1482319/1‑148 (MQ=255)

                                                                                                                                              aCTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTTAAAAGCGtt      >  1:1985597/1‑148 (MQ=255)

                                                                                                                                               ctgctgGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTAGATTCTAAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGTTa     >  2:3238786/1‑148 (MQ=255)

                                                                                                                                                  ctgGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGTTaaaa  <  2:940314/148‑1 (MQ=255)

AGGCCGAAGAGGCAGTGAAAGTGGCGAAAGGCGTTGAAGGGGTGACCTCTGTCAGCGACAAACTGCACGTTCGCGACGCTAAAGAAGGCTCGGTGAAGGGCTACGCGGGTGACACCGCCACCACCAGTGAAATCAAAGCCAAACTGCTGGCGGACGATATCGTCCCTTCCCGTCATGTGAAAGTTGAAACCACCGACGGCGTGGTTCAGCTCTCCGGTACCGTCGATTCTCAGGCACAAAGTGACCGTGCTGAAAGTATCGCCAAAGCGGTAGATGGTGTGAAAAGCGTTAAAA  >  CP009273/4601505‑4601798

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