BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,986,009	A→T	5.4%	D81V (GAC→GTC)	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,986,009	0	A	T	5.4%	353.5 / 8.6	112	D81V (GAC→GTC)	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand): ref base A (50/56); new base T (3/3); total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 118 total aligned reads displayed.

GCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGCC  >  NC_002947/3985862‑3986158

gcTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACa                                                                                                                                                     >  2:389961/1‑150 (MQ=255)

 cTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCTAGGTCTTCGGCGACGACAt                                                                                                                                                    >  2:1933727/1‑150 (MQ=255)

  tCCAGTGTGTTCAACCTCGCCCTCGACTATTCCAACTCGCTGTTCGACCACTTTCCGGAAATGATCTTGCAGGGCTAGGATATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGTCGACGACATc                                                                                                                                                   >  2:353088/1‑150 (MQ=255)

   ccAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCg                                                                                                                                                  <  1:1171996/150‑1 (MQ=255)

    cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGc                                                                                                                                                 >  1:1661627/1‑150 (MQ=255)

    cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGc                                                                                                                                                 >  1:2089220/1‑150 (MQ=255)

       tgtgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGa                                                                                                                                              <  2:174989/150‑1 (MQ=255)

        gtgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                             <  2:1334289/150‑1 (MQ=255)

        gtgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGAAGACATCGcctta                                                                                                                                             <  2:1200913/150‑4 (MQ=255)

         tgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGAAg                                                                                                                                            <  2:818472/150‑1 (MQ=255)

          gtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAgg                                                                                                                                           >  1:622059/1‑150 (MQ=255)

           ttaagacTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGc                                                                                                                                          >  2:1899822/7‑150 (MQ=255)

                ccTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGt                                                                                                                                     >  1:1615813/1‑150 (MQ=255)

                ccTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGt                                                                                                                                     <  1:864377/150‑1 (MQ=255)

                 cTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTg                                                                                                                                    >  1:247927/1‑150 (MQ=255)

                 cTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTg                                                                                                                                    <  1:1899822/150‑1 (MQ=255)

                      ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGcta                                                                                                                                        >  1:1103985/1‑139 (MQ=255)

                      ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTa                                                                                                                               >  1:2119866/1‑150 (MQ=255)

                      ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTa                                                                                                                               >  2:1620831/1‑150 (MQ=255)

                       taCGACTATTCCAACTCGCTGTTCCACTACTTGCCGGAATTGATCCTGCGGGGCCAGGACGTTCGCATCCTCGTGTGCTCGCTGTTCTCGTGCATGAGTTGCGTTTCCGAGTTCTTCGGCGGCTAGAGCGCGGTGGGCGGCGTGATCttc                                                                                                                              <  2:247927/149‑3 (MQ=255)

                       cACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAc                                                                                                                              <  1:1480197/150‑1 (MQ=255)

                          gACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACcac                                                                                                                           >  2:1912360/1‑150 (MQ=255)

                             tATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        <  1:389961/150‑1 (MQ=255)

                             tATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        >  1:1392497/1‑150 (MQ=255)

                               ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGGAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGa                                                                                                                      >  1:2300138/1‑150 (MQ=255)

                               ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        >  2:193562/1‑148 (MQ=255)

                               ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        <  1:193562/148‑1 (MQ=255)

                               ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGa                                                                                                                      >  2:169979/1‑150 (MQ=255)

                                   aaCTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAcccc                                                                                                                  >  1:1677844/1‑150 (MQ=255)

                                     cTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGc                                                                                                                <  1:1912360/150‑1 (MQ=255)

                                        gCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATGGCCGAAGGCGACGTGATCTACCACAAAGACACAGCCTa                                                                                                             >  1:111089/1‑150 (MQ=255)

                                           gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGa                                                                                                                                   >  2:571627/1‑125 (MQ=255)

                                           gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGa                                                                                                                                   <  1:571627/125‑1 (MQ=255)

                                            ttCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAg                                                                                                         <  1:2195567/150‑1 (MQ=255)

                                            ttCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAg                                                                                                         <  1:907035/150‑1 (MQ=255)

                                                aCCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAgg                                                                                                        >  1:1314633/1‑147 (MQ=255)

                                                aCCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAgg                                                                                                        <  2:1314633/147‑1 (MQ=255)

                                                  cACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGc                                                                                                   <  2:1845470/150‑1 (MQ=255)

                                                     ttGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        <  2:1690510/126‑1 (MQ=255)

                                                     ttGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                        >  1:1690510/1‑126 (MQ=255)

                                                         cGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtc                                                                                            <  1:596741/150‑1 (MQ=255)

                                                          ggAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtct                                                                                           >  2:2314439/1‑150 (MQ=255)

                                                           gAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtctc                                                                                          <  1:709268/150‑1 (MQ=255)

                                                             aaTGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGa                                                                                        >  2:428550/1‑150 (MQ=255)

                                                             aaTGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGa                                                                                        <  1:996726/150‑1 (MQ=255)

                                                              aTGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAc                                                                                       <  1:2314439/150‑1 (MQ=255)

                                                               tGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAct                                                                                      <  1:169979/150‑1 (MQ=255)

                                                                gATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGt                                                                                                                                     <  1:2081734/102‑1 (MQ=255)

                                                                gATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGt                                                                                                                                     >  2:2081734/1‑102 (MQ=255)

                                                                gATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGActg                                                                                     >  2:743316/1‑150 (MQ=255)

                                                                 aTCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGCCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGActgc                                                                                    <  2:1615813/150‑1 (MQ=255)

                                                                    cTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACtgctgc                                                                                 >  2:785140/1‑150 (MQ=255)

                                                                    cTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACtgctgc                                                                                 >  2:2064759/1‑150 (MQ=255)

                                                                         gggCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTa                                                                            <  1:743316/150‑1 (MQ=255)

                                                                             cAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACaaa                                                                        >  2:975210/1‑150 (MQ=255)

                                                                                gACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCg                                                                     >  1:1761193/1‑150 (MQ=255)

                                                                                 aCATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCgg                                                                    <  2:1677844/150‑1 (MQ=255)

                                                                                  cATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   <  1:1066287/150‑1 (MQ=255)

                                                                                  cATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   <  2:2119866/150‑1 (MQ=255)

                                                                                   aTTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtg                                                                  <  2:1164368/150‑1 (MQ=255)

                                                                                   aTTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtg                                                                  >  2:1017801/1‑150 (MQ=255)

                                                                                    ttCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGACAGCCATATTCTCGACTGCTGCATGTACAAACCGgtgt                                                                 <  2:1392497/150‑1 (MQ=255)

                                                                                       ccATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGGATTCTCCTACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTAtt                                                              <  1:1620831/150‑1 (MQ=255)

                                                                                           ccACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   <  1:1409145/141‑1 (MQ=255)

                                                                                            cACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAg                                                         <  1:1915128/150‑1 (MQ=255)

                                                                                             aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtctc                                                                                          >  2:157276/1‑116 (MQ=255)

                                                                                             aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAgg                                                        >  1:2160792/1‑150 (MQ=255)

                                                                                                 gggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                             <  2:718326/61‑1 (MQ=37)

                                                                                                 gggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                             >  1:718326/1‑61 (MQ=37)

                                                                                                 gggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGa                                                    >  1:495381/1‑150 (MQ=255)

                                                                                                  ggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   >  1:1937562/1‑134 (MQ=255)

                                                                                                   gCTCGCTTATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGc                                                  <  1:1155255/150‑1 (MQ=255)

                                                                                                      cGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTgg                                               <  1:785140/150‑1 (MQ=255)

                                                                                                           tccccTGCCATGAAGTGCGTTGCCGACTTCGTCGGCGACGACAGCGTGAGAGGCGGCTGTATTGACCACGGCTACCCCGCCGCCAATGGCAGCCATGTTCGCGACTGCTGGATGTATAAGCCGGTGTATTACCAGGACGACCTGCTGTTc                                          <  2:111089/148‑1 (MQ=255)

                                                                                                           aTCCCGGCCATGAAGTGCGTTGCCGAGTTCTTCGGCGACGACATCGTCTACTGTTACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTc                                          <  1:1352731/150‑1 (MQ=255)

                                                                                                           aTCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtctc                                                                                          >  1:679845/1‑102 (MQ=255)

                                                                                                           aTCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtctc                                                                                          <  2:679845/102‑1 (MQ=255)

                                                                                                            tCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCt                                         <  1:975210/150‑1 (MQ=255)

                                                                                                             cccGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTg                                        <  1:2064759/150‑1 (MQ=255)

                                                                                                              ccGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTgg                                       >  2:2099740/1‑150 (MQ=255)

                                                                                                                   cATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGc                                                                                                   <  2:217099/85‑1 (MQ=255)

                                                                                                                    aTGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGtt                                 >  1:2241404/1‑150 (MQ=255)

                                                                                                                     tGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGttt                                <  2:622059/150‑1 (MQ=255)

                                                                                                                           gCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAgg                          <  2:167158/150‑1 (MQ=255)

                                                                                                                                gCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCAc                     <  1:2112431/150‑1 (MQ=255)

                                                                                                                                 ccGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   >  2:1638700/1‑103 (MQ=255)

                                                                                                                                 ccGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                                   <  1:1638700/103‑1 (MQ=255)

                                                                                                                                  cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAc                                                                                       >  1:412852/1‑82 (MQ=255)

                                                                                                                                  cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGTCGACCt                   <  2:2089220/150‑1 (MQ=255)

                                                                                                                                  cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCt                   >  1:1655205/1‑150 (MQ=255)

                                                                                                                                         ttcGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGAc            >  2:415357/1‑150 (MQ=255)

                                                                                                                                         ttcGGCGATGTCGTCGCCAAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGAc            >  2:1898274/1‑150 (MQ=255)

                                                                                                                                         ttcGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGAc            >  1:1122562/1‑150 (MQ=255)

                                                                                                                                            ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATc         <  2:599813/150‑1 (MQ=255)

                                                                                                                                               gacgacATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATcggc      >  1:1684028/1‑150 (MQ=255)

                                                                                                                                                  gacATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCaa                            >  2:2312044/1‑125 (MQ=255)

                                                                                                                                                  gacATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCaa                            <  1:2312044/125‑1 (MQ=255)

                                                                                                                                                   acATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGa                                                    >  1:1561063/1‑100 (MQ=255)

                                                                                                                                                   acATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGa                                                    <  2:1561063/100‑1 (MQ=255)

                                                                                                                                                   acATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAAACGGGGTATTACTAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGcc  >  1:171795/1‑150 (MQ=255)

GCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGCC  >  NC_002947/3985862‑3986158

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