BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,986,007	C→T	5.4%	D80D (GAC→GAT)	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,986,007	0	C	T	5.4%	331.0 / 6.7	112	D80D (GAC→GAT)	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand): ref base C (50/55); 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 = 3.72e-01

Only 100 of 113 total aligned reads displayed.

GCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGC  >  NC_002947/3985862‑3986154

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: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:864377/150‑1 (MQ=255)

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

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

                    gCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATc                                                                                                                             <  2:701399/150‑1 (MQ=255)

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

                      ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGcta                                                                                                                                    <  2:1103985/141‑3 (MQ=255)

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

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

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

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

                        aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAcc                                                                                                                         >  1:1845470/1‑150 (MQ=255)

                          gACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTGGCGGATCTCGGCCATGAAGTGCGTAGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGCTCTACcac                                                                                                                       >  1:771228/1‑150 (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                                                                                                                    <  1:193562/148‑1 (MQ=255)

                                  cAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAccc                                                                                                               >  1:599813/1‑150 (MQ=255)

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

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

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

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

                                            ttCGACCACTTGCCGGAAATGATCCTGCAGGGCCATGACATTCCCATCCCCCTGGGCTCGCGGATCCCGGCCATGAAGTGCGTGGCCGATGTCTTCGGCGACGACATGGCCGAAGGCGGCGTGATATACCACAACGTCGCCTCCTACAAg                                                                                                     >  1:198273/1‑150 (MQ=255)

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

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

                                               gACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGc                                                                                                  <  1:1265208/150‑1 (MQ=255)

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

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

                                                aCCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCa                                                                                                 <  2:1661627/150‑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)

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

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

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

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

                                                                gATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGt                                                                                                                                 <  1:2081734/102‑1 (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)

                                                                                         aTCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTAc                                                        >  1:167158/1‑150 (MQ=255)

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

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

                                                                                            cACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCGCGACTGCTGCATGTACAAACCGGTGTATTACCAg                                                     <  2:2300138/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                  ggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                                               <  2:1937562/134‑1 (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)

                                                                                                           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                                                                                               >  1:217099/1‑85 (MQ=255)

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

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

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

                                                                                                                           gCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCACgg                      <  1:2406558/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                  cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAc                                                                                   <  2:412852/82‑1 (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)

GCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGC  >  NC_002947/3985862‑3986154

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